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INTRODUCTION 


To determine the complete convective heat-transfer distribution 
for configurations flying at hypersonic speeds, one must consider the 
viscous :inviscid interactions associated with the complex three-dimen- 
sional flow fields. Because of the complexity of the viscous: inviscid 
interaction phenomena, many investigators have studied the locally per- 
turbed flow fields using models consisting of basic elemental combinations. 
By varying the relative position of shock-generating elements, Edney^ ob- 
served that there are six basic shock-interference patterns. Rains and 
Keyes^ have categorized shock-interference patterns obtained for a 
variety of space-shuttle configurations in terms of the models by Edney. 

O 

Bertin, et al. , examined surface-pressure and heat-trar sf er-rate 
data for a variety of shuttle-orbiter configurations over an angle-of- 
attack range from 0^ to 60 • It was found that the *typB of shock-inter- 
ference pattern was dominated by the leading-edge effective sweep angle. 

For the relatively low sweep angles of the straight-wing orbiters, the 
interaction between the bow— generated shock— wave and the wing— generated 
shock— wave was a **Type V” shock— interference pattern. For the delta— wing 
orbiters, the shock: shock interaction exhibited the characteristics of a 
”Type VI" pattern for all alphas. 

Theoretical calculations^’^ of the inviscid shock-interaction flow- 
field, coupled with empirically derived surface-pressure :heat-transfer 
correlations of the viscous interaction phenomena, indicate that the local 
increases in heat transfer and in surface pressure associated with the 
Type IV and with the Type V interactions become markedly more severe as 
Y decreases. Since the current shuttle entry-configurations have highly 
swept wings, the present investigation was undertaken to determine the 
effect of the gas-property model on the Type VI shock-interference pattern. 
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Numerical codes were developed to calculate the two-dimensional flow field 
which results when supersonic flow encounters double wedge configurations 
whose angles are such that a Type VI pattern occurs. The flow-field 
model included the shock-interaction phenomena observed in ref. 3 for a 
delta-wing orbiter. Two numerical codes were developed: one which used the 
perfect gas relations and a second which incorporated a Mollier table to 
define equilibrium air properties. The two codes were used to generate 
theoretical surface-pressure and heat-transfer distributions for velo- 
cities from 1167 m/sec (3821 ft/sec), i.e., a wind-tunnel condition, to 
7610 m/sec (25,000 ft/sec), i.e., an entry condition. 



NOMENCLATURE 


a 

s 

dn 

f" 

g' 

h 

H 

M 

n ” 

P 

Pr 

q 

,ref 

r 

R 

Re 

X 

R 

gas 

5 

T 

T 

r 

U 

X 

a ~ 

Y 

6 
n 


speed of sound 
pressure coefficient 

pressure gradient normal to surface 
velocity gradient function 
stagnation enthalpy gradient function 
static enthalpy 
stagnation enthalpy 
Mach number 

transformed p-coordinate 
pressure 
Prandtl number 
heat-transfer rate 

heat-transfer rate to reference sphere 

radius of cross-section for a body of revolution 

radius of reference sphere 

Reynolds number based on local properties and wetted length 
gas constant 

wetted distance along wing leading-edge, also entropy 
static temperature 
recovery temperature 

streamwise velocity (capital letter denotes free-stream velocity) 

axial coordinate 

transformation parameter 

ratio of specific heats 

initial turning angle 

transformed y-coordinate 
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0 - shock angle 

A - sweep angle (see Fig. 1) 

- complement of sweep angle (see Fig. 1) 
y ~ viscosity 

V - Prandtl-Meyer angle 

£ - pressure ratio across a shock wave 

p - density 

Subscript 

e ~ parameter evaluated at the boundary layer edge 

ex - parameter evaluated in the expansion fan 

" “ parameter evaluated in the free-stream region 

n - parameter normal to shock wave 

t - indicates stagnation condition 

t2 - indicates stagnation condition in region 2 

w - parameter evaluated at the wall 



THEORETICAL ANALYSIS 


The complex flow field, which is established when a high-speed 
flow encounters a double-wedge configuration, is dominated by a 
shock-interaction region which imposes a highly non-uniform flow field 
adjacent to the wedge boundary layer. If the two wedge angles are 
not too large, the shock waves associated with the flow deflection 
are attached. The shock- interference pattern which results for this 
case has been designated by Edney as a Type VI pattern. For the 
computer code developed in the present study, it was assumed that 
the Type VI shock- interference pattern for a double wedge configuration 
(as shown in Fig. 1) includes: 

1) the undisturbed free-stream flow, 

2) the flow turned through the angle 6 by a single shock wave, 

3) the flow turned through the angle by two shock waves, 

4) the flow processed by the right-running waves of the expansion 
fan which are centered at the intersection of the two shock 
waves , 

5) the flow which passes through the left-running waves produced 
by the reflection of the waves of the expansion fan, and 

6) the flow turned through the angle A^ by a single shock wave. 
The flow near the juncture of the two wedges, i.e., in region 3, has 
passed through two shock waves. Further outboard on the downstream 
wedge, i.e., in region 6, the flow has passed through a single 
shock wave and is, therefore, at a lower pressure than that in region 
3. Thus, although the flow directions are the same, the gas must 
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undergo an expansion from the root region to equalize the pressure. 

The flow accelerates isentropically through the expansion regions 
so that the pressure and the flow direction in region 5E are identical 
to the values for region 6. The current study is concerned with the 
shock-interaction flow-field phenomena, which directly influence the 
wedge surface flow properties (i.e., pressure, shear, and heat trans- 
fer). Therefore, the flow-field properties along the wall were cal- 
culated including the interactions between the right-running and left- 
running waves in the numerical code. However, no attempt was made to 
model the shear layer which develops between the shock wave and the 
expanding flow, since this shear layer does not interact with the 
wedge surface for the two-dimensional Type VI interaction. 

Perfect-gas code. - The first steps were the calculation of the flow 
conditions downstream of an oblique shock for a given flow-deflection 
angle and for a given gamma, i.e., in regions 2, 3, and 6. The 
flow conditions in region 5E could then be calculated, since the 
static pressure in region 5E is equal to the static pressure in 
region 6 and, under the isentropic-expansion assumption, the stagnation 
pressure in region 5E is equal to the stagnation pressure in region 
3. The expansion process by which the gas accelerates from region 3 
to region 5 e was divided into ten equal steps. To satisfy the physical 
boundary condition that the flow in region 5E be parallel to the wall, 
the total change in the Prandtl-Meyer angle was divided into ten 
equal parts: the five right-running waves of region 4 and the five 

reflected waves constituting region 5. The governing equations for 
these calculations are summarized in ref. 6. Having defined the 
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inviscid flow-field and, hence, the conditions at the edge of the 
boundary layer, the heat- transfer distribution along the downstream 
wedge was calculated using the Eckert-reference-temperature technique^* 

For the wind-tunnel flow-condition, where the perfect gas relations 
accurately describe the gas behavior, the Eckert-reference-temperature 
heating rates compared favorably with the values obtained using the 
nonsimilar boundary-layer code described below* Due to the simplicity 
of the Eckert method, it was used with the perfect-gas calculations. 

The boundary-layer of the downstream wedge was assumed to originate at 
the intersection of the two wedges (point 0 in Fig. 1). 

One way of approximating the high- temperature, or real-gas, pro- 
perties of air is to use lower values of gamma in the perfect-gas rela- 
tions. Therefore, the equations for the perfect-gas code have been 
written so that one can input one value of gamma for regions 1 and 2, 
another gamma for regions 3 through 5E, and a third value for region 6. 

Thus, one can "account for" the varying shock strengths. The perfect- 
gas code was used to generate three "different" types of solution. For 
the present paper, these types of solution are referred to as: 

(a) "perfect-gas" solution for which y = 1.400 throughout the 
flow field, 

(b) "constant-gamma" solution for which y = 1.200 thoughout the 
flow field, 

(c) "variable-gamma" solution for which different values of gamma 
were assigned to the three input gamma parameters. (The 
required values were obtained from the real-gas solutions.) 

To be consistent in the perfect-gas assumption, the specific heat 
of air was held constant both for the constant-gamma and the variable- 
gamma solution. Sutherland’s relation^ was used to calculate the viscosity. 
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Assuming the Prandtl number to be 0.7 uniquely determined the thermal 
conductivity. 

Real-gas code . - Philosophically, the calculation procedure for the 
shock-interaction pattern using the real-gas code was similar to that 
described for the perfect-gas code. However, to account for the high- 
temperature, or "real-gas”, effects all thermodynamic properties are 

g 

evaluated using numerical charts for air in chemical equilibrium . 

Reference 9 was used to define the temperature dependence of the 
transport properties of equilibrium air, i.e., viscosity, thermal con- 
ductivity, and specific heat. Once the oblique shock relations have 
been used to define the static pressure and the entropy for the two end 
regions, the expansion from region 3 to region 5E is divided into ten 
equal steps. The flow conditions in the intermediate regions of the 
isentropic expansion are calculated using the relations of ref. 10. 

The relations require the static enthalpy and the local speed of sound 
for the intermediate regions, which are evaluated using the tables for 
the equilibrium air^. 

The heat-transfer distribution along the downstream wedge is cal- 
culated using a nonsimilar boundary-layer code^^, which was modified so 
that the thermodynamic properties for the viscous flow would also be 
calculated using ref. 8. In addition to accounting for the "real-gas" 
effects, the effect of the acceleration of the inviscid flow is included 
in the nonsimilar code. As before, the boundary layer is assumed to 
originate at the junction of the two wedges. 



DISCUSSION OF RESULTS 


The objective of the current analytical investigation was to 
determine how to apply the shock-interference data obtained in the 

3 

wind tunnel to the reentry of a shuttle orb iter. For simplicity, 
the numerical code was restricted to a two-dimensional flow model with 
the second, or downstream, wedge representing the wing leading edge. 

The initial deflection angle 6 was 5° for all solutions. This value 
was chosen, because the shock standoff distance in the vicinity of 
the wing leading-edge correlated reasonably well with the bow-wave 

trace observed in wind tunnel tests. The deflection angle A of 

s 

the downstream wedge was varied from 25^ to 63°, which corresponds 

to leading-edge sweep angles from 65^ to 27°. The dimensions of the 
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0.009 scale orbit er tested in Tunnel B of AEDC were used to define 
the characteristic lengths of the wedges. 

For this study, numerical solutions were generated to determine 
the effect of gas properties on the flow field and of the wall temperature 
on the heat transfer in the region where the Type VI shock-interaction 
influenced the "wing leading-edge". Flox^-field solutions were 
generated for three free-stream conditions. 

1) a wind tunnel condition where 

U = 1167 m/sec, P = 2.98 mmHg, T = 53°K 
(U = 3821 ft/sec, P = 0.057 psia, T = 95°R) 

2) an orbiter entry condition where 

U = 4330 m/sec, P =0.333 mmHg, T = 273°K 
(U^ = 14,200 ft/sec, P^ = 0.0064 psia, T^ = 491°R) 


9 
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3 ) an orbiter entry condition where 

U = 7610 m/sec, P = 0.0268 mmHg, T = 195°K 

oo '00 oo 

(U = 25,000 ft/sec, P = 0.00052 psia, T = 352°R) . 

Solutions were obtained using both the perfect-gas code and the real-gas 

code at all flow conditions. In addition, the variable-gamma option was 

used to generate solutions for flow conditions 2 and 3. Heat-transfer 

distributions along the downstream wedge were obtained for all three free- 

stream conditions for a wall temperature T of 394 K (710 R) and for the 

w 

two entry conditions for a of 1640 K (2960 R) . 

The effect of the gas properties on the calculated geometry of the 
Type VI shock-interference pattern is illustrated in Fig. 2. The output 
from the real-gas solutions was used to define the gamma distribution 
for input for the variable-gamma solution. Thus, = Y 2 ~ 1.400, 

= Y 3 = 1.214, while = 1.163. The geometry for the variable- 
gamma solution compares favorably with the real-gas geometry. The region 
where the leading edge is influenced by the shock interaction is essen- 
tially the same for these two solutions. The shock layer is thicker for 
the real-gas solution. For these deflection angles, the location and 
the extent of the surface affected by the expansion fan differs little 
between the perfect-gas solution and the real-gas solution. This sim- 
ilarity exists even though, in region 3, y is 1.214 for the real-gas 
solution and is, of course, 1.4000 for the perfect-gas solution. How- 
ever, because the density ratio across a shock wave is higher when the 
real-gas properties are accounted for , the shock layer along the wing 
leading-edge is markedly thinner for the real-gas solution. 

The pressure distributions along the leading edge of a 60° sweep 
"wing" are presented in Fig. 3 for the entry velocity of 4330 m/sec. 

So that the pressure variations could be seen more clearly, the scale 
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has been greatly expanded for the region where the expansion fan impinges 
on the "wing leading-edge". The vertical marks in Fig. 3 indicate the 
various regions of the field, using the same legend as the theoretical 
distributions. 

As noted in Fig. 2, the interaction-perturbed region for the perfect 
gas solution is inboard relative to the real-gas solution, although the 
locations differ only slightly. However, for a given region, the perfect 
gas surface-pressure correlates quite well with the real-gas value (Fig. 
3a), with the perfect-gas solution yielding a slightly higher pressure. 
Over the range of the free-stream conditions of the present study, the 
perfect-gas surface-pressure in a given region of the flow field was 
within 10% of the real-gas value. As noted in Fig. 2, the locations of 
the interaction-perturbed regions for the real-gas solution and for the 
variable-gamma solution are in close agreement. However, the pressures 
for the variable-gamma solution are significantly lower than the real- 
gas solution. The discrepancy between the real-gas and the variable- 
gamma values of surface pressure in a given region typically varied from 
15% to 25%. 

The heat-transfer distributions along the wing leading-edge are 

presented in Fig. 4 for the entry velocity of 4330 m/sec. The 

local heat-transfer rate has been divided by the stagnation-point 
12 

heating rate for a reference sphere, whose radius was chosen to be 
0.0027m and which is at the same temperature as the wedge surface. 

The dimensionless heat-transfer parameter q/q^ (or, since the 
wall temperatures are equal, the equivalent ratio of heat-transfer 
coefficients) is commonly used in shuttle application. Heat-transfer 
distributions were calculated for wall temperatures of 394°K and of 
1640°K. 
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The dimensionless heat-transfer distributions computed using the 
real-gas code are compared with the perfect-gas solutions and the var- 
iable-gamma solutions in Figs. 4a and 4b, respectively. The wall- 
temperature variation had little effect on the heat transfer. The heat- 
transfer distributions for the variable-gamma solution correlate closely 
with the real-gas heat-transf er distributions. The difference in the lo- 
cation of the interaction-perturbed region contributes to the only sign- 
ificant difference between the perfect-gas heat-transfer and the real-gas 
heating (and the difference is magnified by the expanded scale). However, 
for engineering applications, both the perfect-gas distribution and the 
variable gamma distribution are in satisfactory agreement with the real- 
gas distribution. Similar correlations were found for the highest velo- 
city solution. However, near the wing root, the difference between the 
perfect-gas heat-transfer and the real-gas heating was slightly greater 
at the highest velocity (i.e., condition 3 in Fig. 9). 

The local increases in the real-gas heat-transfer distributions which 
are evident at the beginning of each flow region are due to the local ac- 
celeration of the invlscid flow. The waves of the expansion fan produce 
a step-function decrease in pressure and a corresponding step-increase in 
the local velocity at the edge of the boundary layer. Thus, the nonsimilar 
boundary-layer solutions yield local increases in heating due to the local 
velocity gradient. There are no locally severe heating rates, which could 
cause design problems, indicated either in the real-gas solutions or in 
the perfect-gas solutions. However, the flow model for the calculations 
does not include imbedded shock waves or other three-dimensional flow 
phenomena, which might occur near the wing-root fairing. Such flow 
phenomena caused local increases in heating to delta-wing arbiter configura- 
tions^. 



13 


The location and the extent of the interaction-perturbed region of the 
"wing leading-edge" for the real-gas solutions are presented in Fig. 5 as a 
function of leading-edge sweep. The limits of the band represent the inter- 
sections of the limits of the centered expansion fan, i.e., region 4 of Fig. 
1, with the leading edge. The locations are presented as the distance from 
the junction of the two wedges (S) divided by the radius of a reference 
sphere (R, which is equal to 0.0027m). Calculated locations from the real- 
gas solutions are presented for sweep angles from 60° down to the minimum 
sweep angle for which a Type VI pattern exists. The minimum sweep angle 
decreases as the free-stream velocity increases. For the higher velocity 
entry-condition, the Type VI pattern is possible for sweep angles as low 
as 27°. The perturbed region moves inboard toward the wedge junction as the 
velocity increases. 

The effect of the free-stream velocity on the nondimensionalized heat- 
transfer rate is indicated in Fig. 6. Calculations are presented for region 
3 (upstream of the expansion fan) and for region 5E (downstream of the ex- 
pansion fan). Since the interaction-perturbed region moves inboard as the 
velocity increases, the point in region 3 is near the inboard edge of the 
interaction region for the highest velocity. Correspondingly, the point of 
region 5E is near the outboard edge for the lowest velocity. For both the 
perfect-gas solutions and the real-gas solutions and for both wall tempera- 
tures, the nondimensionalized heat-transfer rate increases significantly 
with velocity. The dimensionless heating for the higher entry velocity is 
roughly twice that for the wind-tunnel condition for region 3, somewhat less 
for region 5E. This implies that one should not extrapolate wind-tunnel 
data directly to flight conditions. Instead the wind-tunnel data should 
be used to construct a viable model of the flow field. The flow-field 
model can then be used to generate the required aero thermodynamic environ- 

ment at the conditions of interest. 
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The shock- interact ion geometry, the surface-pressure distributions, 
and the heat-transfer distributions are presented in Figs. 7, 8, and 9, 
respectively, for the highest velocity condition (i.e., condition 3) with 
a sweep angle of 60°. The correlations between the perfect-gas solution, 
the variable-gamma solution, and the real-gas solution are similar to those 
observed in the previous figures for the middle velocity condition. 

As can be seen in Fig. 7, the interaction-perturbed regions are virtu- 
ally the same for the variable-gamma solution and for the real-gas solution. 
The location for the perfect-gas solution is only slightly inboard. The 
shock layer for the real-gas solution is slightly thinner than that for the 
variable-gamma solution and markedly thinner than that for the perfect-gas 
solution. 

The relations between the various solutions for the leading-edge 
pressure-distributions obtained at the highest velocity correspond to those 
noted for the middle velocity solutions. For a given region of the expansion 
fan, the perfect-gas pressure is in reasonable agreement with the real-gas 
value. The difference between the perfect-gas distribution and the real-gas 
distribution is accentuated by the expanded scale of Fig. 8. Even though 
the real-gas solution was used to specify the input values of gamma for the 
variable-gamma solution, the variable-gamma pressures do not match the real- 
gas values as well as the perfect-gas pressures. 

As can be seen in Fig. 9, the wall-temperature variation has no signifi- 
cant effect on the theoretical heat-transfer distributions. The variable- 
gamma distribution closely follows the real-gas distribution. The differ- 
ences between the perfect-gas heat-transfer and the real-gas heating are 
greatest in regions 3 and 4AW, i.e., near the junction of the two wedges. 
Downstream, the correlation between the real-gas solution and the other two 
solutions is similar. 

Theoretical solutions have been obtained for a variety of sweep angles. 
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The solutions of the Type VI shock-interaction for flow condition 2 are pre- 
sented in Figs. 10 through 12 with a sweep angle of 40°. These solutions 
exhibit several interesting characteristics. For this flow condition, the 
shock-interference pattern undergoes critical changes as the sweep angle 
varies. No perfect-gas solution is presented. This is because, for this 
wedge geometry, the perfect-gas relations require a Type V pattern. Even at 
hypersonic speeds, perfect air can not turn through a single, linear shock 
wave parallel to a 50° wedge^. Also note that the curved shock region near 
the reflected waves is not presented in Fig. 10a. This is because the numeri- 
cal scheme used in the present real-gas code does not yield the continuously 
curved shock required to turn the flow from the free-stream direction 
parallel to the flow in region 4E. However, a solution still exists for 
region 6 and, therefore, region 5E. Therefore, the shock-interference pattern 
is still basically a Type VI pattern. The curved shock would, of course, be 
produced in the actual flow. But since the velocity gradients and the pressure 
gradients of curved streamlines are not modeled numerically in the current 
code, the code outputs for this situation; "curved shock not modeled". It 
has been noted that the geometry for the variable-gamma solution is only 
slightly different than the real-gas geometry. However, because of the differ- 
ence, the complete variable-gamma geometry (within the assumptions of the 
numerical code) is generated (Fig. 10b). 

As was noted at the larger sweep angles, the variable-gamma pressure do 
not correlate exceptionally well with the real-gas values (Fig. 11), but the 
heat-transfer rates do correlate well (Fig. 12). 

To simulate the large density changes which exist across a shock wave 
in hypersonic flight, experimental Investigators often make use of wind 
tunnels for whxch the test gas has a relatively low value of gamma, e.g. , 
ref. 13. Using tetrafluoromethane (CF^) as the test gas, the free-stream 
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specific-heat ratio for the Langley facility varies from 1.17 to 1.31 at 
Mach 6 (ref. 14). Therefore, theoretical solutions of the Type VI shock- 
interference pattern for the double-wedge configuration with A = 60° have 
been computed for a perfect gas with y = 1.2 throughout the flow-field. 

These solutions are designated "constant-gamma" solutions. 

The location of the interaction-perturbed region is presented in 
Fig. 13 as a function of the free-stream velocity. Reviewing the legend: 

(a) "perfect-gas" uses the perfect-gas relations with y = 1.400 
throughout the flow field, 

(b) "constant-gamma" uses the perfect-gas relations with y “ 1.200 
throughout the flovj field, 

(c) "variable-gamma" uses the perfect-gas relations with y^, y^> and 

y^ specified from the real-gas solution (y^^ = y^ and Y 3 ^ ^ Y 5 ) » 

and 

(d) "real-gas" uses the equilibrium air properties to describe the 
gas behavior. 

Since y for air is, in fact, essentially 1.4 throughout the flow field for 
the lowest velocity considered, no variable-gamma solution was obtained for 
this condition (i.e. , condition 1). With the exception of the constant- 
gamma solution, the various solutions provide similar locations of the 
interaction-perturbed region over the velocity range considered. At the 
higher velocities, the interaction-perturbed region is relatively inboard 
for the constant-gamma solution. For a given geometry of the double-wedge 
configuration, the location of the interaction-perturbed region is a func- 
tion of the bow shock-wave angle (generated by the initial deflection) , of 
the leading-edge shock-wave angle (generated by the second deflection) , and 
of the Mach numbers in regions 3 and 4 (which determine the expansion waves). 
The shock angle for y = 1.2 is nearer the surface than is the shock angle for 
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Y = 1.4. Thus, since the actual y is essentially unchanged for the first 
deflection angle, the bow shock-wave for the constant-gamma solution inter- 
sects the leading-edge shock nearer the root than is the case for the real- 
gas solution or for the perfect-gas solution (which yield essentially 
identical results). The low-velocity, constant-gamma solution for the in- 
teraction region correlates with the other solutions because both the bow 
shock angle and the leading-edge shock angle are relatively small causing 
the shock: shock intersection to be relatively outboard. 

The pressure distributions along the wing leading-edge are compared in 
Fig. 14. For the wind-tunnel condition, i.e., condition 1, the perfect-gas 
solution (y = 1.400 everywhere) yields flow conditions in regions 3 and 6 
which are identical to the corresponding flow conditions computed using 
the real-gas code. Differences occurred in the perfect-gas solution and 
the real-gas solution for the local pressures and the interaction locations 
of the expansion fan. Because regions 3 and 6 are identical and because 
the differences between the two solutions are attributed to the difficulty 
in using the Mollier charts at these low temperatures, the perfect-gas 
solution with y = 1.400 (Fig. 14a) represents the actual flow. The con- 
stant-gamma solution yields pressures which are significantly lower than 
the other two solutions. As noted previously, the angles both for the 
bow shock and the leading-edge shock are relatively small for y = 1.200 
and, therefore, the interaction region matches that for the other two 
solutions. 

For the middle velocity condition (Fig. 14b), the shock- interaction 
region for the constant-gamma solution is relatively inboard. It has been 
established that the difference is due to the fact that the initial shock 
wave is weak and, therefore, does not significantly alter gamma from its 
actual free-stream value of 1.400. When comparing the pressure from a 
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given region, the constant-gamma value is lower than the real-gas value by 
approximately the same amount that the perfect-gas value is higher than the 
real-gas value. 

For condition 3 (Fig. 14c), the constant-gamma pressure in a given 
region is in very good agreement with the perfect-gas pressure for that region. 
For a given region, the real-gas solution yields pressures somewhat lower than 
the two other solutions. Again, however, because the bow shock-wave generated 
by the initial turning of the flow is much closer to the body when y =* 1.200 
than it is for the real-gas solution or for the perfect-gas solution, the 
constant-gamma interaction region is markedly inboard. 

The heat-transf er distributions for the leading-edge of a wing with 
60^ sweep are presented in Fig. 15. The constant-gamma solution is compared 
with the perfect-gas and the real-gas solutions. Because the free-stream 
values of temperature and of pressure were used to define the flow condition, 
negative heat transfer, or cooling, existed when y = 1.200 for the wind- 
tunnel flow condition. Thus, heat-transfer distributions are not presented for 
the wind-tunnel condition. Heat-transfer distributions are presented for con- 
ditions 2 and 3 in Figs. 15a and 15b, respectively. The comparisons between 
the various theoretical solutions are similar for both flow conditions. The 
differences in the heat-transfer distributions are due principally to the 
differences in the locations of the interaction-perturbed regions. These 
differences are greatest at the higher velocity. The perfect-gas solutions 
(y = 1.400 throughout) compare more favorably with the real-gas heat-trans- 
fer distributions than do the constant-gamma solutions (y = 1.200 throughout). 



CONCLUDING REMARKS 


Using a two-dimensional flow model of the Type VI shock-interaction 
pattern, the aero thermodynamic environment has been calculated for a 
"simulated" wing leading-edge of a delta-wing orbiter. Calculations have 
been made for velocities from 1167 m/sec to 7610 m/sec for perfect-gas pro- 
perties, for constant-gamma gas properties, for variable-gamma gas pro- 
perties, and for real-gas properties. Based on the calculations of the 
present study, the following conclusions are made. 

1. Free-stream flight conditions were found to produce Type VI 
interaction patterns for effective wing leading-edge sweep 
angles as low as 27^, when the real-gas effects were considered. 

2. Perfect-gas solutions for the flow geometry and the pressure 
distribution were in good agreement with the real-gas solutions. 
The use of effective gammas did not adequately represent real- 
gas effects in the surface-pressure distribution. 

3. The correlation between the perfect-gas solution and the real- 
gas solution for the heat-transfer distribution was essentially 
independent of the wall temperature, but depended on the free- 
stream velocity. The heat-transfer distributions for the 
variable-gamma solution correlated closely with the real-gas 
heat- transfer distributions. No locally severe heating rates, 
which would cause design problems, were found. 

4. When the local heat- transfer rates were nondimensionalized 
using a current shuttle design parameter, the dimensionless 
heat transfer increased significantly with velocity. The in- 
crease occurred both for the perfect gas solutions and for the 
real-gas solutions and for both wall temperatures. Thus, one 
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should not extrapolate wind-tunnel data directly to flight 
conditions. Instead the wind-tunnel data should be used to 
construct a realistic model for the flow-field, which can be 
used to generate the required aerothermodynamic environment. 

5. The shock-intersection geometry and its effect on the local 
flow-field was found to be a complex function of gamma. The 
interaction-perturbed region depends on the *‘bow'* shock-wave 
angle , the ”leading-edge'^ shock-wave angle , and the expansion 
wave angles (or, equivalently, the local Mach numbers). Al- 
though gamma often had a significant effect on the locations 
of the interaction-perturbed region, for a given flow condi- 
tion, the location differences were never severe. This con- 
sistent correlation occurred because the shock angle dis- 
crepancies tended to be compensating in many cases and be- 
cause the interaction occurred relatively near the wing root 
so that significant differences in angle resulted in relative- 
ly minor differences in length. 
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Figure 1. - Flow model of the Type VI shock-interference pattern for a double wedge. 
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Figure 4. - The heat-transfer distributions for the wing leading edge; 
= 4330 m/sec, P^ = 0.333 mmHg, T^ = 273°K, A = 60° 
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(b) Real-gas compared with variable-gamma 
Figure 4, - Concluded. 
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Figure 5. - The location and the extent of the interaction-perturbed region 
of the wing leading-edge as a function of leading-edge sweep (real-gas 
solutions) . 
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(a) Real-gas compared with perfect- gas 

Figure 9. - The heat-transfer distributions for the wing leading-edge; U^ = 7610 m/sec, P^ = 0.0268,mmHg, 

T = 195°K, A= 60°. 
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(b) Real-gas compared with variable-gamma 
Figure 9. - Concluded. 


Figure 10. - Calculated geometry of the Type VI shock- 
interference pattern; = 4330 m/sec, = 0.0333mmHg, 
T = 273 K, A = 40 . 
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Figure 12. A comparison of the real— gas and the variable— gamma heat— transfer distributions for the 
wing leading-edge; = 4330 m/sec, - 0.333 mmHg, = 273°K, A = 40°. 
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Figure 13. - The location of the interaction-perturbed region on the 
wing leading-edge as a function of velocity (A = 60°). 
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Figure 14. - Continued. 
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Figure 15. - The heat-transfer distribution along the wing leading-edge for A = 60°. 












APPENDIX A. - GENERAL DESCRIPTION OF PERFECT GAS CODE 


The general solution procedure for the perfect-gas code has been dis- 
cussed previously. This section will provide a more detailed description 
of the governing equations and their place in the numerical routine. Input 
data for the perfect-gas code consists of the free-stream flow-field conditions 
of region 1 and the model geometry. The procedure used to calculate flow 
conditions in regions 2, 3g and 6 is to first call s\±>routine DELTAK, which 
solves for the shock wave angle, and second to call subroutine PTHETA, which 
calculates the flow conditions of regions 2, 3, and 6 using the computed shock 
wave angle. 


Flow-Field Conditions in Regions 2, 3, and 6. The s\±>routine DELTAK uses the 
known turning angle 6 to calculate the shock wave angle 0 in the following 
equation : 


6 = tan 


-1 


tan 0 


(y+i)m^^ 


-1 


( 1 ) 


where the subscript 1 denotes conditions upstream of the shock. A half- 
interval, iteration procedure is used to determine 0. The half-interval 
method is started by assuming 0 to be an average between a lower limit, 
which is equal to 6, and an upper limit, equal to 90^. This average shock- 
wave angle is used to calculate the corresponding 6 from equation (1). A 
comparison is then made between the calculated 6 and the actual 6. If this 
comparison is not within a prespecified tolerance (i.e. , 0.0001 radians) 
either the lower, or upper limit of 0, depending upon the comparison, is set 
equal to the previous iteration’s average value of 0. A new average value 
of 0 is computed using the new limit, and the procedure is repeated until the 
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calculated 6 equals the actual 6 (within the prespecified tolerance). Once the 
shock-‘wave angle has been determined the subroutine PTHETA uses this 0 to cal- 
culate the pressure ratio by the following: 


2YM^^sin^6 - (y-l) 

_ , 5 


( 2 ) 


The subscript 2 denotes the conditions downstream of the shock. Thus, to solve 
for the flow-field conditions in region 3, conditions in region 2 are actually 
the conditions upstream of the shock , and conditions in region 3 are the con- 
ditions downstream of the shock. The density ratio and temperature ratio are 
then calculated as function of 5 as follows : 


^ = (Y+l)g t (y-1) 
(Y“1)? + (y+1) 


(3) 


^ - r (y-1)^ (y+1) 

^ (y+1)5 + (y-1) 


(4) 


The Mach numbers of regions downstream of the shock wave are calculated from 
the equation: 

1/2 


M2 = 


[(y+D? + (y- 1)] - 2(?^-l) 

5[(y-i)5 + (y+i)H 


(5) 


The pressiire coefficient of each region is foimd from the equation; 

P-P, 


'^P = 


C6a) 


0.5y^Pj_M^ 

which, when combined with equation 2, yields the following relation for re- 
gions 2 and 6 : 


4(M ^sin^e-1) 
■' — ^ 

"P 2 

(y+1) 


(6b) 


The procedure to compute the stagnation conditions for regions 1,2, and 



and 3 is to assume the flow decelerates isentropically to zero velocity. The 
equation for calculating the stagnation pressure is: 



where these calculations are carried out for i = 1 (the free-stream) , 2, and 
3 (which serves as the stagnation pressure for the isentropic e^ansion). 

Expansion Fan. The flow in region 3 is assumed to accelerate isentropically 
to subregion 5E. Region 3 is uniquely determined, as described above. Since 
the streamlines in the subregions from subregion 5E through region 6 are 
straight and parallel to the surface, i.e., not curved, 


dn 


0 


Thus , the static pressure in region 5E is equal to the static pressure in 
region 6, which is known. The isentropic assumption requires that the stag- 
nation pressure in subregion 5E is equal to the stagnation pressure in region 
3, which is known. Thus, the flow in subregion 5E is uniquely defined. 

The flow field in the expansion fan is calculated in the subroutine 
EXPAN and the locations of the intersection points of the right running and 
left running waves of the fan are calculated by the subroutine INTRST. Since 
the waves are assumed to be linear, the subroutine INTRST requires the know- 
ledge of two initial points and the angle between these points to the point 
of intersection. The intersection point can then be calculated using linear 
relations. The required angles are the shear layer angles and the expansion 
wave angles (i.e., Mach waves). 
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For the isentropic expansion from region 3 to region 5E, the Prandtl- 
Meyer expansion equations are used. First, the Prandtl-Meyer angle v for 
region 3 and for suhregion 5E is calculated as follows : 



(When using the variable -gamma option of the perfect-gas code, v both for 
region 3 and for subregion 5E is calculated using the gamma for region 3). 
The difference between these two Prandtl-Meyer angles is divided into ten 
equal parts to give the five waves in region 4 and the five waves in region 
5. When crossing the right-running waves of the centered expansion fan, 
the change in the Prandtl-Meyer function is related to the change in the 
flow direction by 

dv = de (9‘-a) 

When crossing the left-running waves reflected from the wall, 

dv = -de (9-b) 

Thus, there is no net change in flow direction, satisfying the condition 
that the flow in subregion 5E is parallel to the surface. The subroutine 
EXPAN calculates the local Mach numbers of the expansion fan from which the 
local pressures can then be calculated. The local Prandtl-Meyer angle is 
calculated by adding one tenth of Av to the v for the previous region. Then 
using equation (8), the local Mach number can be calculated using the half- 
interval interation method. After all the local Mach mmibers have been cal- 
culated, the local pressures are calculated using the following equation: 



( 10 ) 
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The surface pressures have been calculated for the interaction between 
the right-running waves of the expansion fan (region 4) and the reflected 
left-running waves (region 5). See Figure A-1. The flow in subregion 4A 
is directed away from the surface by the angle O.lAv. The expansion of the 
flow from 4A so that it is parallel to the surface in subregion 4AW is ac- 
complished along right-running characteristics. Thus 


^4AW ^4A ^®4AW ®4A^ 


( 11 ) 


A similar procedure is used for subregions 4BW, 4CW, and 4DW. Once the 
Prandtl-Meyer functions for these subregions are known, the remaining pro- 
perties are calculated using EXPAN. 


Wing Leading-Edge Heat-Transfer-Rate Calculations. After the wall pressinres 
in regions 3, 4, and 5 have been found, the final step is to calculate the 
heat-transfer rate on the "wing leading-edge" in these regions. The sub- 
routine ERTQDOT calculates the desired heat-transfer rate. The technique 
used to calculate the heat-transfer rate employs the Eckert's reference tem- 
perature equation for a laminar boundary- layer , as follows: 

q = 0.332(Re (Pr*)°'^^\" (T - T ) (12) 

This heat-transfer rate is divided by the reference stagnation-point heat- 
transfer rate to give a non-dimensionalized heat-transfer for use in corre- 
lation. 

Calculation of the Stagnation Conditions Behind a Normal Shock and of the 
Reference Heating. In addition, the stagnation pressure across a normal 
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shock wave (which of course does not exist in the flow-field) is calculated 
for region 3 and for region 6 using the equation: 



(y+i)Mj_^ 

Y 

Y-1 

r Y+1 1 

2 


12yM^^-(y-1)J 


1 

Y-1 


(13) 


with the proper subscription of pressures and Mach numbers. The stagnation 
temperature is constant throughout the regions and is calculated by the 
adiabatic perfect-gas equation: 


T 


t 



(14) 


These normal shock values of stagnation temperature and of stagnation pressure, 
in region 6 are used to calculate the stagnation point heat-transfer rate to a 
reference sphere. 


The following equation is used: 


,0.4 


, 0.1 


q^ . = 0.64(p u )'^’^ (p y (C T - C T ) 

^t,ref t t w w pt t pw w 


du 

dx 


10.5 


(Pr)'^ (15) 


where 



1 


R T^ 
gas t 


R 


(16) 


and the subscript t refers to stagnation temperature, the subscript w refers 
to wall temperature , and R is the radius of the reference sphere . 


Calculation of the Shock-Wave Angles and the Intersection Points of the 
Shock Wave with the Expansion Fan. Once the flow in subregion 5E is defined, 
the isentropic expansion flows inside of the shear layer in regions 4 and 5 
is defined. It is not possible to match both the pressure and the flow direc- 
tion across the "shear layer" which divides subregion 4ES, and subregion 4E 
(and, sequentially, 5AS, 5A; 5BS , 5B; 5CS, 5C; and 5DS, 5D). The actual flow 
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in subregions 4ES through 5DS is believed to be more complex than the flow- 
field model allows. The required pressure adjustment is assumed to be ac- 
complished by the pressure gradient across curved streamlines: 


dP _ pu^ 
dn ” r 


(17) 


For the flow conditions studied, the radius of curvature is very large, i.e. , 
only slight streamline curvature is needed. The subroutines DELTAK and PTHETA 
are used with the flow angle to generate approximate values of the shock -wave 
angle and the flow conditions in the subregions between the shock wave and the 
shear layer, i.e., subregions 4ES, 5AS , 5BS, 5CS, and 5DS. The intersection 
of the shock wave and the left-running expansion wave is computed using the 
subroutine INTRST. The expansion waves are assumed to be linear from their 
intersection with the shock wave. Because the pressure decreases in the ex- 
pansion fan subregions, the shock -wave inclination decreases for each sub- 
sequent calculation. This results in the "curved" shock wave characteristic 
of the Type VI pattern, when the reflected waves interact with the wing lead- 
ing-edge shock (see Figure 1). 
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NC - number of cases 

FSMACH - free-stream Mach number 

2 

PINF - free-stream pressure (ibf/ft ) 

TINF - free-stream temperature (°R) 

LAMDAS - the deflection angle, i.e. , complement 


of sweep angle ( ) 


EGAS - gas constant (1716 


ft 


2 o„' 

sec R 


GAMMA - free-stream gamma (usually 1.4) 

P2 - pressure in region 2 (Ibf/ft ) 

(optional) 

THETA - first shock wave angle (°) 

(optional) 

DELTA - angle of first wedge (°) 

KTHETA - does not equal zero if DELTA, 6, is known, equals 
zero if P 2 or THETA, 6, are known. 

KN0WN - does not equal zero if DELTA , 6 , is known , equals 
zero if DELTA, 6,. is unknown. 

TW - wall temperature (°R) 

RN0SE - radius of reference sphere (ft.) 

PRN0 - Prandtl number 


CPO, CPI, CP2, CP3, CP4 - specific heat coefficients in the 

equation 

BTU ^ 
slugOR j 


C ^ + C ,T + C ^T^ + C ,T^ + C „T^ 
pO pi p2 p3 p4 




P 

XO - x-coordinate of "nose" (ft.) 

YO - y-coordinate of "nose" (ft.) 

X30 - x-coordinate of wedge intersection point (ft.) 
Y30 - y-coordinate of wedge intersection point (ft.) 



Card # 8 GAMMA3 - gamma in region 3 (can be used to "simulate” real gas 

effects ) 

GAMMA6 “ gamma in region 6 (can be used to "simulate" real gas 
effects ) 

Card # 9 VISSO, VISSl, VISS2 , VISS3, VISS4 - viscosity coefficients in 

the equation 

y = + y^T + y^T^ + y^T^ + y^^T^ 

(Note: Since the Prandtl number for air is approximately 0.7 for most 

temperatiires and pressures, 

K = y C /Pr 
P 

C and Pr are assumed constant for perfect air. Since y is accurately cal- 
P 

culated with the above polynomial, a reliable value for K is obtained). 
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HHOGHAM SH0^lNT ( I NPUT f OUTPUT ) 


PKOGHAM will calculate SH0CK'*'AVE intersections ANO The 
FLOW field conditions HEhImD the SH0C^WAVES GIVEN FrEE- 
stheam Conditions***** 

KFAL MlS» MiJf M3* M6» LAMOAStMgS 
HEAD lilt nC 
111 *'0RM^T(l5) 

UO 99 IC=1*NC 

head 101» FSMaCHi pinf* tinf* lamoas 

101 Format (AE i^.b) 

HEAD 102» HUAS» OAMMA 

102 Format ( 2Eif^»5) 

head 103, pdt THEtA» delta, KTHETA, KNOWN 

103 F0RMAT(3E12.&,2IS) 

HEAD 10A,TWiHN0SE*PRNO 

104 F0RMAT(3E12*b) 

HEAD i05,CP0,CPlfCP2,CP3*CP4 

105 Format (5E12, 5) 

HEAD 106, ''U,YO*A30,Y30 

106 rORMAT (4E12#5) 

HEAD 107, GAMMA3t GAMMA6 

107 FORMAT (2E12.S) 

IHETAsTHETA/57.296 
UELTAsDELTA/57.296 
LaMDAssLAMDAS/57,296 
IF (known, EU,0) GO TO 10 

*****calcdlating FLOW conditions in region ?****« 

Call deltak(fsmach»gamma»oelta,theta) 

10 Call PTHETA(FSMACH,M2,PInF,P2»PHAT21 ,DRAT21,TRAT21,THEtA,OELTA, 
lPCOEFfGAMMA»KTHE| A,KN0WN) 

hhOINFsPINF/ (RGAS*TINF ) 

P2sPRaT21*P1nf 
’>2sTraT2i*T inf 
Hh02sdRAT21*HhOINF 

*****CALCULATING FLOW CONDITIONS IN REGION 3«**** 

UELTA3sLAM0AS-DELTA 
Call DtLTAK(M2»GAMMA3 , deltas, THETA3) 

Call PTHETA(M2tM3,P2,P3,PRAT32,URAT32,TRAT3?,THETA3,l)ELTA3,PC0EF3, 
1GAMMA3,KTHEIA,KN0WN) 

C 

PRAT3lsPRAT32«PRAT21 

rRAT31sTRAT32«THAT21 

URAT3lsDRAT32*DHAT2l 

P3apRAT3l*PlNF 

r3sTRAT3l*TlNF 

HH03*DRAT3I*RH0INF 

PCOEF3s2,0*(PRAT31-U0)/(GAMMA3*FSMACH*FSMACH) 

C 

G3P1sGAMMA3^1,0 

G3M1sGAMmA3-1,0 

PMF3aSQRT{63Pl/G3Ml)*ATAN(SQRT (( GSM 1/63^1 )*( M3*M3-l # 0 ))) -AT AN ( SORT 
1 (m3*M ^-1,0) ) 

PTEaP3*( (1.0*(G3M1/2,0)*M3*M3 )**(GAMMa3/G3m1) ) 
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*****CALcULATlN(j flow conditions IN REGION 
OFLTAftsLAMUAS 

CALL 0ELTAK(FSMACH»GAMMA6,0ELTAb,THFTAb) 

Call PTHETA(FSMACh,M6*PINF,P6»PRAt 61 ♦nRAT61 *TRAT6l*THETA6»0tLTA6, 
1PC0EF6*GAMMA6 *KTHETA#KNOwN) 

P6sPHAT61*P1NF 

FbsTKAT61*TiNF 

HH06=DRAT6i<»HH0INF 


***«*CaLCOLaT ING FLOW CONDITIONS IN REGION 5E****« 


PSEsPb 

PRATSEsPTE/PbF 

tXM5E*SQRT (^.o* ( (PRAT5E) ** (G3M1/GAPMA3)*! .0) /G3M1 ) 

PmFSEsSORT (G3PI/g 3M1)*ATAN(SQHT( (G3Ml/G3Pl)*(£XMSE*EXM5E-i.U) ))- 
1A7 AN(SQRT(EAm5E*EXM5E-1*0) ) 

UEXPsl.,1*(PMF5E-PMF3) 

***»«C0MP0TiNG stagnation CONUIT IoNS*«*** 

MlSsFSMACH*FSMACH 

M?SsM- 

GmIsGAmMA-1 «0 
VjPIsGAMMA*1*0 
EXPONsGAMMA/GMI 

PT1sPinFa( (1.0^GM1*M1S/2#0)**EXPOn) 

PT2=P?*( ( 1.04gM1*M2S/2*0)**EXPON) 

PT3sP2* ( (6P1*m 2S/2,0) **EXP0N) * ( (GPl / (2. 0*GAMMA*M2S-6M1 ) )*♦(!, 0/GMl 
IM 

PTbsPINF* ( (GPl«MiS/2.0)**EXPON)<» ( ( 6P1/ (2* 0«GAMMA*M1 S-Gm1 ) >**(1,0/ 
IGMl) ) 

TTsTlNF*( i.U^GMl *Ml S/2.0) 
fTF-sTT 

Call QOOT (Plo,HGAS*TW*TTfRNOSE*CPo*CPl *CP2,CP3#CP4fOrSPR,PRNO) 


tiOW ANU WING SHOCKWAVE INTERSECTION***** 
aGL3*THETA3+DELTA 

CALL lNTRST(X0»Y0,THETA,AGL3tA30*Y30*Xl»YI) 

*****FIND intersection FOR THE FIRST EXPANSION WAVE*«*** 

*****EQUaTION for the angle of line 4A***** 

WAVE3sASlN(i.n/M3) 

A6L4AslAmDAS-wAV£3 

Call INTRST(XI*Yi,AGL4A»LAMOAS,X30tY30*)^AA,Y4A) 
rMTA4AsLAM0AS*0EXP 

«««««LINE 4b«**** 

Call EXPAN(PMF3fDEXP,GAMMA3#PHF4A,EXM4A#M3) 
waVE4AsASIN(1.0/EXM4A) 

AGL4aaLAM0AS+DEXP«WAVE4A 

Call INTRST(XI.YI*A6L4d,LAMOAStX30tY30*X4B*Y4B) 
fRTA48sTHTA4A4DEXp 



ono ono ooo ono onn non nnr> no 


*«#**LINE 4C*«*** 

Call expan (PMF4A»0EXPt(iAMMA3fPMF4R*EXM4btEXMAA) 

WAVE4rtaASlN(l,O/tXM40) 

AGL4CsLAMUAb*2.0*DEXp-WAVE4B 

CALL INTrST ( AI»Y l,AbL4C*LAMDASfX30»Y30»X4CtY4C) 
IHTA4CsTHT<‘‘»B*UfcxP 

Call EXPAN(PMF4t»»DEXP»(3AMMA3»PMF4C»FXM4CtEXM4b) 

wAVE4CsAsIN(1,0/EXM4C) 

AGL40*LAMDAb*3.0^DEXp-WAVE4C 

CALL INTrST (XI»Y i»AGL4U*LAMDAStX30»Y30*X4D,Y4U) 
•hTaaOs^H^A^C+DEXr 

**«*«(_ 4t^**** 


CALL EXPaN(PMF4C*DEXP»WAMMA3iPMF 4D»FXM4U#EXM4C) 
wAVE4DsASlN ( 1.0/EXM4D) 

AGL4ExLAMt)Ai>*4 .0*DEXP-WAVE4D 

Call INTrST (XI.Y l*A6L4tfLAMDAS,X30»Y30«X4EfY4E) 
ImTAaEsTh rA4U*fJLAP 




Call EXPaN (PmF 40»DEXP »bAMMA3»f^MF4E»EXM4E*EXM4U) 
i«A\/E4EsAsIN( 1.0/FXM4E) 

AGL5AsTHTA4t*WAvt4E 

Call lNTRST(XI*Yi,THTA4Ef AGL5AfX4A*Y4A»X5A,Y5A) 
ThTA5AsThTA4E-UEar 

*«<hh»LINE Sb<HH>** 


CALL EXPAN{PMF4t»DEXPtGAMMA3*PMF5A»FXM5A,EXM4E) 
wAVE5AsAslN(l,0/EXM5A) 

AGL5bsTHTA6A>WAvE5A 

CALL INTRST(X5A»Y5A*THTA5A»AGL5b»X4BtY4b*X5a»YbB) 
rHTA58sTHTAbA-UEXP 

*«««« likjE ^q*«*o* 

Call EXPaN (PMF5A*DEXP »GAMMA3*FMFbRiEXM5b»EXM5A) 
WAVE5MsASlN(l,0/tXM5B) 

AGL5CaTHTA5b*wAvE5B 

Call INTrST (X 5B»Y5B*THTA5B»AGL5C»X4C»Y4C*X5C*Y5C) 
iHTASCsTHTAbb-DExp 




Call EXPAN(PMF5b*DEXP*GAMMA3*PMF5C»EXM5CfEXM5B) 
wAVE5CaAslN(l,0/EXM5C) 

A6L5DaTHTAbC4WAv£5C 

CALL INTRST (X5C*T5C*TMTA5C»AGL5L>.X4D*Y4U*X50»YbU) 
THTAbOaTHTADC-DEXP 

*tt***LiNE 5E**<H»* 

CALL EXPAN(PMF5C*DEXP,GAMMA3»FMF5n,FXM5t)fEXM5C) 
WAYEbUaASiNd.O/tXMbD) 

A6L5EaTHTA5U*WAVE5D 

Call INTRST(X5D»V50f TMTA50»AGL5E»X4E*Y4E»X5E#YbE) 



MAVe5EsASlN(l,0/£XMSe) 

fHTABEsTHTAbO-DExP 

C 

C *****CALC0LATING PRESSURES IN THE EXPANSION FAN***** 

C 

P4AsPTE/( (1.0*(G3m1/2.0)*EXM4A*EXM4A)**(GAMMA3/GJM1) ) 

FaBsHTE/ ( U.O* (G3mI/2*0)*ExM4B*Exm4B)**(6AmMA3/63M1) ) 

FaCsPTE/ ( ( I *0* (G 3 MI/ 2 .0) *EXM4C*EXm4C) ♦* (6AMMA3/63M1 ) ) 

FaDsPTE/ ( < 1 <G3m 1/2.0) *EXM4l)*EXM40) ** <6AMMA3/63MI ) ) 

FaEsPTE/ ( ( I .O* (G3MI/2. 0) *EXM4£*EXM4E) ** (GAMMA3/G3MI ) ) 

F5AsPTE/( U.O* (G3M1/2.0)*EXM5A*EXM5 A)**(GAmMA3/G3M1) ) 

F5BsPTE/( ( l.O* (G3 m1/2.0)*EXM58*EXM5B)**(GAMMA3/G3M1) ) 

FSCsPTE/ { ( i .O* (G3M1/2 . 0) *EXM5C*EXM5 C) ** (GAmmA3/G3M1 ) ) 

FSDaPTE/ ( ( i *0^ (G3M1/2.0) *EXM5U*EXMSn) ** (GAMMA3/G3HI ) ) 

C 

C 

C***«*CaLCULATING flow along the wing In the expansion region***** 

c 

ExP4AWsABS(LAMUAS)-THTA4A) 

CALL EXPAN(FMF 4 A*EXP 4 Aw*GAMMA 3 tPMF 4 AW,EXM 4 AW«EXM 4 A) 
tXp4BWsABS (LAM0 AS-THtA4B) 

Call EXPAN(PMF4B»EXP4bW*GAMMA3,PMF4BW,EXM4BW*£XM4B) 
ExPaCWsABSILAmDAS-THTAaC) 

Call EXPaN(FMF4C,EXP4CW,GAMMA3*PMF4CW,EXM4CW,EXM40 
ExP40WsA6S(LAM0AS-THTA4D) 

Call EXPAN(FMF40tEXP4UW»GAMMA3fPMF4nW,EXM4nW»EXM4D) 

C 

F4AWSPTE/ ( ( 1 .O^ (G3MI/2.0) *EXM4AW*EXM4AW) ** (GAMMA3/63M1 ) ) 
F4BWsPTE/( U.0+ (G3M1/2,0)*EXM4BW*EXM4BW)**(GAMNA3/G3mi) ) 
F4CWSPTE/ ( (X.O* (G3M1/2.0)*EXM4CW*EXM4CW)**(GAMMA3/G3mi) ) 
PAOWsPTE/ ( ( I .O* (G 3 MI/ 2 . 0 ) *EXM40w*EXM4DW ) ** (GAMMA3/G3mi ) ) 

C 

C*****C0MPUTING shock wave ANGLES ANQ INTERSECTING POINTS OF REGIONS 
C 4E AmU b***«* 

c 

Call 0ELTAK(FSMACHtGAMMA6fTHTA4E*THT4ES) 

Call PTHETA(FSMACH*EXM4ES.PlNF,P4EfPR4EI» 0R4E1 t TR4E I , THT4ES, 
1UELTAX,PC0FF,GAMMA6»K THETA, KNOWN) 

Call lNTRST(XI,Yi,THT4ES*AGL5A,X4A#Y4AfX5AS*V5AS) 

CALL 0ELTAK(FSHACH,GAMMA6,THTA5A,tHT5aS) 

Call PTHETA (FSMACh,EXM5AS,PINF,P5A,PR5A1, DR5AI,TR5A1,THt&AS, 
10ELTAX,PC0FF,GAMMA6»KTHETA, KNOWN) 

CALL lNTRST(X5ASfY5AS*THT5AS#AGL5H,X4B*Y4B,X5BS,Y5BS) 

Call 0ELTAK(FSMACH,GAMMA6*THTA5B,THTBBS) 

Call PTHETA(FSMACH,EXM5BS,PINF,P5R,PR5BI, OP5Bl,TR501,THTSaS, 
lOELTAX ,PC0FF,GAMMA6,KTHETA, KNOWN) 

Call INTRST (X5BS*Y5BS*THT58Sf AGL5C,X4CtY4C,X5CS,Y5CS) 

Call 0ELTAK(FSMACHtGAMMA6,THTA5C,THT5CS) 

Call PTHETA(FSMACHtEXM5CS,PINF,P5CtPR5Cl# 0R5CI , TRSCl ,THT&CS, 
1UELTAX,PC0FF,GAMMA6 ,kTHETA, KNOWN) 
call INTRST(X5CStY5CStTHT5CS*AGL5D*X40*Y4D*X5DS*Y50S) 
call DELTaK(FSMACH, gammas* THTA50*THT5dS) 

CALL PTHETA(FSMACH*ExM50S,PINF,P5o,PR5D1, DRSOI * TR 5 DI , thtsos, 
1UELTaX*Pc0FF*GAMMA6*KTHETa*KN0WN) 

CALL INTRST(X5DS*Y5DS*THT50S*AGL5E*X4E*Y4E,X5ES*Y5ES) 

C 

ThT4ES»ThT4ES*57.296 

IHT5As»ThT5As*57,296 

1HT5Bs«THT5»S*57,296 

ThT5CS»ThT5CS*57,296 

ThT51)S»ThT5US*57,296 

C 



non 


c 

UF.LTAsDElTA<> 57.296 
rHETAsTHETA«57.296 
UELTA3 sDELTm3«57.296 
1HETA3sTHETA3»57,?96 
UELTA6sI^ELT''6«57,296 
1 hETA6sThETA 6*57#296 
lHTASEsTHTAbE»57.296 

*«***printing the solutions***** 

HRINT 200 

200 i-oRMt4T(iHl) 

HHINT 299* iC 

299 format < 13^»*CASE a**I 2 f///) 

PRINT 201 

201 EORMaT a 6 X»*FREE-STREAM FLOW CONDITIONS**///) 

PRINT 202 

202 F0RMAT(13X**FSMACh*, 11A,*P1*, 14X,«T 1*, 12X,*RH01*, 12X,*PTi** 12X,*TT 
1 *) 

PRINT 217 

217 FORMAT (26X** (LOF/SQFT) **6X»* (RANKINF)**«X** (SLUGS/CUF'T) *»/) 

PRINT 203» CSMACh* pink* TINL* RHoINF, PTl, TT 

203 Format ll0X*tl2.S»3X,E12.b,AX*£l2.5*3X,E12,5,AX,El2.5,4X*E12.5*//) 
PRINT 204 

204 FORMAT (l 6 X**GAS CONSTANTS**/) 

PRINT 205* KGAS* GAMMA* OTSPR 

205 Format (11X**RGAS s**E 12.S*5X**GAMMA s«*tl 2 ,S*bX**QTSPR S**E12*5*// 
I/) 

Print 206 

206 FORMAT (16X**FL0W CONDITIONS IN REGION 2**///) 

PRINT 207 

207 FORMAT tl6X,«M?** 14X**P2/P1*, lbX*«T2/Ti**14X**Rh02/RH0l** lbX*«p(;OEF 
12 «*/) 

PRINT 208* M2, PNAT21* TRAT21» URAT21, PCOEF 

208 FORMATaoX*tl2,5*6X*El2.b*8X*El2.5*9X,El2,5*9X*£l2.5»//) 

PRINT 212»P2*T2»HH02*PT2 

212 Format ( 11 X*»P 2 =«,E12.5*5X**T2 »*,El2,b*5X,*RM02 =«*Ei2.5*5A,#pT? 
1=**E12,5*/) 

PRINT 209* delta* THETA 

209 FoHMmT aiX**DEUTAl2 s**El2,5,5X**THETAl2 a#,El2.5) 

Print 210 

210 format (///* 16 X**FL 0 W CONDITIONS IN REGION 3 **///) 

PRINT 213 

213 format (15X**M3** 14X**P3/P1** 15X**T3/T1** 14 X**RH03/RH0i«*ibX**PC0FF 
lb**/) 

PRINT 208* M3* PRAT31* TRAT31* URaT 31* PC0EF3 
PRINT 2l5*Pb*T3*HH03*PT3 

215 FoRMAT( 11X**P3 a**El2.5*bX**Tb »* ,E 1 2 ,b *5X * *RHo3 a**E 12.5»bX**PT3 
1=**E12.5,/) 

PRINT 218* OELTA3, THETA3, GAMMA3 

218 format (llX,«OELTA23 s*»*E12.5*5X,0THFTA23 **,E12.5*5X ,*GAmMA 3 a*, 
1E12.5) 

PRINT 211 

211 format (///* 16X**FL0W CONDITIONS IN REGION 6 ***///) 

PRINT 214 

214 Format ( 15X**M6**l4X*«P6/Pl«*l5X*«T6/Tl**l4X**HH0O/RHOltt*lbX**PC0EF 
16 **/) 

PRINT 208* M6* PHAT61* THAT61* URAT61* PC 0 EF 6 
PRINT 2I6*P6,T6*HH06,PT6 

216 format (UX**P 6 »**El2.b*bX**T6 =*,E12.5*5X**RM06 s**E12,5*SX**pT6 
ls**El2,5*/) 



HHInT 219* UELTAb, THETAb, 6AMMA6 

219 Format (UX»*UFLTa 16 s*,E12*5*SXf »THFTAlt> s**El2.6*5X,*QAMMAb a«, 

1 tl2.5) 

PRINT 22p 

220 KORMATI/Z/tAoXf^FOlNTS OF INTEREST IN THE FLOW FIELO AND ON THF BO 
lUY (IN INCHES) *♦///) 

FMINT 2l>n»PrEtTTt 

2b0 Format ( l6X**PTt s*,El2,5t5x,»TTe s**Ei2.5f/) 

PRINT 221 

221 FORM-T (l3Xt«lNlT lAL POINTS ON BOW SHOCKWAVE*) 

PRINT 222*^0»YO 

222 FoRMAT(llX»*X s*»El2.b*5X»*Y =**El2.5) 

PRINT 223 

223 FORMAT</tlJX,»wING INTERSECTION POINT*) 

PRINT 222*XJ0*Y30 

PRINT 224 

224 FORMAT (/» 13X,*dO«» SHOCK-WING SHOCK INTERSECTION POINT*) 

P-InT 222*Xi,Yl 

PRINT 225 

225 FoRMAT</*13A,*FIRST EXPANSION WAVE WInG INTERSECTION POINT*) 

PRINT 222fX4A»Y4A 

PRINT 22fc 

226 format i/,13^»*LINE 4H WING INTERSECTION POINT*) 

PRINT 222*X4h,Y48 

PRINT 227 

227 format U* 13X»*LINE 4C wING INTER9FCTION POINT*) 

PRINT 222*X*C»Y4c 

PRINT 240 

240 FORMAT (/* liX,*LlNE 40 wiNQ INTERSECTION POINT*) 

PRINT 222*X*U,Y40 

print 241 

241 format (/♦ lJAf*LIiNE 4E WING INTERSECTION POINT*) 

PRINT 222»X4E,Y4t 

Print 22a 

228 format (/» 1JA»*LINE SA SHEAR INTERSECTION POINT*) 

PRINT 222»XbA,YbA 

PRINT 229 

229 F0RM«T(/»1JA,*HNE 5B SHEAR INTERSECTION POINT*) 

Print 222»xbH,Ybd 

PRINT 23n 

230 Format {/♦Ua,*line sc shear intersection potnT*) 

Print 222»AbC,YbC 

PRINT 242 

242 F0RMAT(/*13A,*LINE 50 SHEAR INTERSECTION POINT*) 

PRINT 222»XbO,YbU 

PRINT 243 

243 F0RMAT(/*13A,*LINE 5E SHEAR INTERSECTION POINT*) 

PRINT 222»XbE,Y5t 

PRINT 23l 

231 F0RMaT(/»13A»*FL0W IN REGION 5E FROM E)^F»ANSI0N*) 

PRINT 232»THTA5t»EXM5E»PbE 

232 FORMAT(llX**THETA s*»E12.5*5Xi*M5E s*,E12.5»5X**P5E s*»El2.5) 

PRINT 233 

233 F0RMATI/*13X,*MACH NUMdERS IN EXPANSION FAN*) 

PRINT 234* tXM4A*£XM4d»EXM4C 

234 Format (11X**M4A s*»E 12.5»5X,*M4B a«*El2.5*5Xt*M4C a*,El2,b) 

PRINT 244*EXM40»EXM4E 

244 F0HMaT(Ux**M4L) «*♦ E12 ,5 ,5x ♦ *M4E a«*El2.5) 

PRINT 235»EaM5A,ExMSB»ExM5C 

235 Format ( llXt*M5A s*»E12.5»5X,*M5B s«,Ei 2*5,SX»*M5C =*,El2.b) 

PRINT 245*EXM5D 

245 FORMAT (UX»*M50 s**El2.b) 



non 


H«INT 236 

ii36 ^OHMMT(/^lJX,*PHtSSURES IN THE EXPANSION FAN*) 

HPINT 237»P4A,P4b,PAC 

237 F0RMAT<UX**P4A ***El2.5*5X*»P4b =«*Fi2,5.5X»*P4C s**El2.5) 

PRINT 2A6*P*IJ,F*E 

246 I-ORM^T ( 11 X»*P 4 U s*,El2,5t5x»*P4E s*,Fi2.5) 

PRINT 230f PbA*P5b,P5C 

238 F0RMAT(11X,*P5A s*,El2.5f 5X.*P5H a*,n2.5,5X,*PbC a**El2.i>) 

PRINT 247*PbU 

247 I-0RMAT<UX**P5U s*«El2.5) 

Print 4oo 

400 h0RMrtT(lHl*l6X»*I^L0l« ALONG THE WING IN THE FXPttNSlON RFGION*.///) 
PRINT 401 - - * 

4ul FORMAT* 13 A**MACH NUMeEHS*t/) 

Print 402tEXM4AwtEXM4BW«EXM4CW«EXM4DW 

402 F0RMAT(llX*»M4Ai(» s*,E12*5f3X**M40w s**E12*5f 3 X**m 4C* a**El2«5*3X* 
1*M40W s**E12,5»//) 

Print 403 

403 FORMAT *l3X**RPESi>URES«»/) 

PRINT 404*P4AW*P48W*p4CW»P40W 

404 Format ( 1 IX»*P 4 Aw s*»El2.5t3Xf*P48w a*,El2*5t3Xt*P4CW s*tE 12 * 5 * 3 X. 
1*P4DW a*, r 12 , 5 ) 

PRINT 300 

300 FORMnT (///» i3X»*F'LOW IN THE SHOCKED REGION OF EXPANSION FAN*i///) 
PRINT 301 

301 F0RMaT(26X**RE6I0n 4ES**//) 

PRINT 302* THT4ES* X5AS* YSAS 

302 format ( l 6 X**THETA a*,El2.5*5X»*X a* ,El2.5*5X **Y s**E12.5) 

PRINT 303* PfUEI* 0R4E1, TR4E1. EXP4ES 

303 F0RMAT(/*16X**Pr a**E 12.b»5X**0R a«,El2*5.5X**TR s**El2.5*bX, 

1*M 3*,El2.b*//) 

PRINT 304 

304 FORMmT(26X**REGIUN 5aS*»//) 

PRINT 302* THT5AS, X5bS* Y5BS 

PRINT 303* PR5A1* ORSAl, TR5A1* EXM5AS 
PRINT 305 

305 FoRM-aT (26Xt*REG10N BBS**//) 

PRINT 302* IHT5BS, X5Cs, YSCs 
PRINT 303* PR581* DR 5 B 1 * TR581* EXM58S 
PRINT 306 

306 F0RMAT(26X**HEGI0n 5CS**//) 

PRINT 302* THTSCSf X50S* Y5DS 
PRINT 303* PR5C1, DR5C1 , TR5C1, EXM5CS 
PRINT 307 

307 F0RMAT(26X**REGI0N 50S**//) 

PRINT 302* THT5US. X5ES* Y5ES 
PRINT 303* PR5D1, DRSDlt TR5DI* EXM50S 

*«h»*h*calculating heat-transfe« along WING-LEaDING eoge***** 

Call ERTODOI (PTE*TTE*P3,GAMMA3*RGaS *X30*X4A*X4b*X4C*X40*X4E* Y3 o* 
lY4AtY4ti*y4C*Y40* Y4EtP4AWfP4BW*P4Cw*P4DW*P5E*QTSPR*PlNF) 

99 CONTINUE 
Enu 

— subroutine U£LTaK(FSMACH*GAMMA*UELTA*THETA) ■! 

c *****SOLVlNU for theta given delta***** 

c 

REAL MIS 
CONVso.OOOl 
MlSaFSMACH*FSMACH 
GPla6AMMA*1.0 



ooool oon I n o DO 


(hETaLsOELTm 

fHETARsl.S^l 

5 (HETA» (ThETal^THETAR) / i.O 

UELT«C*ATAN( 1 ,0/ (TAN (THETA) * ( ( (GP T «MlS > / <2 . 0* (Ml S* ( (SIN(THET a) )**? 

D-1,0) ) )-l.U) ) ) 

UIFF=UELTAC-UEUTm 

iF (ABS (UIFF) .Lt.CONV) GO TO 2 
iF(DlFf) 3f2,4 

3 (HETALsThETA 

GO TO 5 

4 rHETARsTHETA 

Go TO ^ 

2 Continue 
R t turn 

EnI) 

bUHROUTiNE PTHETA (FSRACH*M2*PINF tPafPRATZl *ORAT2l »TRAT2l « IHETA,"" 
lUFLTA»PC0EF«GAHMA,KTHETAtKN0WN) 

*»<hh»SOLVIng FOR FLOW CONDITIONS BEHIND A wEOGt SHOCK GIvEN EITHER 
P2 OR THETA<hhhh» 

HEAL M1S» M2 
MlS=FSMACH*FSMACH 

oplsGAMMA^ 1 ,0 
Gm1sGAMMA-1.0 
IF (KTHETA.EU.O) 00 TO 2 

HPAT21= (2»U«GAMMA*M1S*( (SIN(THETA) ) ** 2 ) -GMi ) /GPl 
IF (KNOWN, NE.O) GO TO 4 

GO TO 3 

2 PRAT2lsP2/PlNF 

IhETa=aSIN (Sort ( (GPX*PRAT21*GMD / (2.0*GaMMA*M1S) ) ) 

3 UELT«sATaN(1,0/(TAN(THETA)«( ( (GP1*H1S)/(2,0*(MIS»{ (SIN(ThETA) )** 2) 
l-l.O) ) )-1.0) ) ) 

4 M2sSORT( ( (MIS* (GP1*PRAT21*GMD ) -2 , 0* ( rRAT2 1 **2- 1 , 0 ) ) / { PRA 121* ( GMI * 
1PRAT21*GPI) ) ) 

URAT21s(GP 1<»PRAT21*GM1) /(GMl*PRAT2UGPl) 

(RAT2lsPRAT2l*( (GM1*PRAT21+GP1 ) / (GPl*PRAT2l ♦GMI ) ) 

PCOEF=(A,0*(M1S*( (SIN (THETA) ) **2) -1 , 0) ) / (GPl*MlS) 

RETURN 

End 

— SUBROUTINE UUOT(RTREF»RGAS,TW»TT,RNOSE»CPO»CPlfCP2.CP3tCP4*OTSPR» . 
IPRNO) 

*«**«SUBROUfINE TO CALCULATE THE REFERENCE HEAT TRANSFER***** 

HHOWspTREF/ (RGAS*TW) 

HHOTsPTREF/ (RGAS*TT) 

VIST=2,27E-08*(TT**1,5) / (TT*198.6) 

''lSWs2t27E-U8*(TW**1.5)/ (TW^l^e.G) 

UUE0X*( (2#0*RGAS*TT)**0,5) /RNOSE 

Cpw=CpO*CPl*TW*Cp2*Tw*TWtCP3*(TW**3 ) ♦CPA*(tw**4) 
CPTsCP0*CP1*TT^CP2*TT*TT+CP3* (TT*3) ♦CP4*(TT**4) 

UELHscpT*TT-CPW*Tw 

UTSPRs0.fe/*( (RHOT*VIST)**0.4)*( (RH0W*VlSW)**0.i)*OELH*(DutUX**0,5) 
1/PRNO 

return 

End 

I., SUBROUTINE INTRSl (XltYl,Al»A2»X2»Y2»Xl*YI) — 

*****SUBROUTInE To find The intersection point OE the shockwaves 
and the intersection points in the expansion FAn***** 


XTs(Y2-r’l^XiwTAw(Al)-A^i«TAw(A;^) ) / ( T AN. ( A i ) -T A-N ( A2 ) ) 
y laYl. ♦ (XT-Xl ) »TAN ( A] ) 
hK TUkN 
End 

bUHROUT ImE E XPAN (PmF 1 ♦ dFXP » GA 'iP A ♦ py FO , E AMACH . EXM i ) — 
«*<MH»SUbPOU ilNit in CALCULATE *iaCh MIWrEhS tm Tut EXHMjSlO''! T 


Um 1=(,5AMMA-1 . i: 

UP j so A MiA A + • () 

HmFO=PMF I *ULXP 
AsSQHT (U pl/bMl ) 

PCONv=0.0U01 
^fM=£NMl 
Hy=l . 1<*STM 

5 ExMACHs (sTM + RM) /<i,0 
£XMs = EXMaCH<»EXi^ACH 

FMFOf7=A<*A Tain (Su'H I ( {£ X ^iS-l . ) ) * jmI /r,P] ) ) - ATam f SUk 1 (EX^-iS-l .0 H 
Pi'IFF=PMFOG-PmHU 

iF (AHS (Pl)IFf ) ,LE.PC0N\/) <30 TO 2 
iF(PiilEF) 3»2»<* 

3 STMsEXMACH 

Go TO 

4 t^MsExMACH 
LjO to S 

2 Continue 

HfTUkm 

Emu 

— bUBRUUTiNE EHTOUOT ( PTE » T I E » P3 ♦ G » « , X T t X4 A ♦ X4B ♦ X4C » X40 , X4E , YI»Y4A»- 
1 Y4P* Y4C» y 40» Y4E»P4A,pifO,P4C,P^f0»pc;F»QTSPR»PTNF) 

*<hhm»subpoui IMF 10 Calculate elow papaE’ETfp'=: aluing vng-l£ ading 
tUoE 

UT MENS I On S ( 15 ': ) »Oi)OT ( IBn ) « PE U 5.j ) , TW ( 1 bO ) . TF ( i bO ) » £'■« ( 1 5o ) ♦ f ST AR ( 
U50) ,VIS2(lbO) , rcownbO) fRHOSdbO) »AF (IbO) ,mfU5u) ,RFunS(l3U) ,TR( 
2150) »PRG ( ISO) ♦ VISE (ISO) t wHUF ( 1 50 ) . RFS ( 1 bO ) ♦onOTR ( 15fM , pRa T ( I 50 ) 
HEAD *PPtTNALL »UFLS, J»LL 
1 Format (3E12.S, 2110) 

READ 12»VIS50,VISS1 »VISS2»VISS3»VtSS4 
12 Format (^El2,b) 

Xi SX1-X4A 

Xpsx I-XAR 
XAsX 1-X4C 
X4SXI-X40 
X5=X1-X^E 
Y jsYI-YAA 
V?sYI-Y4R 

Yt=YI-Y4C 

Y4=YI-Y40 

Y5=YI-Y4F 

S ( 1 ) *n • 0 

SnELTsDELS 

S4AsS0RT (X1*X1 ♦ YHjYI ) /1?.0 
SAdsSQpT (X2*X2+Y2*Y2)/12.0 
SaCsSORT (X3*X3+Y3«Y3) /12.0 
S4U=SQRT ( X4«X4*Y4»Y4) /12.0 
S4E=SoRT (X5*X5*Y5»Y5) /12.0 
UELTS=0.o2*(S4E-S4A) 

RECsgORT (PFn 
UO b ls2t J 
S(I)sS(I-l)*SnELi 



Tw ( I ) sTw all 

iF(S(T) .LE.i>4A) 1^0 TC *»! 
iF(LL#^U,2) (30 ru 47 
IFCLL.EU.I) (jO tu ^7 
briEL l sDEl rs 
a(I)=54A-n.'34<;i)ELT 
( I) =P3 
LL = 1 

on TO 44 
47 Conti fMiJE 

iF (S (1 ) .LE.b4fl) ^0 TO 
IF (S ( I) .LE.04C) 60 TC ^3 
IF (S(I) .uE.b4n) Of) TC 4b 
IF (S ( i ) .1 E.bf^F) 00 TC 46 
HF(I)=PoF 

IF (LI. .EC1.<^) 00 TO 44 
?>nELT=OEl.. 

LL=2 
on TO 

41 PF(1)=P3 

on To 4A^ 

PF(I)=P<»A 
on T'j 44 
^3 ( I ) =P‘»-H 

on TO 4^ 

45 PF(I)=P4C 

on T 44 

46 PF(l)=P4n 
44 C^NTI^'UE 

IF(I)=TTf<»( (PF ( U /PTF ) <»« ( (0-1 • ) /O) ) 

PpG(l)s(pTE/PFa) )44( (<7-1,)/G) 

Em ( I ) =SOpT ( (PPG ( I ) -1 . ) <»2. / ( G-1 , ) ) 

1 STAR ( I ) = .5* (TE ( I ) ♦Tw ( i ) ) ♦.?2 *kEC<MF 1 ( 1 ) ) * I E ( I ) * ( (0-1*)^?.) 

VISE ( j ) = (VIbSO + ViSSl«Tt (I) + VlSSii<>TE (1)*!^ (I ) ♦ V 1 SS3* I TF ( I ) ♦ 

1 VISS4* (TF ( I) <M»4) ) 

VtS2(i) = (Visso 4 V iSSI^^TSTaH ( 1) 4 VxSS2*TSTAH(T)*T«.TaR(1) 4 
1VtSS3<»T5TAp (l)iHbTAR(I)*T5TAH(I) 4 V T 5i’A*TST H ( I ) <»TST A p ( I ) * T ST A R n 
2)»TSTAR ( T ) ) 

TCON ( I ) =f'.2<*<»32 4l76*Vlb2 ( I) /Pn 
HhOS ( I ) =p ( 1 ) / (H^TSTak (I ) ) 

PHOE(I)sPt (1)/(R*TE(I) ) 

AF(i)=S0pr(o<>F?*Tt (I) ) 

UF ( I ) sEM ( I ) »4F ( 1 ) 

KFNOS (I ) =RHUS (1 ) ( I > / (VIS? ( I ) > 

rtFS ( I ) sWHOE ( I ) «Ut ( I ) ( I ) /V I SE ( 1 ) 

1p(I)=TE(1)«(1.U4h£C«{ (6-l.0)/2.0)« (Fm(1)«»p.O) ) 

gnoT { I ) =.332*S0R • (RENOS ( 1 ) ) « (RR<»'»,3.13) *TCOm ( T ) * (TR ( H -TW ( 1 ) ) /SORT ( 
lS( I ) ) 

UDOTh ( I ) sOOoT ( 1 ) /QTSPR 

PRAT ( I ) =PE ( i ) /PlOF 

6 Continue 

PRINT 60 
60 format (I HI) 

PRINT 77 

77 format (26^»*EL0R PARAMETERS ALONG «l TNO-LEAnINCj EoGE<*«//) 

PRINT T 

7 Ff)RMAT(l4X,^*S(l)*»7A»*PE(I)*»7X»«RES(T)**Rx**UU0r ( I H» , 7X , *(JUr.TR ( I) 

j*, 7X,«PRAT u ) »»/) 

UO Ij I=?*J 

PRINT 8»s ( I ) »PE u ) » reS ( I ) fODO r ( 1 ) «CnnTR 1 1 ) »PPAl ( I ) 

8 format (lr)X»El0.5»E10,R»3A*El2.5»3y,Fl?.i3»3xfF12,b,3x«El2,b) 

10 continue 
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DESCRIPTION OF OUTPUT 


The output for the perfect-gas code includes the flow conditions from 
each region of the flow-field model, the geometry of the shock waves and 
the expansion waves , and the heat-transfer distribution along the second 
wedge, i.e. , the wing leading-edge. The units for a particular parameter 
in any region will be the same as the free-stream parameter, unless other- 
wise noted. The output for the free-stream flow includes: 

FSMACH - free-stream Mach number 

2 

PI - free-stream static pressure (ibf/ft ) 

T1 - free-stream temperature (°R) 

3 

RHOl - free-stream density (slugs/ft ) 

PTl - free-stream stagnation pressiare (ibf/ft ) 

TT - stagnation temperatin’e of entire flow model (°R) 

RGAS - gas constant (ft^/sec^-°R) 

GAMMA - Y for "the free-stream and region 2 

QTSPR - reference heating rate (Btu/ft -sec) 

The following flow conditions are output for each region "I", where 1=2, 
3 , or 6 : 

MI - the Mach number of the region 

PI/Pl - the static pressure ratio 

TI/Tl - the temperature ratio 

RHOI/RHOl - the density ratio 

PCOEFI - the pressin?e coefficient 

PI - the static pressure 

TI - the temperature 

RHOI - the density 

PTI - the stagnation pressure 



DELTAJI - the change in flow-direction between two consecutive regions , 
in degrees. 

THETAJI - the shock wave angle between two consecutive regions , in 
degrees 

GAMMAI - the y of the region 

The output for the expansion fan region includes the following stagnation 
condition and intersection points : 

PTE - stagnation pressure at the edge of the boundary layer 

TTE - stagnation temperature at the edge of the boundary layer 

Intersection points include : 

INITIAL POINTS ON BOW SHOCK WAVE, i.e. , the origin of the coordinate 
system, or the nose. 

WING INTERSECTION POINT, i.e., intersection of the two wedges 
BOW-SHOCK: WING-SHOCK INTERSECTION POINT, i.e., the intersection of 
the shock of the first wedge with the shock of the second 
wedge 

Next Five Points - the intersections of the centered expansion fan 
with the wing leading edge 

Last Five Points - the intersection of the reflected waves with the 
inboard shear layer 

Output for the flow in region 5E include : 

THETA - the flow direction , in degrees 
M5E - the Mach number 

P5E - the static pressure 

The Mach numbers and pressures of each of the other nine subregions of re- 
gions 4 and 5 are the output under the next two headings . 

Output listed for the interaction region between the left running and right 



running expansion fan waves includes : 

EXM4IW - the Mach number 
P4IW - the static pressure 
where "I” = A, B, C, and D. 

Output for the subregions between the inboard shear-layer and the "curved" 
shock wave include : 

THETA - shock wave angle 

X - x-coordinate point of the intersection of the shock wave 

and the left-running expansion wave 

Y - y-coordinate point of the intersection of the shock wave and 

the left-running expansion wave 
PR - pressure ratio with respect to region 1 

DR - density ratio with respect to region 1 

TR - temperature ratio with respect to region 1 

M - Mach number 

The output for the flow parameters and the heat- transfer distribution along 
with wing leading edge include: 

S(I) - distance along wing leading edge 

PE(I) - pressure at the edge of the boundary layer 

RES(I) - Reynolds number 
QDOT(I) - heat transfer rate 

QDOTR(I) - heat transfer rate ratio with reference heat transfer rate 
PRAT(I) - pressure ratio with respect to region 1 
where the (I) index refers to the station for which the calculations are made. 



CASE • 1 


FREE-STHEAM flow conditions 


fsmach 

PI 

(LSF/SQFT) 

T1 

<RaNK1NE) 

RHOl 

(SLUqS/CUFT) 

PTl 

TT 

31000E*ol 

9,29674E-01 

4.90569E*02 

1.1046oC-06 

2.43547E*05 

1.732792*04 


gas constants 

R6AS ■ 1.71bO0E«03 GAMMA ■ 1.40000E*00 OTSPR * 5.G927SE*02 
FLOW conditions IN REGION ? 

m2 P2/P1 T2/T1 RH02/RH01 PCOEF2 

1.02AAOE*ol 4.03662E+00 1.606*AE*00 2.51277E*nO 2.52784E-fl2 

P2 s 3.75355E*00 T2 • 7,880712*02 RHo2 ■ 2.77562E-06 PT2 b l.87l00E*05 
0ELTA12 B 5,OOOO0E*OO THETA12 ■ 8,3311bE*00 

FLOW conditions in REGION 3 


m3 

P3/P1 

T3/T1 

RH03/RH01 

PCOEF3 

3.39620E*00 

1.376092*02 

l.06815E*01 

1.28829E*ni 

1.137202*00 


P3 ■ 1.279S9E,02 T3 b 5.240022*03 RH03 a 1,423052-05 PT3 a 5.08890E*02 
DELTA23 a 2.50000E*01 THETA23 ■ 3.193612*01 GAMMA3 a 1.400002*00 

FLOW conditions in region 6 


M6 

P6/P1 

T6/T1 

RH06/RH01 

PC0EF6 

2.83044E*00 

7.562032*01 

1.357352*01 

5.571192*00 

6.21179E- 


p6 a 7,03l74|«01 


T6 a 6.6S873E*03 


RHo6 a 6.15396E.A6 


PT6 a 2.05891E.02 



DELTA16 ■ 3«00000E*01 


THEtAl6 • 3,79703E*01 


QAMMA6 ■ lt40000F«00 


POINTS Of interest in the plow FIeLO and on the body (IN INCHES) 


PTE ■ 8.Al422E*03 TTE • 1.73279E*04 

initial POINTS ON BOW SHOCKWAVE 

X ■ 0* Y ■ 0» 

WING intersection point 

X B 9,00000E*00 Y ■ 8<I0000e»01 

bow shock-wing shock intersection point 
X ■ 9t839o9E*00 Y > 1.A40S3 e*00 

FIRST EXPANSION WAVE WING INTERSECTION POINT 
X » l,02588E*0l Y ■ l.53677E*00 

LINE 48 WING intersection POINT 
X ■ 1.02787E*01 Y ■ l.5«827E*00 

LINE 4C WING intersection POINT 
X a 1.03007E*01 Y ■ 1.56097E*00 

LINE 4D WING intersection POINT 
X a 1.03252E*0I Y ■ 1>S7508E«00 

line 4e wing intersection point 
X » 1*03526E*01 ^ ■ l*S9093E*00 

line sa shear intersection point 

X ■ 1.06104E*OI Y ■ 1>9A68SE*00 

line SB shear intersection point 

X a li06574E*0l Y ■ 1.97692E*00 

line sc shear intersection point 

X B l.OTOaOE^Ol Y ■ 2.00B46E*00 

line 5d shear intersection point 

X a lt07626E»0l Y ■ 2.04l66E*00 

line 5e shear intersection point 
X B 1.08216E*0l Y B 2,07662E*00 

FLOW IN REGION SE FROM EXPANSION 
THETA B 3,00000E*01 MSe b 3,82344e*00 pSE a 7*03174E*81 

MACH numbers in expansion Fan 

M4A a 3,43600E*00 mAB a 3,47627E*00 MAC m 3.SI700e*00 

MAO B 3.SSB90e*00 MAE ■ 3,6013IE«00 

MSA B 3.6AA21E*00 MSB m 3.68763E*00 a 3.73228E*00 

MSD ■ 3*777A8E*00 

pressures In the expansion fan 

P4A a I.208S8E*02 p4B a 1.1A106E*02 PAC a 1.07689E«02 

PAD a 1*01A9AE»02 PAE a 9,S618AE*01 

PSA a 9.00A70E*0l P5B a 8,A787lE*01 P5C a 7.96868E*01 

PSD a 7«488IAE*01 



Flow alonO the mIng In the expansion region 


MACH numbers 

MAAW ■ 3.4.7627E*00 MABW b 3,S59i0E«00 H4CW b 3.64A08e«00 M40W 

PRESSURES 

P4AW B 1,j4106E*02 P4BW b 1,01A67E»02 PACW b 9,Oo833E40I P.AOW 

FLOW IN THE SHOCKED REGION OF EXPANSION FAN 


REGION AES 

THETA B 4.2A924E401 
PR a 9.il872E*01 

REGION 5AS 

THETA B 4.15671E401 
PR B 8,797I5E*0l 

REGION BBS 

theta a 4.06615E*01 
PR B e,48364E401 

REGION 5CS 

THETA B 3.9.7480E401 
PR B 8.1689oE*01 

REGION SDS 

THETA B 3.8B555E401 
PR B 7.86319E«0I 


X B 1.14125E*01 
DR B 5t63987E«00 


X B 1.I6193E*01 
OR a 5.6275SE*00 


X B i.I8523E*01 
DR B 5.61A69E400 


X a 1.21U9E«0I 
DR a 5.60086E400 


X B 1.2Al07E*0l 
OR B S.S8644E400 


Y B 2.8822AE400 
TR B 1«61683E«01 


Y a 3.065S9E400 
TR a 1 i56323E*6I 


Y B 3,26571e* 00 
TR a ltS1097E401 


Y B 3.48414E400 
TR B 1.4SeSlE*0l 


Y a 3«72243E*00 
TR a 1.40755E«01 


3.73268E*00 


7.96434E40I 


B 2,43377e*00 


B 2.509S4E400 


a 2.58621E400 


a 2,66626E*00 


a 2.74724E«00 



FLOW PARAMETERS ALONG WING-LEADING EDGE 


S(l) PE(I) RESii) 


OOOT(I) QOOTRII) pRATd) 


.00400 127.95674 
.00600 127.95674 
.01200 127.95674 
.01600 127.95674 
.02000 127.95674 
.02400 127.95674 
.02600 127.95674 
.03200 127.95674 
.03600 127.95674 
.04000 127.95674 
.04400 127.95674 
.04600 127.95674 
.05200 127.95674 
.05600 127.95674 
.06000 127.95674 
.06400 127,95674 
.06600 127.95674 
.07200 127.95674 
.07600 127.95674 
.OSOOO 127.95674 
.06400 127.95674 
.06600 127.95674 
.09200 127.95674 
.09600 127.95674 
.10000 127.95674 
.10400 127.95674 
.10600 127.95674 
.11200 127.95674 
.11600 127.95674 
.12000 127.95674 
.12104 127.95674 
.12122 114.10592 
.12140 114.10592 
.12156 114.10592 
.12176 114,10592 
.12194 114.10592 
.12212 114.10592 
.12230 114.10592 
.12246 114.10592 
.12266 114.10592 
.12284 114,10592 
.12302 114.10592 
.12320 101.46675 
.12338 101.46675 
.12357 101.46675 
.12375 101.46675 
.12393 101.46675 
.12411 101.46675 
.12429 101.46675 
.12447 101.46675 
.12465 101.46675 
.12463 101.46675 
.12501 101.46675 
.12519 90.06329 

.12537 90.06329 

.12555 90.06329 


4.32403E«02 

8.64607E«02 

1.2972iE*03 

1.7296iE*03 

2.16202E*03 

2.59442E*03 

3.02662E«03 

3.45923E«03 

3.69163C«03 

4.32403E«03 

4.75644E«03 

5.18684E*03 

5.62124E*03 

6*0536SE«03 

6.46605E«03 

6.9i84SE«03 

7.35086E*03 

7.78326E*03 

8.2l566E«03 

8.64807E*03 

9.06047E*03 

9.51287E*03 

9.9452|E«03 

1.03777E»04 

1.0810iE*04 

1.1242SE*04 

1.16749E*04 

1.21073E«04 

1.25397E«04 

1.2972iE*04 

1.30843E*04 

1.23735E*04 

1.23919E*04 

1.24103E*04 

1.24268E*04 

1.24472E*04 

1.24656E*04 

1.2484oE«04 

1.2502SE*04 

1.25209E*04 

1.25393E*04 

1.2557SE*04 

1.16532E»04 

1.16706E*04 

1.18866e«04 

1.19053E«04 

1.19227E*04 

1.1940iE«04 

1.19574E«04 

i.l9745E*04 

1.19922E*04 

1.20096E*04 

1.20269E«04 

1.13369E*04 

1.13532E«04 

1.13696E*04 


3.93434E*02 

2.76200E«02 

2.27149E+02 

1.96717E*02 

1.75949E»02 

1.60619E«02 

1.48704E«02 

1.39100E*02 

1.31145E*02 

1.24415E*02 

1.18625E«02 

1.1357SE*02 

1.09119E*02 

1.05150E«02 

1.0l584E«02 

9.83584E*01 

9.54217E«01 

9.27332E«01 

9.02599E«01 

8.79745E»01 

8.58543E*01 

8.38804E«01 

8.20366E«01 

8.03093E«01 

7.86867E*01 

7.71587E*01 

7.57164E*01 

7.43520E*01 

7.30588E*01 

7.16308E*01 

7.15222E*01 

6.77579E*01 

6.77075E*01 

6.76572E«01 

6.7607oE*01 

8.75570E«01 

6.75070E«01 

6.74572E«01 

6.74074E*01 

6.73576E*01 

A.73o63E«01 

6.72589E«01 

8.36267E«01 

6.35601E«01 

6.35336E*01 

6.34873E«01 

6.34410E*01 

6.33949E«01 

6.33488E«01 

6.33026E*01 

A.3257oE«01 

6.32112E*01 

6.31656E*01 

S.96954E«01 

5.96524E«01 

S.96095E«01 


6.67658E-01 

4.72105E-01 

3.65472E-01 

3.33829E-01 

2.98586E-01 

2.72570E-01 

2.523S1E-01 

2.36053E-01 

2.22553E-01 

2.11132E-01 

2.01306E-01 

l#92736E-0i 

1.6S175E-01 

1«78439E-01 

1.72388E-01 

1.66914E-01 

1.61931E-01 

1.57366E.01 

1.53171E-01 

1.49293E-01 

1.4S695E-01 

1.42345E-01 

1.39216E-01 

1.36285E-01 

1»33532E-01 

1.30938E-01 

1.28491E-01 

1.26175E-01 

1.23981E-01 

1.21897E-01 

1.2l373E-0i 

1.14985E-01 

1.14900E-01 

1.14814E-01 

1.14729E-01 

1.14644E-01 

1.14560E.01 

l.l4475E«0i 

1.14391E-01 

1.14306E-01 

1.14222E.01 

1.14138E-01 

1.0797SE-0i 

1.07696E.01 

1.07817E-01 

1.07738E-01 

1.07659E-01 

1.07561E-01 

1.07503E-01 

U07425E-0i 

1.07347E-01 

1.07270E-01 

1.07192E-01 

1.01303E-01 

1.01230E-01 

1.01157E-01 


I.37609E*02 

1.37609E*02 

1,376o9E*02 

i.37609E*02 

1.37609E402 

i.37809E*02 

i.37609E*02 

i.37<>09E*02 

i.37809E402 

1.37609E*02 

1.37609E*02 

i.37609E*02 

I. 37609E*02 

J. 37609E*02 
1.37609E*02 
1.37609E*02 
1.37609E*02 
i.37609E*02 
i.37609E*02 
i.37609E*02 
i.37609E*02 
1 .37609E*02 
1 .37609E*02 
i.37609E*02 
i.37609E402 
1.37609E402 
i.37609E*02 
i.37609E402 
i.37609E*02 
1.37609E*02 
I.37609E*02 

i.227UE*« 

1.22711E402 
1.22711E*02 
1.22711E402 
1. 227116*02 
1.227116*02 
1.227116*02 
1.227llE*02 
1.227116*02 
i.227llE*02 
1.227llE*02 
1,6’1196402 
1.09119E402 
1.69119E*02 
1.69119E*02 
i.69ll9E*02 
1.09lT9E*02 
1.69119E402 
1.09li9E*O2 
1.09119E*02 
i.0’1196*02 
1.09119E402 
9.66554E«01 
4.66S54E«01 
9.66554E*01 



.12573 

90.06329 

1.138592*04 

5.956662*01 

1*01085E«01 

9.6-554E*01 

.12591 

90.06329 

1.U023E*04 

5.952392*01 

1*010122«01 

9.68554C*01 

.12609 

90.06329 

1 . 141862*04 

5.948132*01 

lt009402-0i 

9.685542«01 

.12627 

90.06329 

1.143502*04 

5.943882*01 

1.008682-01 

4.685542*01 

.12645 

90.06329 

1.145132*04 

5.939632*01 

1*007962-01 

9.685542*01 

.12663 

90.06329 

1.146772*04 

5.935402*01 

1*007242-01 

9.685542*01 

.12681 

90.06329 

1.148402*04 

5.931172*01 

1*006522-01 

9.685542*01 

.12700 

90.06329 

1.150042*04 

5.926952*01 

1*005802-01 

9.685542*01 

.12718 

90.06329 

1.151672*04 

5.922742*01 

1*005092-01 

9.685542*01 

.12736 

90.06329 

1.153312*04 

5.918542*01 

1*004382-01 

9.685542*01 

.12754 

79.64335 

1.084642*04 

5.582862*01 

9*474122-02 

8.584962*01 

.12772 

79.64335 

1.086112*04 

5.578912*01 

9*467422-02 

6.584962*01 

.12790 

79.64335 

1.087712*04 

5.574972*01 

9*460742-02 

8.584962*01 

.12808 

79.64335 

1.089252*04 

5.571042*01 

9*454072-02 

8.564962*01 

.12826 

79.64335 

1.090782*04 

5.567122*01 

9*447412-02 

8.564962*01 

.12844 

79.64335 

1.092322*04 

5.563212*01 

9*440772-02 

8.584962*01 

.12862 

79.64335 

1.093852*04 

5.559302*01 

9*434142-02 

8.564962*01 

.12880 

79.64335 

1.095392*04 

5.555402*01 

9*427532-02 

8.564962*01 

.12898 

79.64335 

1.096922*04 

5.551512*01 

9.420932-02 

8.564962*01 

.12916 

79.64335 

1.098462*04 

5.547632*01 

9.414342-02 

8.564962*01 

.12934 

79.64335 

1.100002*04 

5.543762*01 

9.407772-02 

8.564962*01 

.12952 

79.64335 

1.101532*04 

5.539902*01 

9.401212-02 

8.564962*01 

.12970 

79.64335 

1.103072*04 

5.536042*01 

9.394662-02 

8.564962*01 

.12988 

79.64335 

1.104602*04 

5.532192*01 

9.388132-02 

8.564962*01 

.13006 

79.64335 

1.106142*04 

5.528352*01 

9*381622-02 

8.564962*01 

.13025 

70.31739 

1.039072*04 

5.210232*01 

8*841762-02 

7.582032*01 

.13425 

70.31739 

1.070982*04 

5.132022*01 

8*709042-02 

7.582032*01 

.13825 

70.31739 

1.102892*04 

5.057232*01 

8*582122-02 

7.562032*01 

.14225 

7 .31739 

1.134862*04 

4.985612*01 

8*460592-02 

7.582032*01 

.14625 

70.31739 

1.166712*04 

4.916962*01 

8*344092-02 

7.582032*01 

.15025 

70.31739 

1.198622*04 

4.851072*01 

8*232272-02 

7.582032*01 

.15425 

70.31739 

1.230532*04 

4.787752*01 

8*124822-02 

7.582032*01 

.15825 

70.31739 

1.262442*04 

4.726852*01 

8*021482-02 

7.582032*01 

.16225 

70.31739 

1.294352*04 

4.668222*01 

7*921982-02 

7.582032*01 

.16625 

70.31739 

1.326272*04 

4.611722*01 

7*826102-02 

7.582032*01 

.17025 

70.31739 

1.358182*04 

4.557222*01 

7*733612-02 

7.562032*01 

.17425 

70.31739 

1.390092*04 

4.504612*01 

7*644332-02 

7.582032*01 

.17825 

70.31739 

1.622002*04 

4.453782*01 

7.558072-02 

7.582032*01 

.18225 

70.31739 

1.4539j2*04 

4.404632*01 

7.474662-02 

7.562032*01 

.18625 

70.31739 

1.485822*04 

4.357072*01 

7.393962-02 

7.562032*01 

.19025 

70.31739 

1.517732*04 

4.311032*01 

7.315822-02 

7.582032*01 

.19425 

70.31739 

1.549642*04 

4.266412*01 

7.240102-02 

7.582032*01 

.19825 

70.31739 

1.581552*04 

4.223152*01 

7.166692-02 

7.562032*01 



APPENDIX B. - GENERAL DESCRIPTION OF REAL-GAS CODE 


As has been discussed previously, philsophically the real-gas code 
is approached in the same way as the perfect-gas code. As with the per- 
fect-gas code, the free-stream conditions and the model geometry serve as 
input. The real-gas code uses a Mollier gas table to find the flow-field 
conditions in place of the perfect gas relations of the perfect-gas code. 


Calculation of Flow-Field Conditions * The flow conditions of regions 2 , 

3, and 6 are calculated in subroutine RELGAS. The procedure of this sub- 
routine is to assume an initial value for the shock wave angle 0, and solve for 
the density as follows. 


p^tan 0 
^2 tan(0-6 ) 


(18) 


The velocity behind the shock wave can then be determined by using the 
conservation of mass equation: 


U 




2n 


(19) 


The pressure is then calculated from the conservation of normal momentum 
equation : 


"l " "An' - f2“2n' 


( 20 ) 


Then the enthalpy behind the shock wave is calculated using the conserva- 
tion of energy equation: 


^2 = 


h, + 0.5U, 

1 In 


- 0.5U 


2n 


( 21 ) 


Once the static pressure and the static enthalpy are known, the remaining 
flow properties behind the shock wave can be found by using the subroutine 
M0LIER. M0LIER, which incorporates the gas table, needs only two properties as 
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input to find the remaining flow field properties. The density found from 
M0LIER is compared to the density behind the shock wave calculated from 
equation (18). If the comparison is not acceptable (to within a prespecified 
tolerance), a half-interval iteration method is used to redetermine the shock 
wave angle . The procedure to calculate the flow condition behind the shock 
wave is then repeated until successive values of the densities converge. 

Calculation of Stagnation Conditions . The stagnation conditions are calcu- 
lated after the flow conditions in regions 2, 3, and 6 are calculated, by 
first calculating the stagnation enthalpy by the following equation: 

H = hj_ + 0.5 (22) 

Using the stagnation enthalpy and the entropy of regions 1, 2, 3, and 6, 
the stagnation conditions for an isentropic deceleration of the local flow 
are found from M0LIER. 


Calculation of Expansion Fan Flow . To calculate the expansion fan flow, the 
flow is assumed to expand isentropically from region 3 to region 5E. Thus, 
the entropy of region 3 defines the entropy of the expansion fan. Therefore, 
the two properties serving as input for the M0LIER subroutine , are the local 
static enthalpy and the entropy. 

As before, the streamlines for subregion 5E through region 6 are parallel 
and straight for the assumed flow model. Therefore, 



(23) 


and the static pressure of subregion 5E is equal to the static pressure of 
region 6. The enthalpy of region 5E is found by using the static pressure of 
region 6 and the entropy of region 3. To calculate the enthalpy of each of 
the other nine regions in the expansion fan, the difference between the 
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enthalpies of region 3 and 5E is divided by ten to get ten equal increments 
The enthalpy of region 4A is, therefore, the enthalpy of region 3 plus the 
above increment. The enthalpies of the remaining subregions are gotten by 
appropriately incrementing the previous enthalpy . 

As a check to the above procedure the Prandtl-Meyer angle for the sub- 
region of interest was calculated using the equation: 

1 f f2(H-h) n d(H-h) ^ 

Reg. 3 

When crossing a right-running wave, the change in the Prandtl-Meyer 
angle is equal to the change in the flow direction, i.e. , 

dv = de (25a) 

When crossing the reflected left-running waves, 

dv = -de (25b) 


Thus, since the flow in region 3 is parallel to the wall, the sum of the 
change in v through region 4 should equal to the sum of the change in v 
through region 5. For all cases computed thus far, the difference between 
the v-sums has been small. Thus, the net change in flow angle has been 
essentially zero, as it should. The intersection points of the shock waves 
and expansion fan waves are calculated the same way as in the perfect gas 
code using the subroutine INTRST . 

The flow properties at the wall account for the waves of the centered 
ejqiansion fan, i.e., region 4, and of the reflected left-running waves, i.e.. 


region 5 . As noted before , 

^4AW ■ ^^4A ■ ^®4AW " ®4A^ 


(26) 


SO that 


^4AW ~ ^4B 


(27a) 
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Similarly , 


'^4BW " ^4D 


(27b) 


^4CW 


5A 


(27c) 


and 


^4DW "^50 


(27d) 


Once the Prandtl-Meyer functions are known, the remaining properties are 
readily calculated, since flow in this region is isentropic. 


Wing Leading-Edge Heat-Transfer Rate Calculations. After the inviscid flow- 
field conditions have all been calculated, the heat-transfer rate along the 
"wing leading-edge" is calculated. The procediare for calculating the heat- 
transfer is to first use the numerical routine as described in ref. 15 to 
set up the initial boundary-layer profile. 

This initial profile serves as an input into the NONSIMBL code (ref. 11). 
Several modifications were made to this numerical routine for a laminar 
bo\mdary-layer. One modification was to eliminate the need for user "exper- 
ience" in establishing the initial, input guesses for the wall values of the 
shear function, f"(0), and of the heat-transfer function, g'(0). Instead, 
since the temperature of the edge of the boundary layer and at the wall in 
region 3 were known from the previous subroutines , the Initial value of 


g'(0) was assumed to be 


g’(0) = 



0.6 


The initial guess for the shear function was assumed to be: f"(0) = 0.47, 
which seemed to provide reasonable results over the entire velocity range 
(in the absence of gas injection of the wall). 

Another modification was to incorporate into the routine a proced\jre 
for calculating alpha (i.e. , the coordinate transformation parameter). 
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The above procedure is designated subroutine PIGYBAK and transforms the 
distribution of the dimensionless velocity function F and the temperature 
ratio 6 (THETA) into the new coordinate using alpha (a). The transforma- 
tion is 

, -ari 
n = 1 - e 

where n is the newly transformed y-coordinate (ref. 11) and ti is the 
transformed y-coordinate using the standard Lees-Dorodnitsyn transformation. 
To calculate the appropriate value of a, the value of n where u = 0.99u^ 
in the similar solution for the first station was identified. Then, the 
value of a was calculated to be the value which makes the produce an equal 

to three at this point, i.e. , n = 0.95. 

This viscous-layer profile for a similar boundary-layer which is sub- 
ject to the flow conditions at the initial point in region 3 is used as an 
initial condition for the subroutine EJ0YCE, which solves the nonsimilar 
boundary layer. This subroutine is actually the numerical code described 
in ref. 11 with modifications. Modifications were made to this code so that 
one could obtain real gas values of density, thermal conductivity, viscosity, 
and specific heat of air. The M0LIER subroutine was used to calculate the 
density in the boimdary layer. The tabulated values of Hansen (ref. 9) for 
thermal conductivity, viscosity, and specific heat were input as a function 
of temperature. Then, the subroutine SPLNTRP was used to calciilate the 
transport properties at the desired temperatures. SPLNTRP is a curve fit- 
ting subroutine. The final output of EJ0YCE gives boundary layer profiles, 
heat-transfer rates , and other data for five points in each of the six 
regions and siobregions on the "wing leading-edge". 

Calculation of the Shock-Wave Angles and the Intersection Points of the 
Shock Wave with the Expansion Fan. As was the case for the perfect-gas 


code the flow direction in subregions 4E, 5A, 5B, 5C, and 5D of the expansion 
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fan are used as boundary conditions for computing the shock waves in this 
region. The subroutine RELGAS is used to generate the shock wave angle 
and the flow conditions in the subregions between the shock wave and the 
shear layer, i.e., subregion 4ES, 5AS, 5BS, 5CS, and 5DS. The intersection 
of the shock wave and the left-running expansion wave is computed using 
the subroutine INTRST using the same procedure as discussed in the perfect- 
gas section. 



Input Cards 


Card # 1 
Card # 2 


Card # 3 


Card # 4 


Card # 5 


Card # 6 
Card # 1 

Card # 8 
Card # 9 


NC - number of cases 

XO - x-location of "nose" (ft) 

YO - y-location of "nose" (ft) 

X30 - x-location of wedges intersection point (ft) 

Y30 - y-location of wedges intersection point (ft) 

U1 - free-stream velocity (ft/sec) 

2 

PI - free-stream pressure (Ibf/ft ) 

2 2 

HI - free-stream enthalpy (ft /sec ) 

DELTA - angle of the first wedge (^) 

LAMDAS - complement of sweep angle ( ) 

N0PT - option in M0LIER for the two input properties enthalpy 
and pressure ( ) 

TWALL - wall temperature (^R) 

RADIUS - radius of reference sphere (ft) 

PR - Prandtl number 

E - convergence criterion on boundary condition (typically .0005) 

DELT - step size in y-direction (typically .05) 

EPS - convergence criterion when variable step size is used. 

KEY - equals 1 if using variable step size (equals 0 for the 
fixed size routine) 

IN - number of y-points 

ALE - coordinate transformation parameter 

PIGYBK - if equals 1.0 edge properties are not read in as input 

MM - number of x-points for calculations 

N - number of y-points in boundary layer 

KK - equals 0 for two-dimensional flow, equals 1 for 
axi-symmetrical flow 

N0PRINT - if equals 0 output will be printed. 



Card # 10 


WMI - molecular weight of injectant 


WMS - molecular weight of stream (WMI = WMS for this routine) 


Card # 11-13 TEMP - temperature (°R), fifteen temperatures on three cards 

Card # 14-16 VISC - viscosity ) fifteen viscosities correspond- 

ft^ ing to above temperatures. 

BTU 

Card # 17-19 TC - thermal conductivity 5 fifteen thermal con- 

Tx sec ductivities correspond- 

ing to above temperatxires 


Card # 20-22 


CP - specific heat 


_MT^) 

'■slug-OR-’ 


fifteen specific heats correspond- 
ing to above temperatures 


Card # 23 equals Card # 3 if NC > 1 
Card # 24 equals Card # 4 if NC > 1 
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HROGKAM REALGAS aNPUT f OUTPUT ) 

UIMEWSION H**(50) ,P4(50) ,TA(5U) ,hA( 50) f04(5o) tAASuiSO) ♦Y(50) ♦ 
IunU( 50) */4(t>0) 

UIMENSIOM WAVE4 (SO) f AGL4(50) ♦ A4 ( 50 ) t Y4 ( bO ) ,THTA4(50) fU4(50) t 
1EXM4 (bO) tSiGUNU (SO) 

dimension T>iTS4(ll)» AS4(11) *YS4(11) , US4 ( 1 1 ) iPS4 ( 1 1 ) , TS^MU)* 
1HS4(U) *hS 4(11) *r,S4(ll) fZS4(U) *SS4(11) 

COMMON/KONSI/TV*A| l*pauius*pr 

C0MM0N/WlNG/PT3»TT3»U3*T3»P3f X30*Y30»X4»Y4»P4*Pb»T4»U4»T5t 
COMMON/CaSE/NCASp 
t:OHMO^/ShOC^/NF LllBl ISMUCK* JSMUCt^ 
keal L^*MDAs 
HEAD Z2i2* nU 
222 FORMAT (Ib) 

HEAD l0i»X0tY0fXl0»Y30 

101 Format ( 4El2.b) 

UO lUifiNC 

HEAD 200»Ul»Pl,Mi , DEL T A i LAMOAS,NO pT 
200 F0RMAT(£>El2.b,I5) 

HEAD 102. T«ALL» radius* PR 

102 FORMAT (3El2»b) 

C 

HlsHl/(32.i /b*77(J,0) 

Call MOLitR(Mi*pi * nopt*ti*zi*si*rhoi.gammah 
H lsHl*32.l7b*778,0 
C 

c *****CALCULmTING FLOl# CONDITIONS IN REGION ?««*** 

CALL KELGAS (UI»p1*RHOI*H1*OELTA,n0PT*U2*P2*RHO2*H2*T2*Z2»GAmMA2. 
1IHETA*S2) 

C 

c ****4CALCULATING FLOIn conditions in region 3***** 
DELTA3sLAMDAS-DElTA 

call HELGAS(U2»P?»«H02*H2*DELTA3»NCPT*U3fP3*RHU3»H3»T3»Z3»OAMMA3» 
11 hETA3*S3) 

C 

C *****CALCULATIN6 FLOW CONDITIONS IN REGION f,*<»*** 

DELTAbsLAMDAS 

NFLUBsi 

CALL REL6AS(U1 »Pi »RHC1 f HI fDELTA6«N0PT*Ub*P6«RHub* H6«T6«Z6*GAMMAb* 
irHETAb*S6) 

••♦♦•computing Stagnation conditions***** 

MSTAGsHl*0.b^Ul*lJl 
HSTAGaHSTAb/ (32 . ] 7b^778 ,0 ) 

Call molier(hstar*pti»3*tti*zti*si*rti»gti) 

Call MOLI£R(HSTAg,PT2»3*TT2*ZT2*S2*RT2*GT2) 

Call M0LlER(HSTAr,,PT3f3*TT3.ZT3*S3.RT3fGT3) 

CALL M0LlER<HSTAG*PT6*3*TTb»ZTb*Sb*RT6*GT6) 

C 

PRAT2l»P2/Pi 

PrAt31»P3/PI 

pRATbi«p6/pl 

TRAT21»T2/T1 

TRAT3UT3/TI 

TRAT6l»Tb/Tl 

URAT21«Rh02/RH01 

URATJURH03/RH01 

DrAT61«Rm06/RH01 

c 

Ml«Ml/(32*l7b*77a.O) 



non nnn n n non 


H2«M2/(32.1/6*77fl,0) 

H3aH3/(32«W6*77fl,0) 

M6»H6/ (32»l76*77fl»0) 

**#**coMpUTiNG Expansion fan flow***** 

IF( ISHOCK.EO.l) Go TO 703 

CALL MOLlER(H5E*p6,l,T5E*Z5EfS3*RH05E,65E) 

A3SQaP3*GAMi1A3/RH03 
A5ESQ*P6*G3E/RH0RE 
A3SQaA3SQ/ (32,17^*778,0) 

A5ESQ»a5ES&)/ (32.) 76*7 78.0) 

£TAlaHSTAG-H3 
EtAnsHST AG-HSe 

Nall 

Y(l)sSQRT (2.0 *ETa1/A3SU-1.0)/ (2.0*ETA1) 

K (N) =SQRT (2.0«ETaN/A5ESQ-1.0) / (2.0«ETAN) 

MsN-1 

MSTEPs(ETAN-ETA1)/M 

ETAsETaI 

00 98 I=2»M 

ETAsETA^HSTEP 

M4(I)sHStAG-£TA 

CALL MOLIEr(HA(I) ,P4(I) ,3,T4(I) ,Z4(I) ,S3,R4(I) »64(I) ) 

A4SQ ( I ) *P4 ( I ) *G4 ( I ) /R4 ( I ) 

A4S0( I) sA4SW( I)/ (32,176*778,0) 

Y(I)=SUHT(2.U*ETA/A4SO(I)-1,0)/(2,0*ETA) 

98 Continue 

SIGONU(1)=0.0 

bIGDNU(2)sslGDNU(l) ♦(¥ (1) ♦Y(2) )*HSTEP/2*0 
UO 97 jx3»N»2 

5IGDNU( J)=S1GDNU( J-2) ♦ (Y( J-2) ♦4. 0*Y ( J-1) *Y ( J) )*HSTEP/3,0 
IF ( (J*l) .GT.M) Go to 97 

SlGDNU(j4l)asiGUNU(J-l) ♦ (Y(J-l)*4,0*Y(J)*Y(J*l) )*HSTEP/3,U 
97 CONTINUE 

00 9U jPs2,N 

ONU ( JP) =SlifUNU( JP)-SIGUNU ( JP-1) 

90 CONTINUE 

0NU4aUNU (2) ♦ONU (?) *0NU (4) 4DNU (5) ♦ONU (6> 

ONUSaUNU ( 7 ) 40 nU ( h ) ♦UNO (9 ) 40NU ( 10) 4CNU (11 ) 

UNUDIF 3 ONU 4 -UNUS 

*#***FINDING VELOCITIES AND MACH NUMBERS IN THE EXPANSION FAN***** 

HSTAGaHSTAG*32. 176*778,0 
00 9b Ka2tM 

H4 (K)»h 4 (K) *32. 176*778.0 
U4 (K) aSORT (2,0 * (hSTAG-H4 (K) ) ) 

A4SQ(K)sA4SW(K) *32.176*778,0 
EXM4(K)=ij4(N)/SQrT(A4S0(K) ) 

H4(K)aH4(K)/ (32,176*778,0) 

96 continue 

H5EaM5E*32. 176*778.0 
A5ES0aA5ESQ*32. 176*778,0 
U4(N)aSQRT(2,0*(HSTAG-H5E) ) 

EXM4(N) =U4 (N) /SQPT(A5£SQ) 

HSEaHSP/ (32. 176*778,0) 

*****FIN0ING intersection POINTS OF SHOCK WAVES AND EXPANSION FAN***** 
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lHtTA=THF.TA/b7.2'56 
lMtTA3aTHETMj/57.296 
LAMDaS*LAMUMi,/5 7.2'Sb 
utLTrisDEL U/57»2<^6 
UfcLTAb=UtLTMb/57.296 
C 

«(3L3 = 7HETA3 + UELT A 

tr.LL INFi^ST lXn»V ♦ T rit T m » a(5L3 » a3U ♦ y3 0 ♦ X I »YI) 

C 

A3SQ = A3bO*32. I 76^^778, i) 

aaVE 3 = A 5 J ( 1 . n / (ii 3 /S(;h r ( A 35 Q) ) ) 

A^,L4AsLA^'!DAb“'•'A Vp 3 

^ALL INTpSI lAl » Yt ♦ AbLAA»LAMl)A8*AJo»Y30*X4A» Y4A) 

C 

tXMA ( J ) ( ^-3SU) 

A4 (2J sX4A 

14 ( 2 ) =Y4c 
AGL4 ( 2 ) -aGL 4 A 
t-L = l 

^Vfvi= ((M* 1 ) /2 
UU 9b L=2»M 

•*AV£4 (L) sASl.'J ( I • /EXN4 (L) ) 

IF {N.^i-U 

92 A(-L4 (L+n =LAi'10AS*SI0aiviU (L) -Wa 'Ea (U 

I- ALL INTpST ( XT ♦ Y r » AbL‘+ (L + 1 ) » L AMUAS » X 30 * Y30 »X4 (l*1)»Y4(l+1)) 

u() TO 9b 

94 lL=LL*1 

I mTAa (L) = I H I m 4 (L- 1) -LNtJ (L ) 

A'iL4lL*-l)=TnTA4(i )+wAVt4(L) 

LaLL INTRbTiX4(L) ♦ Y4 ( L ) » TH T A4 < L ) » AGL4 ( L 1 ) , X4 ( LL ) » Y4 ( LL ) ♦ A4 ( L ♦ 1 ) ♦ 
lr4(L+1)) 

vr-o TU 9b 

93 i H r A4 ( L ^ = LAriL)Ab* I GU\U ( L) 

LL=LL*1 

AI.-L4 (L + 1 ) =Ti^TA4 (I )+aAVL4(L) 

LaLL INrpST(Al,YT,TrlTA4(L) »A&L4(L*1) » X4 ( LL ) < Y4 ( LL ) , X 4 ( L ♦ 1 ) ♦ Y 4 ( L ♦ 1 ) 
1 > 

95 LOntinuE 
C 

c**<hh»calculating FLLk»f coNoniONs Shocked hegions 4t anc 

c 

i-sl 

t'*FLUb=2 

Hi = H]*32.l7o<»778,o 
UO 77 I=^♦M 

I hTA 4 { I ) =T n I A4 ( 1 ) *57. 296 

call HELOAStUl»Pi THTA4 (I) ,NOPl fUS4 ( D *PS4( 1) ,RS4(I) ♦ 

1»)S4U) » lS4(i) ,ZSi (I) »GS4(I) »THTS4(I) *SS4(I) ) 

ImTSa (I)=lHl54(n/b7.296 
LsL*l 

IF(I.GT.F) go to 75 

CALL INl«STtXl»YT»THTS4(i) , AGL4 (I+l)*xA(L)»Y4(L) «XS4 (I*l)»YS4(i+i) 
1 ) 

GO To 76 

75 call INTpST US4 ( r ) » YS*^ ( 1 ) ♦ THTS4 ( I ) tAGL^ (I + l)»A4(L)tYA(L) ,AS4 (!♦!)♦ 
] Y S4 ( i ♦ ] ) ) 

76 lHrS4(I)sTHIS4(I)*57,296 

77 CONTINUE 

MlsHI/(32,l 76*77 h,0) 



lHETA3-Th£rAJ4»!37.?9b 
LAM0AsiLAMUAS«57 , ?96 
U&LTA sUI^LTAD*5 7.?9fc 
ueLTAs0tL.TA<»57*c'J6 
lH£TM = THfcTA*b7.i?96 
7CJ CuNTiNUt 
i^kINT 100 
100 r ( ihi ) 

KWINT -iOOf U 

300 K)WMAT (33X,*CASt 
f^klNt 20] 

2nl roHMMT(l6Af«(-ReE-STKEAfl FLOW CCiND j T I OnS* * /// ) 

PRINT 202 

202 KuRM«r‘U^»" U1 «♦ I 1X,*P1#, 14X»»T1*» i2X,«RHUt*» 12X»*pTl*f 12X,<»TT 
1 )*) 

PRlNl 217 

21? t- OKMAT^26X,<* (I.BK/'SUh T) (HANK INF) « tax, <» ( SLUbS/CUF T) * » / ) 

PRINT 203» ul, Pl» Tit HH(]1, PTlt TTl 

203 PoRMmT (l'iXtt-l2.b»3XtEl2,b,4XtEl2.5t3X,£i2,5tAAtci2.b,4XfEi2#bt//) 
print 20a* 'jAivtPAi, Hit /I 

2r<* PORMaT (1 IXt^oAPM.*, I =-»ttl2.5tSXt«Hl =» , £ 1 2 , 5 t bX t *Z 1 s* t E 12 • b ♦ /// ) 
PRINT 2 Ca 

206 P ORMA T ' lhX» A(.|_o,^ CONCITIONS IN HEOIOM ?.**///) 

PRINT 207 

2a 7 P OHM at ( l5Xt*U?« t ’ 4Xt«P2/Pl<*t 1’3X»<»T2/T1** 14X»*HH02/Rpi01« t /) 

PRINT 20rft u2t P-AT2lt TRAT21f URAT21 
2C6 PORMAT UoXtLl2.b.6XtE12,btHX»£l2,'5t9Xtt. 12.5t//) 

PRINT 212»Pc;»T2t >H02tPT2 

2]2 PORPAT ( UXt*»P2 s^vtEl2tbtbXt*T2 =* t E 1 2 ,b »5X t *Ht^u2 s«tEl2.5»bXt*PT? 
l=*tE12.5t/> 

PRINT 205* l i 2t^':^t22tt3AMlVtA2 

2a 5 PORMAT (UXt*1 T2 » E 12 .5 t bX t « h2 s» t E T 2 .b t 5 x » *^ 2 ***£’ 2.5,bA»*GAMMA 

1 C s«t£l2.bt/) 
print 209» utLTA. IHEIA 

2''9 PORMaT (UXt*iJELT ■v =<» 1 1 12 . 5 1 5X * *THET A =*»£12.b) 

PRINT 21' 

210 PORMAT V//t ioX.*f.*LOW tONUlTIONS IN REGION 3*t///) 

PRINT 213 

21 3 PORMAT (lbXt*iJ3«ti 4Xt*PJ/Pl<»tlbXf*T3/Tl*»Ux»*HPi03/RH0l«t/) 

PRIM 208t U3t PwAT3lt IHAT31* UrAT3) 

PRINT 215»PJ»T3» ^h03»PT3 

215 P ORMaT ( 1 lXt*p3 s^t£12.btbXt<»TJ a* t E 1 2 ,b » 5X t *Wrtu3 a*tF 12«S»6Xt<»PT3 
l=*f £12.5t/) 

PRINT 220» IT3»H’5t23 

220 PORMAT (llXt<*TT3 s* t E 1 2 . 5 1 5 A t *H3 a« , F 1 2,5 1 5x » *X3 a*tEl2,5,/) 

PRINT 21B» R£LTA->t TPETAJt GAMMA3 

218 PORMAT (llXt^^OELTp a* , E 1 2 . 5 1 5X ♦ * THFT A » E 1 2 ,5 1 5X t »QAMma 3 a^»,El;?,5) 

PRINT 211 

211 PORMAT<///t lbXi*FLOw CONDITIONS IM RFGlUN 
IFdSHOCK.ER.DGo TO 70l 

print 214 

2l4 H)RMAT ( 15Xt<>U6*t l4Xt*Pb/Pl*t 15X»*T6/Tl<»t 14 x»*HPi06/RH01<h/) 
print 20h» D6t PpATblt THATblt UHaT 61 
PRINT 21b»Po»T6f.-HObtPTb 

216 PORMAT (llXt*P6 att,ei2.5tbXt*T6 a» ,£ i 2 ,5 *5X t *RHOb a« t F 12 ,5 *5X t *pTb 
l*<»*El2,5t/) 

PRINT 22l» 1TbtHP,,Zb 

221 POKMATdlXt^^TTb -*tEl2.5tbXt*H6 a* ,F 1 2« S *5x a#,E12.5,/) 

Print 219 * uELlAxt TnElAbt GAMMA6 

219 FORMAT ( llXt*UELT, 5 a«,Ei2,St5X»«ThETA x**El2t5*5Xt*GAMMA6 =**£12, 5) 
PRINT 301 

301 PORMAT (///t loXt*FLOw CONDITIONS IN THE EXPANSION FAN*t///) 



UY (1 ) sY (?) 

UY(2)sY(3) 

UY (3) s- ( 1 ./Ai) « ( fY (1) ♦X2)*Y (3) ) 

UY(4)=Y(5) 

UY(5) =- ( 1 • / I A3« at ) ) « ( ( Y ( 1) ♦A3*X?) *Y (5) ♦ ( (X J*Y (3) **2)-Y ( i) *Y (2) «Y(3 
1 ) ) *A4) 
rUHi'i 
tiMO 

■ " ■ sJhRUUTInE -'•OMlN f N*DElT» r » 0 V»P»TE ,ER,YPH»OEP»£^S»DPS,^VE»Kt YtUERIV, 

liiOX) 

C— — — (iENEKAL I'IUi'IERICaL IMfc&KATIUN WITH SINtiLE-SfEP EHROH (\NaLY 1>IS 
UOUBLE KHElriSION UEPfSiXf TwOtUPbavE 
l^IMEivblO! , A 14) (5) »Cy (N) »P (N) f TE (N) »EP (N) 

UImEinSIUN ykk(K»a ) tUEPCN) »OPSaVE (N) 

C- 

C— — 

ExTtf^NAL UEnlY 

C — --------- 1„>ATA 

ou Tu (lu»li>»?v)) TNOX 
10 bIX = 

XEXO = O.ODo 

A( J ) S -V.Oi^U/SIv 
M(2) = 3/.Ouu/f5l» 

«(3) s -bi'.uJO/SrX 
A (4) s 5b»0UU/SI< 
t3(l) s i,UUu/$IX 
b(2) s -b.0t^')/SI» 
b(J) = I'^.Oi-'O/SIv 

0(4) zi -A ( 1 ) 

RATIO s 1V.V//270.0 
aTX = 6.0UU 

I WO s 2.000 
('.1 s 4 
I'^2 = .! 

•"13 s a 
•'14 a 3 
hi-TUkn 


SE T OH FGH INTEGRA T I ON) 

C— — tNTHY reset 


15 mssign 2ot;o rn ipls 

NOP a 0 
XOUNT a f) 

OELbYft a OElT/SIy. 

OELBY« a DEL T /Two 
00 lb jalfN 
UPSAVE(J) a OEP(J) 

1 b UV ( J) a L)EH ( J) 

Call OERI'' ( ' »DV»P,6»N) 

17 continue 

00 18 Oal*N 
rpR(0»Ml) a p(J) 

18 ER(J) a ZERO 

19 continue 

MFTUriN 

C— — — ENTHY-POInt for one numerical ImTEGHaTIOn step 

2 U V a r 

T * Y ♦ OEli 
NOP a KUP ♦ 1 
00 To 1 Pl5 

C— — RUnGE-NUITA PdoCEDuHE 
2000 VV a V ♦ 0ELHY2 
00 25U Jal»N 



250 u\/(J) = oth>(J) ♦ YHW ( J,Mi) «DELBY2 
^.iLL UEKiv (vv,UV.P»3*n) 

25 i CutMTlNUt 
Uu 2b 

260 - pEB^j) ♦ P(J)*PELBY2 
^/VLL UEBIV ( vv»0V,T£»4»r^) 

261 CuNTiiMUt 

op 2?u Jsifi'f 

uv (J) = pEpij) ♦ T£(j)*UELT 
270 IF (j) = 2.0*(TE ( I) ♦H(J) ) 

Pt-RIV ( I ,Pv»f ,5t^) 

?71 COiMTlNUE 

u o 2 b 0 J — 1 * ' '< 

uEP(JI s UEH(J) * UELbY6*(P(J)*TE(J)4Y^^<J»'^l^> 

2ho Uv(J) s ijEplJ) 

i-ALL Ut-Hiv ( I ♦L»v*p»5»N) 

2Pl t-ONTlNljE 

R(.) Jslfl'l 

290 tpR(j,M4) s P(j) 

2m1 ConTIi'-uE 

r\0UNT = kUUNT ♦ ' 

C- CHECl^ ImE ^■UMbER OF INTEGRATION STEPS NAPE BY R-K 

IF (KGUNT ,LT.i) GO TO 5000 
292 MSSIOM lOOu TO IpLS 
^0 TO 5000 
So 00 MU = M4 
'•lA s m3 

m 3 = m 2 
M2 s FI 
Ml = MO 

IF (KEYlEO*Oi KFTUPN 
IF (KuP.LE. Jl RF. TUwN 
i^INCsC 
p() 50 0 

IF (ER ( J) .GT.EPS) go to 515 
IF (Eh ( J) *iOuOO.O,L1 .EPS) MNC*MNC*1 
boo CONTINUE 

IF (KING. £0. IV) GO TO 5b0 
up 51 C J = l»iv 
uPSAvt ( J) = UFP ( i) 

510 CONTINUE 

return 

515 continue 

1F(KOP.Eo.4) go ro 5 ao 
1 S T - utLI 
520 PELT s 0,5po<»[)ELT 
PO 5io Jal.iv 
PEP(Jl s DP3AVE(J) 

530 continue 
oO TO 15 

5a 0 I a T " 4.0uO«DEi T 
oO TO 52 0 

550 PELT = 2,Ouo*OELt 
GO TO 15 

C— ...... mDAi'^S pREnlCTOH-CORHECTOH PROcEOUPE 

lOOo PO 710 Jsl 

IE(J) = H (3) *YPH ( J.Ml ) ♦ y(2)*YPR(J.M2) ♦ R (1 ) *YRR ( J.M3) 
UV(J) a PEP(J) ♦ 0ELT»(A(4)*YPR(J*P1 ) ♦a( 3)*YPK(J*M2) ) 
PV(J) a pV(J) ♦ nELT<MA (2) »YPH ( J*m3) +A ( 1 ) *YPR ( J*m4) ) 

C- SAVE PHtPICTcD VALUE 

710 Er(J) = DV ( J) 

CALL DERIV ( r.nv.p.i.N) 


PROCEDURE 



30^ 

302 M)KMAT 4«»//) 

HKINT 3U 

314 fOHMAT (5r:X,*VEL0rlTIES«»/) 

^'R INT 3l5f 04 (^) ,ii4 (3) »U4 (4) *U4 (5) ,U4 (6) 

315 M)HMAT 4iX,"uA *4>f tl2.bfbX»*Ub s4,El2,5»5X»*uC =* » E I 2 ,5 *5At*U0 
ltl2,5»5x,*Ut s«»ri2,5»//) 

RWINT 303 

303 >-OHMAT 'lsnA»*HW£SSllH£S*»/) 

PRINT 304 t pH (2) »u>4 (3) »P4 (4) ,P4 (b) ,P4 (6) 

304 cortMAT (ilX»*PA , E 1 2 . b » bX t *P-J ,E 1 2 , b » bx » *PC a* , £ 1 2 ,5 » bX » *PO a* 
ltl2.b»5Xf*Pt a«f Fl2tb»//) 

PRINT 305 

3p5 M)HmaT (4gX,o IENPr-KA lURtS*»/) 

Print 30btT4 (?) ,t4 O) » T4 (4) ,Th ( b) ,T4 (6) 

3pb M)RMAT (11X»<»T A a*»El2.b*bX»*lri ** * E 1 2 • b » 5X » « TC ** » E J 2 , 5 * 5X » *TD a* 
lP-l2,bt5A,*i- a«»c-12.S,//) 

PRINT 307 

307 t-ORMAT (boX,#L)ENSTTl£S<»»/) 

PRINT 30A»';jH (2) » v4 (3) fH4 (4) tR4 (b) ,R4 (6) 

3oe t-ORMAT UiX,hrhOA a* » E 1 2 . b » 5X » *PH0r3 a*,Ei2,5.5X»*RH0C s4,Ei2#5*5X* 
l*RHOU a«,Eu.5tb.'f <»KhUt s#,El2,b»//) 

PRINT 309 

309 I-ORMAT (boX,*tNTH^LHIEb*»/) 

P‘--INT 3l0»H4(?) * ^4(3) »H4<4) *HH(b) ,H4(6) 

310 KJRMAT (iIX,H^,^ a<nei2,b,bX**Mri a* , p ] ? , & , 5x , «HC ** » E 1 2 . 5 » 5X » *H() a* 
ltJ.2.b»5X»«'Jt a«»r. 12.5,//) 

PRINT 311 

311 l-ORMAT (47X,*tF^Er7IVE GAMMAS**/) 

Print 312»'^4 (2) *'=!4 (3) »G4 (4) ,04 (b) ,g4 (6) 

312 CORMAT U 1 X»*bA ai»tEl2.b»bXt*Gr) ** ♦ E 1 2 • b * 5X » *GC **» » E 1 2 . 5* 5^* *GD ** 
lLl2.b»5X,*GL sAfc I2.5f///) 

Print 313 

313 t-URMAT ‘ l3X«*REGInN 5*»//) 

PRINT 3l4 

PRINT 315*Uh (7) »h 4 (B) »U4 (9) f U4 ( 10) ,U4 ( II) 

PRINT 303 

PRINT 304 ♦ PH (7) »P4 (8) »P4 (9) ,P4 (10) ,pf, 

PRINT 305 

PRINT 306*TH(7) »t4(8) »T4(9) ♦T4(10) »T5E 
print 307 

PRINT 308*Ph (7) »m4 (8 ) i K 4 ( 9 ) t R4 ( 1 (j ) ,RH05t 
PRINT 309 

PRINT 3iO*HH (7 ) »r 4 (8) iri4 (9) *H4 (10) ,H5E 
PRINT 311 

PRINT 312»G4 (7) »<;4 (8) ,G4 (9) ,G4 ( lO) ,65E 
PRINT 401 ♦iJNUnlF 

40l t-OWMAT (UX,*UNUUtF a« , E 12 .5 ♦ /// ) 

PRINT 3lR 

316 format < IfeX^HROiNTS OF INTEREST IN THE FLOW FIELO ANU OM TFit BODY 

iiN Inches)**///) 

PRINT 317 

317 format (13X,h1nITtAL POINTS ON BOw SHOCK WAvF*) 

PRINT 3la*X0»Y0 

318 FORMAT 'llXi*A. s*,E12,5*bX»*Y «**El2.5) 

PRINT 319 

319 Format (/» i3A,Hwi\iO intersection point**) 

PRINT 318*XR0,Y3.m 
PRINT 32ft 

320 format (/♦ i3x»**aow SFiccN-wiNG Shock intersection point**) 

PRINT 3l8tXl»Yl 

print 321 



32i ► 0 ><MaT (/f 1 r'AN I Nl tWSECT lON P 0 INI 6 *) 

UO lLr 2 ,,< 

8ti KPiNr 3l8 » ( TL ) , Y4 ( IL J 

kkInT 400 

4PU r 08MAT ( AHi ♦ ibX ♦«(- ALOMG THE WIiJG IM THE EXPANSION REGION**///) 
PHlNT 4U5 

40 5 hOHMAT ' 1 3 X ♦* VEL Ur If I E S* »/ ) 

PRINT 4 02*U* (3) ♦'.* (5) ,U4 ( 7) ,U4 (9) 

4oa hOHMAT ( 1 ]X**U4 Aw s*»E U.S*3X**Uahw r* , E i2 , 5 » 3A » *U4Cw a» , E • 5 ♦ 3X * 
1 *1, 4DW s*»El2.5»//) 

PhJNT 403 

403 hORMAT 4.3X,*^^PtSSUHES*»/) 

PRINT 404 jp* (3) *r4 (b) ,H4 { 7) ,p4 (V) 

404 l-ORMaT (liX**P4Aw a*,E 12.3*3X»*P4hw a* , 1 1 2 . 5 * 3A » «P4C w a* , E 1 • 5 , 3X * 
l*P4Dw =**EIa.5) 

PRINT 500 

500 l■uRMAr(///*ibX,*^LOw CONUITIONS IN SHOCKED R£UiUN 01- EXPaHSION FaN 
1 **///) 

IF ( JSHOCn.ER. 1) Gr TO 702 
PRINT 501 

5ni ^ORMAT (<i6A,*KFGl 'IN 4ES**//) 

PRINT 502» MT S<* ( -) » XS4(?)* rS4(7) 

502 K)RMAT (13A,*ThE Ta s* » E 1 2 . 5 » 5A ♦ *A a*P T ? . 5 * 5 x » *Y s*»El2,5»/> 
print 503» 034 (o) *PS4 ( 0 ) ♦ TS4 (5) »rS4 (6) »hS4 (R) »G54 (6) 

5p3 roWMAT (13A,*.J a*,E;12,5»JX**P s*tFl2.5*3A**T a*ftl2.5*3Xi*H S**F12. 
l5.3x»*H s**c.l2.5.3x»*G s*»El2.5»//) 

PRINT 5U4 

504 roRMAT (2bA»*REGIrN 5 as**//) 

PRINT 502» ! nTS* { ■») » XS4(b), ^54(8) 

RRINT 503» OS4 ( 7) *PS4 (7) *TS4 (7) fPSA (7) »HS4 (7) »GS4 (7) 

PRINT 505 

505 TORMAT (2AA,*RF6lrM bpS**//) 

PRINT 502»lHTS4(O» XS4('^), YSAI*;?) 

PRINT 503» U 54 (8) »PS4 (d) »TS4 (B) »RS4 (8) »HS4 (H) »GS4 (8) 

PRlNl bOo 

bob l-OKMAT {26X,*kEGI->A' bCS**//) 

print 5C'2»THTSa (R) » AS4(U»), YS4(lo) 

Print 503 * 0a4(R) *PS4(R) *TS4(R) »RS4(9) tHSAfR) »US4(R) 

PRINT 50 7 

507 TORMaT (26A,*kEC-I ;N bCS**//) 

PRINT 502» lr«TS4 (1 0) » XS4(11)* YS4(H) 

PRINT 503* OS4U«) »PS4(lu) *754(10) ,RS4U0) *HS4lIu) »GS4(lo> 

IF (JSHOCk.nl. 1 )Go TU 704 
7o2 PRINT G02 

602 format (26A,*CUHVp|) SpOuK NOT 'MODELED*) 

704 continue 

NCASEall 

OALL PiGYttAN (HST^G) 

GO TO 99 
701 PRINT 601 

60l t-uRMAT (26 X,*nOT /i TYPE VI PATTERN*) 

9S continue 

LNO 

— -bUdROUTlNE REL6A«;(U1,PI *HH01 ♦HI ♦OEL T A * NUPT ♦ J2 * P2 ♦ RhO? * r2 ♦ T2 * 22 * — 
|OAMMa2*Thlta,s2) 

COMMON/ ShOCn/NFL' I d* I SHOCK *JSmOCK 

OELTAb0ElTa/57.2q6 

IHETALaUELTR 

IhETARsI.571 

5 ( HETAa ( T hETaL ♦ tH£TAR)/2.0 



O O o o 


IF (Theta, LT.1.57-79)GU to 7 
If (NFLUti.OT. i)<30 TO 8 
IsMOCKsl 
HfTUHN 
8 JSHOCi^sl 
HF TUHivi 

7 Continue 

I anT=tan ( Tmc-Ta) 

I AiMTHUsTAlM META-OELTA) 

8J nTsSJN ( Thr. r A) 
bl nTmUsSIN( IhETA. DELTA) 

COSTsCUS ( I Md A) 

EhSsT anTno/ I ANT 

HHU2bsRH0l/c-PS 
01NsUl*Sl''fT 
U1 f sui*'-ObT 
»-'2NsEPb*UlN 

o?r*u t 

P2=P l+PH0i*ulN<njl M-WF'02 (j*U 2N*02N 
H2=Hi •5*U2 N*u2N 

u?=U2N/SllMTi'ii.) 

h2 = h 2/ (32. 1 /b«77fJ.O) 

Call .'iULItp lti?.P?»NUP r .T2»Z2*S2»HHC?.GAMMA?) 

HijiFF sRH02b-HM0? 

LONVsAbS (HDiF F/KH02(i) 

IK ( COwV"(' • 0 1 ) 2 , 2 . b 
t J,F(ROIFF) 3»2»A 
3 InETARsTHETA 
ot) TO ^ 

A IhEToLsIhLIm 
«0 TO 

c LONTiNUE 

M2aH2*32, 170*778,-' 

UELTAaDELTA-»57.2 

IHETa=THFVa^57.2oH 

Hh f UHN 

LNU 

— bijHfiUUTlNE iNTNST(Xl»Yl»Al»A2»X2»Y2»XI»ri) • 

v****SUoPOU I iNt rO FIND THE INTERSECTION POINT OF THE SHOCKWAVES 

ANu THE Intersection points in the f:xpansiom fan***** 

Ais ( Y2-y l-»Xi*TAN f Al) -A2*TAN (A2) ) / (TAN (Ai ) -tan (A2) ) 

YlsYl-f (Xl-xi)*TAKi (Al) 

return 

tlNjU 

•bUHROUTlNE HlGYd'K (HS) » 

uJMENSIOn F(20U) ,0(200) »FN(200) ♦bN(?f*0) »YN(200) »ETaO( 200) ♦ETA(200) 
lfTHTTA(2oO) 

UIHENSION Y (b) »FpA (3) .ljA(3) 

UOUHLE PpEClblON YtAAfrib.CC.UO.EE 

ut MENS I On Dy(5)»P(5)»TE(b).YPR(b.4) <5) 

OIMEnSIOn ER (b) ti'PSAvE (5) 

dimension XA(50), Y4(50)» P4(bO)» T4(5U)* U4(bu) 
double precision DPSaVE»OEP 
external OERIV 

CoMM0N/C0NST/a1» a2»A3»A4»E»EPS 



COMMON/KONSI/Twai L*HauIUS»PH 

l.OMMUN/WING/PT3»TT3»U3tTJ*P3»X30*Y30»X<t»y4*PAtFt>»T4»U4«T‘St 
L0MMUN/CASE/|\iC4SF 
if (NCASt.GT. 1) bn TU 1100 
KE:«0 It tt utLT* EPS 
1 »■ (iWMmT ( 3t. ii? • 3 ) 

1100 CnNTlNLlfc. 
m)=u3*u3/ 

A3 = 1.(:/PP 

m 4=U3*U3« (PH-l ,0 W (hS*PH) 

IF (A3) 

?U COi'fTli'^Ut 
Art = 0 • 

dH = 0 . 0 
CC=0.A7 
ui.) = TwaLL/TT J 
tfs ( 1-Ob) *0.n 
IF (NCASE.l’f • i) Tu liOO 
HEAD l3»KtY 
]3 t-ORMAT (lib) 

HEAD lA»INfALF 
lA f OHMAT ( 1 l<JtEi2,b^ 

1200 LONTlNUt 
Ml) sAA 
» (2) =HB 
» (3) sCC 
t (4 ) sQD 
T (b)=EE 

3 bMSAVE = Y(3) 
r PPSAV = Y ( 3) 
r PPUi = Y (3) 
oPOl = Y(b) 

EPPOb = Y(3) 
bPOS = Y( ) 

II = j 

30 call NUMJim (!N.tnELTt T tQVtPt rE»EH* YPPtYtEPStDP$AVE»KEY»iiERlO* 1) 

1 = 0 , 

I'isb 


1 = 1 

C.\LL NUMliM (NtDELTt T » C V t P t TE ♦ £R t YPR t Y t fRS » opS a V t » KE Y • UER I V ♦ 
call NUi'iI ) C'<»f)£LTt T tCV tPf TEtEH, YPr, y tFPStDPSAVEt''.EYf uERlO* 
IF (ABS (Y (3) ) ,L I .E) fjO TU 33 
1 = I ♦ I 


uu To 5 

33 IF (aBS(1»0oo 0-Y(2) ) ,bl ,£) GO TO 34 
IF (ABSd .Oooo-Y (4) ) .LI.E) go TO 9<5 

34 ERA (I I) = ,M2) 
oA(Il) = Y(4) 

II = II + 1 

IF (1I.EQ.2) GC '0 35 
IF (II.NE.3) GO TO 36 
T ( 1 ) = Aa 
T (2) = Db 
Y(3) = FPPsHV 
Y (4) s CD 

Y(5) = 1.0b*bPbA,/E 
EPP03 s Y (3) 
up03 s Y(5) 

EPPOS = Y(3) 
oPOS s Y(5) 
bO TO 30 


U 



35 n 1 ) = AA 
» (2) = db 

VO) = l,0S<»rpPS'.v 
» (4) = UO 
rib) s OPbAvc. 

rpp02 s y(i) 

Pp02 s Y(b) 

(-PHOS = Y(3) 

OPUS = Y(") 

^0 TO 30 

Jt UFPA s FPA ( 1 ) -FPa (2) 

UtiA s GA(1)-«A(2> 

UFPAl = Ff^A ( U -Pt-A (i) 

UGAi = bA ( i » -GA (3) 

UKPPO = PHPU1-Fpp<.i2 
UbPO = GpOi-bP03 
= OFPA/OPPPr 
cs| s uFPa i/UoPU 
Oi = 1 ,0000-PPA ( 1 ) 

Ul s UGA/UFPPt 
tl = uGAl/Uop^i 
PI s i.UOU0-GA(l^ 

UfcFPP s (Cl^t l-F 1*01 ) / ( Afe*El-01*Bl ) 

UtGP s (Ab*P l-Dlu-Cl)/ (A6*EI-01*H1) 
i-pp a FPPOi ♦OfP P'-’ 

OP s tiPCl+OcoP 
T (1) s AA 
Y (2) = PJb 

' O) s FpP 
T (4) = UO 

T(5) = bp 

OO TO 3 
9*^ CONTiNUt 
T (1 )sAA 
r (2) sHO 
M3) s FPPOi) 

Y (4) = Up 
T(5) = bpOS 

call NuMItMNtntLTOtOVtPtTetErtfYPRtYfPPbtDPSAVtf KEYtDEPlV*!) 

1 * 0 . 

W a b 

A = 1 

ULSTI a ( Y (- 4 ) -A1 » ( Y (2) **2) ) /A2 - Y(2) 

IhTAI a ( Y (^) )«(i .0-Y ( 2 ) ) 

IhFEI a (Y (2) )*M .0-(Y (2)**2) ) 

111 * 1 
P ( 111) a Y (2) 

0(111) a Y(*) 
tTAO(Ill) a r 

Call NUMitM(iM,0ELT*TfCV*P*TEtE«*YPR,Y*EPi»»0pSAV£*i<EY*DEHlV*2) 
ULST2 a (Y(‘*)-A1*(Y(2)**2) )/A2 - Y(?) 

IHTA2 a (Y (2) )*(i ,0-Y (2) ) 

IHTE2 s ( Y (2) ) « (1 ,0-( Y (2) **2) ) 

AOELST a (o.b«CL«:Tl ♦ ULbr2)*0ELT 
mTMETA a (o.b«THTAl ♦ THIA2)*0ELT 
ATMETE a (o.b*THTEl ♦ THTE2)*DELT 
55 Call NUMIN (N,0£lt* r*DV»P»TE*tH»YPR»Y»EP5*DpSAVfc*GEY»UERlV*3> 
UOELST - ((Y(4)-al<»(Y(2)**2))/A2 - Y(a))*D£l.T 
UTHETa a ( ( Y(?) )*(1,0-Y(2) ) )*UELT 

OTHETE a ( ( Y (2) )*(1.0-(Y(2)**2) ) )*CFLT 

aOELSI a AOtLST ♦ UUELST 
ATHElA a ATritTA ♦ DTHLTA 



'^THtTt = ATrltTt + UTHtiE 
ill a 111 ♦ 1 

tiADdll) = 7 

1“ ( 1 1 1 ) a r ( <i ) 

0(111) = ^ (‘♦) 

IF (AtJS < y i J) ) .LT .E) GO TO 990 
1 a 7 ♦ 1 

«0 TO 

990 Continue 

N 1 S i i 1 

MA] s A1 

'.njM a 0,c. 

C CROSS SIhEaI'I station IRaNSFORM with -in- no. of stations 

UO 751 i = 1»N1 
YNd) = 

lF(ALF)92»9^f91 
92 iF (F ( I ) -.99) 91 »9'' .90 

90 AL.Fs3.o/fc. 1 AiM I ) 

91 CONTIimUE 

a oui^i ♦ ,i)5 
/5i Continue 
0 s 1 

ou 752 1 a i,iiv 

CTAd) = ALUOd, Vd.O-YN(I) ) )/mLF 
00 753 JJ a 1,M 
IF (El A (l/-ElAn(J) ) 76 0.761,762 

762 Continue 

J s J+i 

00 TO 75"^ 

76 1 continue 

( 1 ) = b ( J) 

CNd) a F(J) 
oo To 

7feO J a j-7 

IF (ElA(l) -El mHCJ) ) 770*761,772 
770 KwlNI 7B0 
7H0 TOWMAT (2A* funny *) 
oo TO 7^9 
772 continue 

AF s (tlAd)-ETA!i(J) )/(ETAU(J*1)-ETAD(J) ) 

FNd) s aF*MF ( J+1 ) -F ( J) ) ♦ F(J) 

UNd) 5 aF » 1(5 ( J + d -b (0) ) ♦ G(J) 

7b3 Continue 
752 CONTlNUt 

00 7b i 1 a 1 , IN 

(HTTA(I) a uiM(l)»AAlwKWd)*FN(I) 

7B1 continue 

F N (I N ♦ 1) a 1 , () 

W(dlN ♦ 1) a 1,0 
IHTTA (IN ♦ A ) s 1 ,0 - AAl 
799 continue 

Call t JOYCE (FN,T hTTA,ALF) 

999 CONTINUE 

return 

EnO 

MM* SUhROUTInE OEHIv (T,Y,UY*L0C»N) ■!! 

0IMEl^iSI0^ Y (5) »Oy (5) 

UOUBLE precision X,X1,62 
Common/Cons ( /A i»a2, A3* A4,£ 

AS ( ( Y (4) - (Al*Y (2) **2) ) /A2) 

Xla(A** (-0,45) ) 

X2a(-o,25)*(X*« (-1 .25) )* ( (Y (5)/A2)-2,*(A1 /a2)*Y (2)«Y (3) ) 



on I ^ 99 


7?0 CONllMUt. 

• — DOUHLt H^<£CTSIu^ AUAMS MUULTQN PfCtC 

U() 

730 uv(J) = U£P(J) ♦ I)ELT<» ( H ( 4 ) ( J) ♦!£ ( J) ) 

'-ALL UEHT V ( I ♦nv ^ ) 

/31 UiNTii'^Ut 

uo 74 c J = 1 ♦ "i 
tPP(J»M'+) = H(J) 

weP(J; = Ue.r'(J) *■ Ut-L (iM4) ( J) ♦!£ ( J) ) 

lm< ( J ) = u L H I J ) 

SI''-PL<^ S 1 1 £HHUh t:.SIIM-47t 
/•4U = K A T i Arte ( tH { J ) -iJ V ( J ) ) 

/A i t-ONTi'VUt 

oi) lU 5U0 0 

C-i-jl) 

— ailHHUUTlNt C-JOYCtr ( t- N ♦ IhT I 
I i 1 M rt 1 Ui'j f" ( S I! ) . T ^ r T A ( b 1) ) 

'o i ME'ajS lUf'; HF ( b-' ) » rtC V I N ( bo ) f I life. 7 A w ( SO ) » PAH ( SO ) ♦ i-'HG ( SO ) ♦ fe I'lAtH ( SO ) * 

) If (SO) , V 1 SE ( ■30) (30) ♦ 7HETAE (SO) ,UF ( SO) ,t)rtth«Ut (SO) ,«ExiiOl (So) , 

,>n(S0) tUSOO) fRtI A (So) ♦a(S0) ,FL.CW (SD) iFINJCSO) »t.A(5l)) 
uIMEwSIOIm fe (SO) f rHE7A (bO) fCSCSO) »Fl (S0> ♦7 hETa 1 (bO) »CS1 (b(j) »fe2 (SC) « 
lUFLCUs (St; ) ,uM CS’ (bO) » 1 (SO) ,FF (bO) ,FL.C (bO) ,Fl-Cl OO) » tHFTA^ OO) , 
?LS2 (5(1) »F LC^: (bo ) . YN (SO) 

uIMENSIOn CPI (bOi fCPS (SO) »\/IS t (bO) »VTSb (bO) » ICuiol (S') t7CoNb (SO) « 
]UTFKC (S(^) ♦CP-l(bU) (bU) ♦G?l (SO) *XS(5M) *XJ (So) fWlSMISn) ♦ ICUNM(SO) 

PfrtHOM (bw ) tC-AF (SO^ t SCfU (Sn ) ♦PRi-iO (So ) 
uImEinSIun E I a (b(M ,GE7 a (bo ) , Y (bO) »DEl. ST (bO) , THMUM (bO) 
ulMENSIUiv DbTAH ('^^^•0) »C0b7H (SO) 

f i ('IE '<b 1 Gf'-' A ( .31' » 3 . 3 ) » F ( bO » 3 ♦ 3 ) « C (bo » 3 » 3 J » AH ( S .i ♦ 3 » 3 ) » H i 50 ♦ 3 » 3 ) » MAH ( S 
1 il 1 3 ♦ 3 ) » M h I i V ( S U . 3 * 3 ) » A IJ J (i A H ( b 0 * 3 , 3 ) « A G ( b 0 , 3 » i ) * ia ( 5 0 » 3 » 1 ) ♦ ') A G ( S 0 ♦ 3 
2»] ) ♦iJ(Sti»3<i) »Fi*^ (5U»3»i) »w(SU»3»l) 
u iMENblUi. 7 aw (Sr ^ .'3UC I 1 (SO) tSTNOT (5o) ,bHECHT ( bu ) ♦ A T C ( SO ) ,WUU7M (So) 
uiMEivSiUf'i (b ) f (bO) 

uIMLhSIUim HCoFF (-i(j) »H()y/lNrt (SO) « tw (So) 

olMtiNblOKi 7t- IP ( 1'=^) » ^ IbC ( lb) , 7C ( lb) »CP ( ib) ♦ vl ( lb) ♦ V2 (IS) » Y'i ( ib) , 

1 ICX ( ib) » tCe US) » fP3 ( lb) ,LP1 ( irt) ,Cp2 (IS) »bp3 (lb) 

u I HENS I UN Xi(S;)).Yl(So)'»A4(S0)»Y4(50)»Si(b0) tiJc-Lx (5U ) »p4 ( bO) » 
ll4(bl), u‘»(bO) 

COMMUM/AOINS I /TwAi L»HAL)iUb»PR 

U)MMOlM/w 1NG/H73 » t 73,U3» 73 »P3» X3U tYSO »X4 »Y4 ♦ PA »Hb » 74*04 » 7bE 
COMMUN/CA -SE/ i iCASr 
IF ((MCAbE,G7 . I ) Gr. 70 ilOU 
kFAD 15* Pl^sYHK 

15 F0F<HA1 (E^£:.S) 

Hfe. AO 1 »MH»M»NK*NnPF(('J 7 

1 F OHM A 7 (4lW) 

HEAU c»WMif'*OiS 

2 roHMAj (2tl2.b) 

HEAD 10»TEMP*VlSrt7C»CK 
10 F0HMA7 (bFlS»b) 
iloO 00N7 (nUE 

call SPLn 7PH ( m , tEHP» V ISO* Y1 » Y2»V3) 

Call SPL'\i 7 HP ( ) 4 * r Emp « fC»1Cl»l C2 ♦ 7 C3 ) 
call SPLM I HP ( 14 ♦ tEMP*Cp»CP 1 »CP2«CP3 ) 

iSENTHOPiC FLOP oF A PC.REEC7 OaS JC COPPUFE EL)Ot PpOPERTtES 
FTESP73 
I 7£=I73 
on n I = J »5 

Al (I)=X3 '-X4C(*1) 

Y 1 ( I ) SY30-YA ( !♦! ) 



, ( I) =S‘Jh 1 (Ai (I)*Xl(t)^Yi(I)iM'l{I))/l?,U 

]1 CoixjTXKUt 

A 1 i\i T s . i o S 1 U ) 

A ( 1 ) sXINT 
>JFL=U .?^S1 ( 1 ) 

( 1 ) =Pi 
( 11 ) =P6 
i ^ ( 1 ) =T3 
14 ( 1 1 ) sT‘:5t 
U4 ( 1 ) sUJ 
1=1 
L = l 
U() i 

IF (PlGYbn ,C 1 •0 1 tjo rO o3 
Mfc AU i?fpt {C^) , TF. fM) *Lt (M) »TW (M) »KaC ('‘I) 
li! f-oHMAT ^5t:i<eo) 

*.’0 TO 1 J 

63 1F(KK,£0. J) on I i 64 

Al) J , H A^ I !S ) 

3 f OWH A T ( Gt 1 P • 3 ) 
bO To IJ 

6A kam (i'l) sRAUlub 
13 t.il'jTlMUt 

hh (M) spA (U 
It (M) =14 (L) 

UF (M) =Ua (L) 

O'-LL M0Li£H (H,Ht iM) , Tt (M) , Z # ENT , H HO ♦ 0 ) 
t'hOt. { M ) =lvH0 

V f SE IM) =SPL’'t'' AU < ' 4 ♦ IF MH» Vise » 1 » V?, V3 , TL (M) ) 

I ^ { M ) s ( iV ') L L 

lnET“r. (h) = TE(M5/I Tfc 

I hET AW (i'll =-Tw U'')/tTE 

owEMUFIi'^) = HHUF (M) itoF (M) /VISE (M) 

iELX(M)=I'El 

/. (M + I, ) =A (lA) +UELX rM) 

IF (X (M+ 1 ) .(iu.Sl ( T ) ) GO TO 61 
bO TO 14 

6i It (L,EQ»Vl bu TO bi 
! = !♦ I 
L = L + ? 

uEL = 0.2* (i>i U ) -5i ( I-l ) ) 

X (M+ 1 ) =bi ( i-i ) ♦ o,5*OEL 

bO TO 14 
b'^ u=il 
1 = 11 

i| ( I ) sSl (3) ♦X (t + 1 ) 

OFLsDtLX ( 1 ) 

J4 bONTlNUt 
KRINf 9 

9 format ( ] Ml ♦2X»*nUTPU T UATA«///) 
uELT = 1.0/(i^i-l) 

ivfM = N-1 

c'iMAX = 

KEXIiST ( 1) =UMtr.iCE ( 1 )*.X1M 

U(j 1*- M = 2iHr'i 

RF X IM ( M) = Rc.A IK T ( M-1 ) + ( OREnOE ( m- 1 ) ♦iJhENOF ( m) ) auELX ( M) / 2 • 

16 COMTii''UE 
PRINT *7 

1 7 F OHMa T ( 2 a , ♦ 6 y ♦ *PE ( M ) » , 4 X » «THE T A F ( M 1 * » 6 X ♦ *UE ( M ) » ♦ 

] 4x »«URENOt. (M) #»3y »<*HExlN 1 (M) */) 

UO 19 MS1»MM 

PRINT 1 H ♦M,F'E (P) . THET At (M) ,UE (M) ♦URFMOt (M) ♦REXiN f (p) 



o o n n 


Ifi I- OHM AT ( 1 l5»bEi<?.5) 

1 'i LOI'JT li'iUt 

tVALUAT lOl^i uh S 

5(1) = HhUE U ) AV t SE 1 1 > *Ut ( 1 ) * ( (HAH < 1 ) ) *A InT/ (:3•0**^K) 

L'U 3l' M = ?tMi'iAX 

5(M) = S(M-i) •*• (khUE (M)*VISE(M) Al)E (M)<>( (HAn(M)*<>KK)<»*?.) 

1 ♦ HhOF (M+i ) I SE (^♦i ) *UE (M*l ) * ( (RAH (t'Hl ) ) ) *DELA (tA) /2,^ 

30 C(jNT l^Ub 

5 (MM) =s (MMAA J (HnOE (MMAX) (MMAX) «ut (MMAA) »( (HAi)(MMA*) **K) 

1 «A^) ♦ HhOE (hM) *^/ISE (Mi"i) *iJE (MM) « ( (RAO (MM) **KA )**<?) ) *' . 0 

i-s ( 1 ) = S (if ) -S ( 1 ) 

hETA(I) = (A,o«(9(2) ♦ S(l))/(UE(2) ♦ UE(1) ) )*( (UE(2)-iJE(i) )/(S(2) 

) ) 

U() 4J Ms <i».'1MAX 

US(M) s (b(i'^l)-<; (M-1) )/^,0 

beta (M) = (<r,(i«S (M) * (ut (M+1 ) -UE (M-1 ) ) ) / (UE (M) * (S (M+ I ) -S (M-1 ) ) ) 

41 LOI'JTInOE 

HkIiNT 43 

4 3 I- OHM A I {/EA »*X (M) «»8X»*S (M) 7X »<»US (M) *»bAt*d£TA (M) <*»t)X» 

1 ’♦ELO'm (M) A/ ) 

UO 4b M=l»Mi-)AX 

HHINT 44»i'»A(H) »-;:{M) ,JS(r1) .HETATM) «FLCw(M) 

44 t-uHMAT ( I lb»bR12.S) 

45 COivlT 1 Mlt 

kklNT “tBl ♦P I L » T Tf »ALF 

451 i- DKM-T (///ioy I E=« »t. 12,5 ♦! OX* 4TTF=4»t.l2,S flux ♦*aLFs«»el<J»5//) 
f, s 1 

UO 33 I = 1*N 

rod) = ='»0 ♦ (Ui)AijtLT 

IF (PIGYH^4t0.ij. ) OG lU 35 

M 1 ) s F U ) 

hitT/'(I) = InTfAd) 

CSd) s 1*0 
GO To 3 fc> 

35 UONTINUE 

kEAU e2,F(i) .THt TA(1) ,CS(X) 

32 FOKMAT (3tU.5) 

36 continue 

M ( 1 ) s F ( I ) 

I hE7 1 d ) = j HE r •' ( 1 ) 

u s 1 (I I = CSd) 

33 continue 

F 2 d ) s F I ( 1) 

F 2 ( 2 ) = F i ( 4 ) 

00 34 I s l»i«i 
UFLCUSd) = u.O 
OFLCSld) = DFLrUSd) 
id) = TT£»friFTAd) 

34 continue 

UO 5 ' 1 s 1 ♦ NN 

FFd) s F ( I) / ( ALf» d , 0- ( 1-1 ) *OELT) ) 

50 CONTINUE 

EVALUATION OF Tmp stream FUNC T lOfNi ♦ F L C 1 
FlC(1)=FLCwIm) 

FLC(2) s F'LC(l) * (F(1)/ALF ♦ F <2 ) / ( ALE* ( 1 . 0-UtLT ) ) ) 4DELT/2. 0 

UO 51 I s 3»NN»2 

FLCd) » FLkd-2> ♦ {FFd-2) ♦ 4,G*FFd-l) ♦ FF ( I ) ) *DELT/3 , o 

IF ( ( J *1 ) .G ( ,NN) GO TO 51 



>-LC(Ul) s fLC(l-l) ♦ (KF(l-l) ♦ 4.0 *fF(I) ♦ FF ( uu ) *n£Lr/J.O 
6 1 CONT INUt 

uo 1 s ItNN 
FLCI (I) s FLC( I) 

LONTlNOt 

FRINI 

5j f ORMAl ( IHI »^Xf * T*»bX »<»F ( T ) »»6Xt*THFTA ( i ) *«<^A t*Cb ( I ) «/) 

L'O 5b 1*1 ♦M 

rRiNI 54,X»F (I) tlHtTA (I) »CS(I) 

5A FORMAT (lib OF 12, 5) 

5b CONTl'MUt 
1 * 1 

tPSU) sSPLMML ( 14* TtMPtCF,CPl »CP2,CP3» f ( 1 ) ) 

I CONb ( T ) rbKL.W4L f 1 4 * T tMP , T C * T C I » TC2 * T (^3 » I ( I ) ) 

V IbS( 1) =SPLi'i'/AL ( 1 4* TF.fiP» Vise* VI » V?, V3* n 1 ) ) 
upl ( i ) =CPb { i ) 

I CON 1 l T ) * • '(S ( I 1 

V isi ( I ) =v rss ( 1 ) 

CALL MOL] tp (H,PF f M) ♦ I ( 1 ) *Z*FNT »i3pn*0) 
khOM ( '! ) SKnu 

CkMU) s CPld) * Cb(l)*(CPS( D-CPI (1) ) 

IF' (Cb(l)-l»t/) b=(0*b5i*5bl 

550 AS(1) = (CSU) * '*MI)/(WMS ♦ CS(l)«(WMl ” WmS)) 

A r ( 1 ) s 1,0 - A b f 1 ) 

IF ( A 1 ( 1 ) ,tw. 0. (d Ou 10 b51 

012(1) * ((i.O ♦ ( ( V Ibb ( 1 ) /ViSI ( 1 ) ) <H»n-b) * ( ( wMl/viiMb) **0*?P) ) «*2, ) / 
1 (2.82H4<» ( ( 1 • w ♦ >'b/ivMl) <»*n,5) ) 
o21(l) = ((i.O ♦ ( (VI bi (1) /VI bS(l) )**o,b)#( (F.‘Mb/i-vMI)«*i) ,2b) 

1 (2.8284* (( i ,u ♦ ^MI/v»MS) **0,5) ) 

viSMd^, = vibs(n/(i.') ♦ 6i?(i)*Ai(n/Ab(i) ) * visi(n/(i*o ♦ c 2 i( 

I i ) *Xb ( 1 ) /A I ( I ) ) 

ICONM(l) = ICGFb(l)/(1.0 ♦ l,v)8b*Gl?(l)*Al (1 )/AS(l) ) ♦ TCUiVKl)/ 

id.O ♦ 1.06b*C21 d)*Abd ) /AI (1) ) 
uO To 552 

551 vjSMd) = ViSS(d 

icONf'di) = icnixSd) 

55c CUNTlNUt 

I aU s RiIoM ( 1 ) »V ICM d ) *Ut (M) *UF. (M) (PAP (M) *«KN) *ALF* (F (2) -F ( 1 ) ) / 

1 (mELT* ( (2,0*b (M) 1 **0,5) ) 

CF * TAU/ (0*b*HMoE (M) *UE (M) *UF_ (M) ) 

npAT * - TCuNM(l)*TTE*OE(M)*pH0Hd)*(RAU(M)**NK)*ALF*(THF:lA(2)- 

T I HtT A ( 1 ) ) / (utt T* f (2. 0*b (M) ) **'j,b) ) 
aiDUT = “hLa I 
J a N 

CPS(I)=SpLNVALd v*TtFP»CP»CPl*CP2,Cp3,T (I) ) 
lC0Ns(I)sSPLNVAL(14,TtMp,TC»TCl*TC2»TC3*l (I) ) 

V ISS ( I ) aSPLNVAL ( 14* 1£MP*VISC» V1*V2,V3*T (I) ) 

ICONM(i) = IcnObfl) 

V ISM ( , ) = VISS ( 1 ) 

CPM d ) a CPb ( I ) 

PRNO(i) K VibM ( I ) *CPM (1 ) /TCONM ( i ) 
hFC = SuhT (PHnC ( ') ) 

IREC = 8E'-*( TTE-TC (M) ) ♦ TE(M) 

alNO a Oi'iOT/ (Phor (M) «Ut (M) *CPM (N) « (TREC-THETaW (M) *TTF) ) 

PRINI 556»T4u,CF 

556 r OHMAT (//’^A **1411-** E 12.5* 1 OX **CF a»,F12*b) 

PRINT 557 * I pEC, QDUT *bTNO 

557 FORMAT ( //V A » *THEca* » E 12 .5 ♦ S X ♦ *QUOTa* *E 12 .5 * 9 A » *b iNOao 1 2 • 5// ) 

Efc CoNTli'^Pt. 

59 continue 

UO 5/ 1 s X»N 

c evaluation uf TMr thekmo properties of the components 



LP5( i) sSPLNVAL( TEf^H,CH,CPl»CH2,CP3, T U) ) 

I CUNb ( T ) =bPL,>IVAL I 1 A»TtMP» rc» TCl » TC2 * 1C3» I U ) ) 

V ISS ( I ) =SPL ( I ^♦TEi'iPt V ISC# A » V;>» V3»T { i ) ) 

SPI ( 1 ) sCPb ( i ) 

I CON I { I ) S I ::UNS ( 1 1 

V isi ( I ) sviss n ) 

C-----tVALuA I Ion Uh TP'- IHEhi'IU PWOPEHTIES of fHF MlA(Uwt 

hi =(wMi//::'b.u)<‘*.Abl 
hS =(w M S/2 f, •!')■***. 4bl 

I'lhfCdJ = <;.b‘*n/(10.**B)*(T(I)«*l,S)/(PE(M)*Pi*hS) 

F7 Ci)iVlT I l\iU£ 

U() SB I = l»Nh' 

IF (CS(I)-l.u) /P»71#hi 

71 cijNTiNUt 

CPM(l) = Cpb(T) 

V I SM ( ) = ViaS ( 1 1 

I rONM { I ) = i 0(■.^ S f I ) 
wC TO f3 
12 oontINUE 

oPM(i) = CPi(T)+ CS( [)<»<C^S(I)-CPT (I) ) 

.712(1) - (U.O ♦ ( (VlSS(i)/VlSI (I) )**n.b)«( (*#Mi/wMS)o»0. 2b) )**?,)/ 

.U2.ti2B4<*( d»o ♦ .MS/i^i-H) #n»y ,5) ) 

P21(i) = (d.C ♦ ( (V iSI ( i ) /VISS ( I )) «( (^'"^b/wMI ) <»«0. 2b) ) *#2* ) / 

H2#H284* ( ( 1 • V ♦ -I'll/wOS) ««0,b) ) 

.^S(1) (CS(l)<>^-< d)/(wKiS ♦ Cb n ) * (WMI-.#MS) ) 

Af(l) = l.:^-AS(l) 

IF (>I (1 ) .fc>.0.0) Ou fu a 

vtSMlD = ViS-(l)/(l,^ ♦ 612 U) <»XI ( T ) /XS ( I ) ) ♦ 

1 V I SI ( 1 ) / (1 . J ■Mi2i ( I) « ( AS ( I ) /XI ( 1 ) ) ) 
iroMvi(i) = ( CONS ( T ) / ( I •() * 1 . i)6b«Gl? ( I ) * (XI ( 1 ) /XS ( I ) ) ) 

} ♦ lCOMfI)/(1.0 ♦ 1 .Obb«s?l (T>*(XS( D/XI (I) ) ) 

73 OOwTlNUt 

oaLL MOLltWlrifFt fiv!) »2# I ( 1) »Z#1 .nI »PH0*6) 

PHOM(i)=HHO 
Fti CONTlNUt 

c- — ---COMHUTmIIOi'* of kLOw PNOPEPTIES at the d. L • t^OOE 

i = 

ICONM(I) = iCONSfl) 

V I SM ( ' ) = ViSS ( I ) 

CPM(I) = CPb(T) 

C ALL I'iOL I LR ( n , Ft f M ) ♦ 2 » T ( 1 ) ♦ Z » t.N T » RHO ♦ 0 ) 
khOM ( I ) =RHO 

RUMUt = RHOu (-vl) <*■/ iSt ( W) 

00 6u T = i#i'l 

CkF(I) = HHUM ( 1 ) I SM ( 1 ) / ROMUE 
SCNO(I) s V ASM ( 1 ) / (HHOM ( I ) *DIFFC ( I ) ) 

HRNO(I) = V IShM n ^’•CKM U ) ZTCONM ( I ) 

60 continue 

C computation of I-E MAIhIa elements 

Of) 9-7 I = 2 # Nit- 

«LE r = ALP« ( 1 . j- ( Ul ) <^CELT) 

M ( I# 1 #1 ) =- (CkF ( I ) a (AL h r««2. ) / (2.0* (nFLT**2 , ) ) ♦ ( alF T/ (4,0*UELT) ) 
1 ^(ALF<#cHF<i)-FL(,:( I ) -ALh T<» (CRF ( 1 + 1 ) -CRF ( I-l ) ) / (2 . > *OELT) - 2 * 0 *S (M) « 
pUFLCUS ( I ) ) ) 

M(l#l»2)=0»u 

A(I#i#3)=u«i^ 

M ( 1»2# 1 ) = WRF ( 1 ) * (UE (M) **2.) (ALFT**H. ) * ( F ( I ♦ 1 ) -F ( I- 1 ) ) / (CPM( I) *T 
1 l€<»4,0*(oELi<»*2, )*777.y) 

A ( 1 #2#2) =- (CkF ( I ) * (ALF T*»2. ) / (2 . ■'* ( DELT **2 . ) *PhNO( I ) ) ♦ 

1 (ft LET / (4,a*UtLT) )*(AlF*ChF ( 1 ) /PrNO ( I ) - ( <^CFT/CPM ( I ) )* 

2 ( (CRF (1*3 )*cPM ( I.l ) /PRN0( !♦!) -CHF ( I-l )*CPM ( 1-1) /PHM.' ( I-l ) ) / (2,0 
S«OELT ) ) -F' - ( 1 ) -2 . 0*S ( M) *UF'LCOS ( 1 ) -CRF ( 1 ) * ( CPS ( I ) -Cp I { I ) ) * AlF T* ( CS ( 



<n4.n-csn-i) )/(Ci.M(i)4»sCNO(i)«2.o»cfri t) ) ) 

A ( I »3) s (ALH ( 1 )* (CPS ( I ) -CPI (i) ) » ( THlTa ( I^l) -THt lA ( 1-1 ) ) 

1/ (CPii I 1 ) *SC''0 ( I ) <> (H.O) * (uELT**2. ) ) 
m(1,.5,1)sO,u 
A ( I ♦ Jt2) sO.v/ 

A(l,J»3)s-( IAI.F I ) *CRF (I)/(2.0*(nFLT**2.)*SCNu(I) ) 

!♦ ( ALF TV (a.qAuEL ^ 5 ) » ( alF*CPF ( I ) /SCN'C ( I)-FLC { T ) "2 • u*S ( m ) «OF UCUS ( I ) 
?-rtLF T*(CkF(l + l ) /<--CimO ( 1*1 ) -CRF ( I-l ) /SCNt) ( I-D )/(2.0«UtLT) ) ) 

» (I ♦ 1 1 1 ) = ( Ji-KT*« ’. ) «CkF (I ) / (DtLT<»<»2 , ) ♦bt T A ( m) *F ( l) *2 , 0*S (M ) *F ( I ) /O 

1 6 ( M) 

( I » i » 2 ) = U . 

( I 1 1 ♦ 3 ) = 'V " 

o( I,2»l) s riFTA( t)*(bF(M)**?.)*i-<H0E(M)/ (CPm(U<^I l E*HH()m ( I ) « T 77 ,9) 
c5(I*2»2) = (ALFTA*r,)frCWF ( 1) / ( HP(nU ( T ) * ( pt L I *#? . ) ) +2,0AS(M)*F (l)/r)S(M 
1 ) 

o ( I ♦ 2 ♦ 3 ) s 0 . 'J 
» ( I , J , 1 ) = v) , 

( I » 3 ♦ 2 ) = 0 , U 

o(I»3t3)=(ALFTAA>,)«c«F(I)/(SCNU(I)#{pELT*»?.))^2»'*S(M)*F(l)/OS(M 
! ) 

t ( I ♦ It 1) =“A ( 1 , 1 * 1 ) - ( ALF I*A2. ) AC.-(F { I ) / (ULLTtt'Op.) 
b ( 1 ♦ 1 »2) =0, J 
C(l»l»3)=*-.»u 
C(I*2tl)=-A( 1,2*1 ) 

U ( I =-A u ,2» p) - ( ALF T«»2. ) <*CPF ( I ) / (PkNO( I ) * (DEL ) ) 

b(l,2»3)=“AlI,2»0 

C( I»3*l)sU,u 

b(I*i,2)=b,U 

C ( 1 , 3,3) =-A ( 1 ,3* 0 - (ALF T *»*2. )«C><F ( I ) / (SCNO ( I ) * ( UtLT*<»2 . ) ) 
u ( 1 , 1 , 1 ) = ( Ai-F T«A;>. ) «C»^F' ( i) A (F ( 1 + 1 ) - ?,*F ( I ) + F ( I-i) ) / 

1 (2.0* (ntLT*<>2.) ) + (ALF r« (F ( I + D-F ( I- ]) ) / ( A . 0 «u£L T ) )*(FLC ( l) ♦ 
?mlFT*(CHF ( !♦! ) -C 'F ( I-l ) ) / (2.0<*OtLT) -ALF*CWF (1> +2 , ( M) *|)F LCDS (I ) ) 
3 + hETa (M) »PHUc. (i“ ) /Phuy U ) *2. 0*3 (ri) « (F ( I ) **2, ) /Ua (M) 
i'(I» 2 , 1 ) = (ALFT**7.)*CPF (n*(TriEU(I^l)-2.0*TMEU(I)*THETAU-l) )/ 

I (2.0*(DELr**2.)*pDNUa) ) ♦(ALFT'»nHETA(l*.l)-rHt (A{l-l) )/(A.')*Otl T) ) 
2* (2. 0*S {'■') *DF L CD« ( 1 ) + FLC ( n -ALF*C;3F ( T ) /PHNO ( I ) ♦ aLF T* ( ( CRF U ♦ 1 ) *CPM 
3(I^l)/PKf'jD(i*l )-, KF (l-l)*CPMa-l) /PRMpU-l) ) / (2..’*IJtLT) )/CPM( I ) ) 
4*2.0*S (M) *F U ) *TmF:TA ( n /uS (M) 

U( 1,3*1) s(ALFT**P.)*CHF (1 ) * (CS ( 1 ♦ 1 ) -2, U*CS { I ) ♦LS ( I-l ) )/ 

1 12.0* (Dtl-1 *^2, ) *'<CfNU(l) ) ♦ (ALF T* (CS (I + l ) -CS( I-l) ) / (4, ^ *OELn ) * 

2 (FLC { 1 ) ♦2.C*S (F' ) t^OFLCUS ( I ) -A1 .F*ChF ( T ) /SCi'^o ( T ) ♦alf T* (CpF (?♦!)/ 
3sCi^iO( 1*1)-CHF (1-1 ) /SCND( l-l) ) / (2,0*PEL N ) *2 . '^*3 (M) *F ( I ) *C^ U ) /DS ( M 

<4 ) 

99 LONTliviUE 

liMS s (HDVIiv. (M) *^CIHU (1)*((2.0*S(M)) **0»3) / ( ALF*HhOM (1 ) *VlSH U ) » 

1 ut (M) * (RaD (•'•t) **tVK ) ) ) * (uELT ♦ .‘3*UFLT*DtLr* ( 1 . - v SCNO ( 1 ) /CRF ( 1 ) ) * 

2 (CPF (2) /SCNO (?) -rRF (1)/SLNG(1) )/uFLT ♦ HUVIN(H)*l (2.' *''(M) )**0.5)* 
3i>CN0 ( 1) / (mLF*RHU.' ( 1) * VlSM ( 1 ) *UE (R) « (RaU (M) **«KK) ) ) ) 

C 

C tVALUATlO'^ DF THF CuEFFICIEMTS 

c 

P(l*l*l)sO,u 

n ( 1 * 1 , 2 ) a 0 , 0 

n(l*l*3)sU»o 

n(l,2,l)=0,u 

H ( 1 ♦ 2 ♦ 2 ) a 0 , 'J 

n(l*2*3)alj,D 

n(i»3*l)=lj,y 

n( 1*3*2) aO.u 

n(l*3*3)a-l.u/ (1 .O^ENS) 

o(l*l*l)aO.U 

0 ( I *2 * 1 ) aThtTArt (•«) 



1^4 I = f ♦ '"N 
IH) 111 J=i*J 

uo 111 ^=i»J 
( 1 ♦ J ,^) =0 • 

L»0 111 L = ltJ 

111 CONTliNUt 

UO lit! J=1»J 
>JQ ll^ ^=1*J 

t3AM{l*J*K)=b(T»J,K)-6n(I,J,K) 
lie LONTlNUfc 

M|)JbAM(I,l,i)=bAt-(I,2,d)«BAH(I»J»3)-BAH(I,3,>')*bAH(I,2,3) 

»«UJHAH ( 1 1 ) =-b^H( 1 ,2, 1) «dAh( 1»3,3) tPArtl I ,2, J) *HAH ( t ,3f 1 ) 

MUjB«H(U3»i)sWAf (I»2»l)*BAH(l»3f?)-HAB(I,2.2)*brtH(l.,3»l) 
Al;JbArt(l,l♦2)s-H^H(l,l,2)«dAH(I»3*3)♦RAM(I,l ,3)«HAH(1 ,1*2) 

MliJBAH ( 1 ,2,2) =BAm ( I* ] »l)OdAH(I»3»3)-BAbU»l,3)*bAH(i»3»U 
«i)JbAM( 1 1 Jfe) s-b.\ri( 1 ♦ 1 » 1) «rtAHn»3,2) ♦RAH( 1 , 1 »2) <»bAH U , 3t l) 

«lJJbAh( 1 » 1 » J) sBA^. ( I» 1 »2) *HAH( 1*2,3) -BAH (1,1 , 3 ) *bAH (1 , 2 , 2 ) 

AUJbAH ( 1 ,2, J) s-R/.H ( 1 , 1 , 1 ) «BAH( I »2,3) ♦PAH{ I , 1 ,3) *«Ah ( 1 , 2 , 1 ) 

AOJbwH ( 1 , J, j) =bAu ( I, 1 , 1 ) *HAH( I ,2,2) -RAb ( I ,2* 1 )*bAH( I ,1,2) 
ub. rbAHsrtan( 1 ,1,1) *BAR ( 1 , 2 , 2 ) *dAh ( 1 , 3 , 3 ) *13 AH ( I ♦ 1 ,2) ( 1 ,2,3) «HAH ( 

li,3,l)*ba'i(l,l ,3) *BA|- (1,3,2) *BAH (1,2,1) -bAH ( I,l,3)*h»»H(l,3,i) «HAH ( 
?1 ,2,2) -bAh( 1 ,2,3) ttBAR ( 1 ,3,2) *BAH (1,1, 1 )-bAH ( I ♦ i ,2 ) *HAh ( I ,2, 1 ) «rAh ( 
3i,3,3) 
ui.) 113 J=i,J 
MO 113 ^sl,3 

B,\^llNV(l.J,^)sAU^BAH(l,J,K)/DtTuaH 

113 continue 

00 llA Jsl,J 

on ii<f ^si,J 

'^'(I,J,K)sU* 

OO 114 L=i,3 

H ( 1 , J ,K ) =H U , J ,^ ) ♦BAH IN V ( I * J,L) «C ( I ,L ,r. ) 

114 l.QNH.iUC 

0(J 121 Jsl,J 

aR ( I , J , 1 ) — 0 • 

OO 121 L=i,3 

( I » J, 1 ) *A0 ( I , J. 1 ) *A ( 1 , J,L) *'J ( 1-] ,L» 1 ) 

121 continue 

UO 12t J=l,3 

o <:■ U ( I , J , 1 ) = u ( T , J . 1 ) - a 3 ( I , J , 1 ) 

122 continue 

UO 123 J=1,J 

0 ( 1 , J , 1 ) = U , 0 
L-0 123 L = l,3 

3(1, J,l) s 0(1, J.l) ♦ bAHiNV ( 1 , 3,1.) ( I ,L» 1 ) 

123 CuNTl ^-Ut 

124 continue 

bOLUTlON For TFF- w VECIOB 
*<(N,l,i) = i,n 

w(N,2,l) * IHET4^ (M+i ) 

'^(N,3,l)~i»n 
r ( N ) s w ( N , i , 1 ) 

1 beta (N) s w(N*2,n 
CS(N) s ’A'(('J,3»1) 

UO 134 11=1, NN 
1 s N - I 
UO 132 Jsl,3 

ri»j (1,3,1) = 1. 

UO 132 Lsl,3 

nW(I,3,i) s hW( 1.3*1) ♦ n ( 1 ,3 »L ) ( I ♦ 1 *L » 1 ) 



non 


132 Continue 

UO U3 J=l,3 

133 CONTiNUt 

»■ ( I ) = w ( 1 , 1 ♦ 1 ) 

tMETAd) = MI»2,1) 

CS ( I ) S w ( 1 »3» 1 ) 

13^ continue 

Ci)NVO s (K(2)-e (i)-h2(2)*F2(l) )/(F(?)-M1) ) 

Lf- (nOPHnT.oI , 0) uo TO Vd 
d-dl'll 9VH #CC(’«V(- 

998 ("OHMaT ( 1 0x*«rCN''G=**E12,3) 

98 COwTiNUt 

IK (Ahs ( (K (2) -F ( li -F2(2) ♦t-2 (1) ) / (F (2)-F ( i) ) ) .LE.O.OOO's) Go 10 150 
00 1 J6 i ~ 1 ♦ N 
K2(I) = F (1) 

t^tETA2(l)sT'^tTA(T) 

Cs2(lJ - CsU) 

138 continue 

UO lJ7 1 = 1*N 

t- ( 1) = 0.5<MF?(I) *81 ( I) ) 

IhETa(I) = .3*( ThETAc ( i) ♦ ThETAI(I)) 

1(1) = ITE*»MFTA(I) 

CS(I) = i'.5*(CS2(l) , ♦ CSl (I) ) 

13? continue 

t\/ALUAT lOi'l uK THf STGEmm f UNC T IOn ,FCL ? 

OO 1A(., 1 s 1»NK 

FF(I) s K (I)/(ALF*(1,0-(1-1 )*uELI) ) 

140 continue 

KLC(l) - 0.i3*(FLr “(M) ♦ FlCw(M + D) 

•■LC(2) s FLC(D ♦ (K(i)/ALK ♦ F (2)/ (AL'‘*tl.''-OELn ) )«D6LT/2*0 
UO l«fl I S j»NK»v 

'FUC(l) * Flc(I- 2) ♦ (FF(I-2) ♦ 4,.,«FF(1-1) ♦ F> ( 1 ) ) «OEl. I / 3 * u 

Vr ) ( . (1 1C 1 A 1 

KLC(I + 1) = FLC(I-l) ♦ (FF(I-l) ♦ 4.0*FF(I) ♦ F K ( i ♦ 1 ) ) *!)El I /3 . 0 

141 continue 

UO 1A2 1 = i*MN 
FCC2(1) = FCC(I) 

OFLCUbd) = (FLC'Ml)-FLCi (I) )/DS(m) 

142 Continue 

IF (nOPKNT.OT .0) GO TO 90 
PRINT V99 

999 FORMAT (3X,ArtFJOTCE»SALVMTION IS NEAR*) 

90 Continue 

go TO 58 
15U CONTINUE 

.V| S fv| 4 ^ 

00 160 I=1*NN 

ETAd) s ALOCMl.'/d.O - yN(I)))/ALF 
fd) s lHETA(D*TTt 

Call mol 1 C.R (H ♦ PE (M) * 2 » T ( 1 ) * Z *£N f ♦RF'O »g ) 

HHOMd)aR!-l(.) 

160 continue 

Factor ~ ( (2,c*S(m) )«*i)*S)/(H!Mot (m)*uf(m)*( (rau(m) )**kk) ) 

T ( 1) a O.U 

UO 161 i a 2»NN 

ok: I Ad) a El Ad) - ETAd-1) 

t(I) S FaCTOK*RH''.E (U) * ( 1 .0/RMUM ( I-l ) ♦ 1 • O/RHOm d ) ) *oetA ( 1 ) /2 , 0 ♦ 
iMI-i) 

161 continue 



I s N 

Uf TA ( 1-1 ) sUtLT/ (-LF*(i,0-YN{I-l) ) ) 
t TA ( i ) =tl A ( 1-1 ) ♦ 'tTA ( 1-1) 

Y ( I) sKhOK (M) «FACrOH* ( 1 , 0/«HOM ( I- 1 ) ♦ 1 , ( I) ) *UETA 

1 1 ) 

UFLSI (1) = '.-,0 

L»tLSI(2) - t-ACTO.v« ( («n()E («^) /hhOM ( ] ) -F ( 1 ) ) ♦ ( HhOt ( N: ) /WHOM ( 2 ) -F ( ? ) 

1 ) / ( 1 ,n-utLT) T/(2.0 *aLF) 

<J0 162 I s -3 » f 2 

UtLSMi) = utL6Tf I-<i) ♦ FACTuH*( (PM-)E(i^')/«HOM(l-2)-F (1-2) )/(1.0 
l-YNU-2)) ♦ A, i)* (KMUE (M) /PMOM l I-l ) -F ( I-i) ) / ( I . v-YN( 1-1 ) ) 

(HhOE {iv ) /r%hOYi { T ) -F I I) ) / ( 1 ,0-YN (T ) ) ) *(jtL T/ ( J . v*aLF ) 
iF ( (l*l ) .(ir.iJfM) ;U TC 162 

i-'tLSid + D = (lELsr(I-l) ♦ FAC(OH«( (SMOt(M)/PMUM(l-l)-F(I-l) )/(1.0- 
ITiM(l-i)) ♦ O* ( 'HOt (H) /KHOM ( 1 ) -F ( I ) ) / ( 1 , 0-YM ( 1 ) J * ( RhOE < M ) / WHOM ( 
21* 1 ) -F ( 1+1 ) ) / ( 1 . -YN( l+l) ) ) ^ucur/ (3,r«ALF) 

162 LO''JT1hUE 

USTAK(M) = UUStKfM) ♦ iJELST(NN) 

IhmOk;(1) = ^,0 

IHMOM (2) = ^ ACTUo* ( (F ( 1 ) -F ( I ) *F ( 1 ) ) ♦ (F (2) -F (2) *F (2) ) / ( I.O-uElD ) * 

1 OELr/(2.0«ALF) 

L**) l63 l=jfi'i|\i»2 

iHMOM(i) = lriM0MfI-2) * KACTOh* ( (F ( 1-2) "F ( 1 - 2 ) *F (1-2) ) / ( 1 *U-YN ( 1-2 

D) ♦ A,0*(F (1-1 )-F ( I-i) *F ( I-l ) )/ (1 ,0-YN U-D ) ♦ (F(n-F(I)«F ( 1 ) )/ 
2(ltU-YN(I) ) ) «r;EL TV (3. 0<»ALF) 

IF ( (I+l) .u| ,NN) oO TO ib3 

lHMOM(I + l) s TFM^'M(I-I) ♦ FACTOH* ( (F ( I-l ) -F ( I-l ) *F ( 1-1 ) ) / (1 .0-YN ( I 
1-D) ♦ A.J4»(F(1).F(I)*F (I) )/(l,(J-YN (I) ) ♦ (F(1 + 1)-F(U1)*F (1 + 1))/ 

2 ( 1 +0-YN ( 1 *1 ) ) ) «upLT/ (3.0*ALF) 

163 tOfYllNUE 

IK (ivOPH|NT.or.C) GU TO 1/Y 
KM IN I 1 + A (F ) 

175 FuWMttT ( IHI * lr,A*'in'1s«» l6*bX»#X (M) a«',Fln»5/) 

K w 1 N I 1 7 b 

176 FOKMm I ( / a , *»1« » 2'- ♦ *YN ( I ) « » 2X * » Y ( I ) « ♦ 3X * *£T a { 1 ) *f *F ( I ) » » bA , t HET 

1 A( 1) *»5A,*C6 ( T) 

UO 1 Hi) 1 si 

HR I NT l/7»ltYN(l) fYd) tETAd) iF(I) ,THFTa(I) ,CS(1) 

177 rOHMAT (5X, 1 lF,Ft ,3»2Fa,b, jEl2,b) 

180 LDNTliMUt 

179 CONTlNUt 

Call MOLltH ( m,FE <m) , 2 , 1 ( 1 ) ,2»£NT ♦WKn*G) 

HHOMd isR'-'O 

CpMd. = Cpid) * CS( l)*(CRS( D-CPI (1) ) 

IF (CS(l)-l.O) 3f'0»301»3ul 
30U ASd) = (CSll) » AMI) /(WHS ♦ CS(l)«(WMl - whS)) 

AT (1) s ) .0 - ASd) 

IF ( XI ( 1 ) • E'J • 0 • i* ) (30 10 301 

ICONS d ) aSPUNiVAL ( 1 A ♦ TtMp ♦ TC ♦ TC I » TC2 » TC3 f T (] ) ) 

VlSS(l)asPLt^YAL (1 A*TEmP,V 1SC»V1»V2»V3, f (i) ) 

ICONI ( 1 ) aTCOiMS ( 1 ) 

VISI (I)=vISp(I) 

wl2(l) = (U.O ♦ ( ( V I S5 ( 1 ) /VISI ( I ) ) • ^) * ( (>»'M1 /wMS) «»0»2i3) ) **?. ) / 

1 (2.82H4* ( ( 1 .U ♦ «MS/WM1) **0.5) ) 

621(1) a ((1,0 ♦ ( ( V I SI ( 1 ) /V ISS ( 1 ) ) **0.5) * ( ( .(MS/WMI ) **0,25) ) <h» 2, ) / 
1 (2.82H4* I d .0 ♦ 'Ml/wMS) **0.b) ) 

VfSMd) a ViSS ( 1 ) / ( 1 . t) ♦ G12 ( 1 ) *XT ( 1) /AS ( 1 ) ) ♦ VlSI(l)/d.o ♦ g21( 

ii)*xsd)/Ai d) ) 

ICONM(l) a I CONS ( 1 ) / ( 1 .0 ♦ 1 . ')6b«G 1 2 d ) *A I ( 1 ) / AS ( 1 ) ) ♦ TCONi(l)/ 
id.O ♦ 1.06b*G21 (1 )*XS(1)/AI (1) ) 

60 TO 302 
30i viSMd) a 


VlSS ( 1 ) 



ICONi'i(l) = ICOKSfl) 

3nc C(,)i>lTlNUt. 

I AU S WHOM (?)-»■ (i) )/ 

1 iOtUT*( (?.0^b(K) >«*U.bn 

LF = TAU/ (o.b«KH. E (t'l) *ue (M) *ue; (M) ) 

hEAT = - Tcui■jM(l)^»TTE*ue('H)*WHOM(l)^^(RAu(M)»*^^)<*ALr»(THElA(<?)- 
1 • -lETA ( 1 n / (ut LT* Mi. (M) ) ) 

U()UT = -HtAl 

tPt*SPLNvML( 14»T-MP,CP»Cr'l,CF';>»CH ^,TF 
iCONtaSHLN\/«L ( lA. TtHPf I C» TC] » iC2»tC3«Tt. (M) ) 

PRWOIN) s VisF (M> «CPE/ rCUNt 

Htc = S>UP r (HHNC ( ■') ) 

m£C S KEC« I I TE- TE (M) ) ♦ TE(M) 

allNO = U!)0 i / It-ihUt-' (M) «Uc (I'l) « ( TPpC-THF T /'W (M) *n t- ) ) 

iAW(M)slpEi 
unUTT (M) sUUUT 
i f NOT (M) sSTI'^U 
bPtCH I (M) =^PE 

IF ((\wPRM.vjT.0) (iO TU 'J<Ji 
PPliMl 30()»TAij,CF 

306 t-oWMAT (//VX»«TAU-o*tli.b» 1 OX ♦ *»CKs* , E 1 ? • S>) 

PKlwr 30 7 f rwtC.iJUUT fbTiMO 

30 7 roPMA r (//yx*»TNtr3tt,Eli.S3iyx»<»auOTs<»»E.l<!»5»RXf *S rN0 = «»El2*i//) 
KPlNl 30P»'^C.LST (^lN) »THi>oOM{NN) *UST / W fM) 

30fc t-OKMAT (VX**UtLSTr:»»Lli.b*dX*<»ThMUMs#,Fli.5» AXi«UsrAHs*,Eli.b//) 

30*^ tOtMTi.^UE 
UO ?«0 

CF (I) s F (l)/(ALF-<»(l,i;-(i-l)«iJEU ) ) 

2«0 CoiMTlNUt. 

M-.C(1) = O.b* (FLr>*' (M) ♦ F LCW ( i')+ i ) ) 

FLC(2) =FLp 11) * (F(i)/ALF >F ( 2 ) / ( ALF* U . r -I /t L I ) ) ) *l»EL T/i . 0 

UO 2oj IsJ,iViN,c 

FLCdi = FLC(T-2) ♦ (FF(i-2) ♦ A.(i«FF(I-l) ♦ F> ( Ji ) > *<’EL T/3 . 0 

IF ( ( l-»i ) .bl .NK) '-0 TO «ibi 

ri,C(i*l) = FLC(l-l) ♦ (F>(I-1) ♦ A.nopFd) ♦ F> ( I + I ) ) «OEl I /3.0 

281 continue 

Uu 2«c' 1 = i»(nJK 

UFLCUS(l) s IFLC ( 1) -FLtl ( I) ) /US (M) 

282 continue 

UO 2^0 IslfN 
M ( I ) s F ( i ) 

IHETAI (1) = THET ' (I) 

csUM = csm 

FlCI ( I) = FCC( I) 

1(1) = iTt^TnETAM) 

2po Continue 

F2(1)*F1 (1) 

F 2 ( 2 ) = F 1 ( 2 ) 

IF (M-MM) 5:*»900.900 
900 CONTINUE 
HP TURN 
EimO 

subroutine spun rRP (n tX»F »R»/»T) — — » 

C * * * * * « « « iHH* ^ 

C* ACCFlPTb AND LtAVES UNALTtREO,,. <» 

C* Ns NUMBtR OF POlMS-i (CURRENTLY OlMENSIONtU FOR 30 pTS) « 

C* X 3 1^DEPENDE^ T VARIABLE ARRAY * 

c* F 3 deperufm variable array « 

C * HOTh X and F must be UImEnSIONEO as TmFY arc. In this RUUTKJF-* 

c* Produces and rf TU h ^ s the parameters r»/«t » 

C* parameters of SPLNtRP characterize the in FERpOLAi INB natural spline* 
c* Function S(X) * 



c 

ulMt^lSIO^ Xt:30)»r(30)»H(jl),Z(30),T(3n)»HH(?'^)»AA(29)»GG(<i^)» 

1 

C 

MJsN-1 
UO 1 1=1 »M 

1 K(I) = (FU*l^-F(I))/lX(Ul)-K(I)) 

I'O 2 Jsi»i\)U 

UU (J) *6» J) -F ( J + 2) )/(;( (J) -X ( ,)♦?)) -K (J) )/ ( A (vJ + 2) - X ( J* i ) j 

UtJMsX ( J*?) - A ( J) 

»*M ( J) s (X { J* i ) - A ( ,|) ) /CUM 

bH(J)=2.'^ 

2 OG ( J) s ( X ( J + 2 ) -A ( , I* 1 ) ) /UUM 
t.uLL yWiSOuv (t.i, A/' ♦bbtGG»uO,Z) 

I'jpaN*! 

X (inIH) sO.O 
X ( 1 ) =0,0 
uo 3 Ksl.N 

■i lu^;)s(z(K♦l)-z(^))/((xl^♦l)-xlK))^^6.^) 

'1 + 2) =K(l)-r (1)<K (X(2)-X(l) )**2) 

C H ) =r? ( “1 ; 

H(NP)=R(MJ♦^A(^F)-X(^))«(Z(N)/'2.04.T(^i)*2•0*(X (HH) -X ('.))) 

WEIUHm 

L.-4L) 

■— » SUbHOUriNE TMlbuLV {^^ALHHA,BETA»^,/iMMA»UtLTAt2) 

C 

DIMENSION ALPHA ( >s>) »Bt TA (29) »UELTA (?9) »GAMMA (29) tZ (3; ) 

C 

iP(2,6T.NU) GO rO 3 
PO 1 1(J*2*NU 

wUUMs (ALPHA (lH))/(bETA(IU-D) 
beta ( ID) sbE I A ( ID) -(JOU.'1*GAMMA ( ID”! ) 

1 UFLTm ( ID) SOC.LTA ( tfj) -UPUM*OELTA ( lU-1 ) 

3 X (N) =UELTA (iMU) /bi-'TA (^U) 

IF (2,6T,IMU> (;U TO 4 

00 2 Ju = 2,fuj 

JI = N-JO 
JI 1=J1+1 

2 X ( J n ) s ( I.) tL I A ( J i ) -GAPma (JI)*Z(J 11 + 1 ))/HETA(JI) 

4 CONTINUE 
HFTUKM 
tiMlJ 

function SpLNVAL rNyXtF «HtZtT* ARC) 

C* Hl^iuHNS THE V-aLUE 'iF 3(AKG) » 

c 

dimension X (jo) *F (30) »K (31 ) *Z (3u) »T (30> 

c 

C MUST first check to see ini WmaT InTfRVAL Apr, LiES 

nN*N*1 

IF ( A ( 1 ) ,LE.AP(i) GU To 4 
SPLN\/ALaF ( 1 ) ♦H (N + 2) * (ApG-X ( 1 ) ) 

GO TO 99 

4 if ( X (Nn) ,ul ,Ak.-;) go TU -3 

SPLNVAL* F' (nN) ♦P /NN) « (AHG-X (N‘1) ) 

GO TO 99 

b U() 1 Is2,nN 

ICNTsI 

1F(X( 1) .gT.AKG) go to 2 



I UJNTlNUt 

d IVAL s ICNT 

Iv/ALlsI^/AL -i 

V AL=Z ( I VALI ) /;?.04T ( IVALl ) A (AH^i*X (IVAL) -<:.0*X ( IVALl ) ) 

V aLsVALA ( ARO-X ( 1 ,/AL) ) ♦H (I VALl ) 

w ALsVAL^ ( Ako-X ( I ■/ ALl) ) (I VALi ) 

SMLN V AL*\' ML 

9V ui,)NTiNUE 

Kt- rUKN 

S>l)kROUT li\it iiOL It.t' (H»F » WOk f »T»^ 

C nCRTsO LUU^ UH uWOPS bASFU ON k AND H 

C iMOktsl L^ON liP MWUPS bASEu ON k AND S 

C iSijkTs? LOOn up PkOkS bAbtO ON k AND T 

C iNiikTsi LOUK OP ku.OkS bAStO ON H ArviC S 

oiMfcNSION rLP(J’»£a) »MZ(33,<i:)) »TT(33,20) ,7T(JJ^<^0) ♦ OAvE ( J J» ?i) ) , 
Ic.mIRU (33»<i0^ »F LkM (6bC ) ♦HZO (6b0) » FTC (660) »ZTO (boO) i G«'”'EO ( bOO ) ♦ 

2tf' FROn (bbO) »I-LP/. f20) »nTBL(33) »ENTRCV (2(^»33) *FLrv (20*33) 
t.NUlVALtlvCE (CLkO.FLk) , (h-iOfH^) » ( T TC . TT ) » ( ZTO , Z I ) , ( fj A^'EO * QA-lt ) , 

1 (?■ NTHOOfPNfHO) 

UA FA kO*CPOH*HC * t OH, G »H, ALE fR rOtCk/?) !*^* *3. ABlbb* H 7 , 3Ab,23.b, 

1 J2.2*b3.3b,2,30 2sa5,33. /U5, ,23866/ 
cOMMON/Pl.Ati/ lCEr,L*IFk 
POMMUN/CHttN/ 1/ 
intALsO 

/: = 1 , il 

oAi1MA=1 *4 

iF(U.E‘J,J3) GL T!) b 

00 2 ks1,2() 

LL=33*(K-1) 
uo 2 Ls1»33 

r LkO (LL + L) sl-LP/ (« ) 
iF(K.cU*l) Ou TO 2 
iVO(LL*L)=H40(L» 

2 U'NTInUE 

1 /s 33 
o/s2;j 

oi) 1 1 = 1 , IZ 
nTbL ( i ) shZ ( 1 * 1 ) 
ug i J=1,JZ 

LnTROV ( J, 1 ) stNilR- ( I , JZ-J* I ) 

1 I-Lkv ( J, I ) =FLk ( I , iZ-J+1) 
b IF (NOP I ,F0, J) (iu FO 4U 
kL*AL061n (p/211b.) ♦ 10. 

IF(NOPT-I) 10,20,30 

10 IF (H.LT, 100, ) (:U TO 100 

CnLL 0 IN F ( H, mZ, Z . Z T , PL,F LPZ , IZ, JZ , T, TT, S.tjjTkU, u , Gamf'a, OaNl ) 

IF (Z-,2222£^30 ) 1 1 , 113* l«J 

11 6 hO= P/( 32 , 2 a 5 J,-, 54 ZaT) 
uo TO 50 

loo IF (H.LT**”^ • ) 00 Tn 12 

laH/Cp 

i^HOsk/ (bARA I ) 

In I b= (CkOkAALyulO (M/HO) -ALObl 0 (k/PO) ) AALE*bOR 
bO TO 50 

12 1FP=1 
iOEAL=l 
bO TO 50 

20 lF(PL.LT.FLkZ(l) ) bO TO 22 
00 21 J=2*J4 

IMkL-FLPZ(o) ) 2 '3*202*21 

21 continue 



2p 2 It-' (S-eNrK0(2, J) ) 20b, 204, 204 

203 iK (S-KNTwo (2, J-1 ) ) 20S>t204,204 

2 04 C/iLL UI^iT ( S* ti'll Wr',H,h y » PL»FLP^» I JZ« T» i T* ;?» / T f 1 » GflHMA,GAMfc. J 
IF (H-,«:222t^J0) n »22,22 

205 fisWTO*CPOH«t KP ( (s-SUfi*ALt«ALU'->lU (P/PO) ) / (CpoW*ALt) ) 

I =H/CP 

hhOsP/ (GAP* I ) 
bu Tu 50 
22 iFPsl 
ll)£AL = l 
wo lU 50 

30 IF ( 1 .LT ’Aly. ) GO TU 301 

Call UINT ( T» IT«n.nZ»PLtKLPZ» J2 , Z » ♦ S ♦ £N ( KU » ,l .OAPMajUAmE) 

IF (H-,2222t-*j,;5) 1 1 ,31,31 
301 1F(T.LT«U«) oO T- 31 
h=CP*T 

KHOsH/ (G«H« I ) 

•o() To 101 

31 1FP=1 
li)£ALsl 
W() TO 50 

40 IF (M.LT? 100» ) t-U TO 401 

uall uint 1 (p»f:m-oo,pl»flpv ,h»htbl* jz, U) 

If (PL-,2222t*30) 1 » 42 *42 

41 lF(PL.Gl.l2»P9t!y7) GO TO 44 
Psl 0.4* (PL-1 U , ) «?1 lb, 

I.M) TU 10 

40l lf(h.LT»0*) wO (■' 42 

ALPS (SOP-S) /ALE*rPOP*ALUblO (H/HU) 

p=po*io»«*'-:lp 
I =H/CP 

HHUsp/ I f 
GO TO 50 

42 IFPsl 

49 lOtALsl 
ui) TO 5(1 

44 10EALS2 

50 PFTUPN 

UflTA (FLPZ ( J) , J = l ,20) /6. 30103, b,6q«q7, 7.0, /,J01U3, 7,0Siby7* R. 
lU, b,30iii3» tt,b9A97* 9,0* 9.30103, q,69b97, 10. 0* lU.30l03» 10.69H 
297, 11 *0, 11 .301 ^3, U,69G97* 12.0* 12 .30 1 o3 * 12 .f>9897/ 
uaTA TTO (L) »L»1 ,204) / 0.» 419,* 834; , 16 16 .* 2339 ,* 1032 ,* 3478 , * 

1J848, *4025, *4175. ,4320. *4468, *4b00. *4880. ,514b. * p 730. ,623b. ,6565. , 
pbRib, *7005, ,715b. ,7270, , /ibO . * 7693 , , 7911 , , H 1 00 . , 8262 . , 84 1 u , , 8568 , , 

36723. . 8900. .9140. .9450., 0., 4l9., 834 , , 1 61 6 , , 2339 , , 3032 . , 35l 7 , , 

43907. . 4130. *430 0. ,446 0 . ,4600 . »4 /bo . , 4940 . ,5260 ., 5760 ., 6300 .» 67n5 . , 

56984.. 7 195., /35o., 7490., 7596., 7927., 81 72,, 0374., 8550., 87 12., 8878., 

69054. . 9250. .9490. .9700., 0.» 419., 834 , , 1 61 6, , 2339. , 3032 , , 3546 , , 
7 J964, *4224, ,4390. ,4563. ,4708. ,4865. ,50^0. ,534b. ,577b. ,6350, ,6795., 
8^092,,7305.,7495.,764b.,7772,,8125,,838b,,859b,,a784.,q950,,9135,, 

99324.. 9520..9750..10030., 0.» 419., 834, , I6lb, ,2339, , 3032 , ,3571 . , 

14035.. 4295 . * 4500 •* 4669*, 4830»*499o.,6lbO«, 5450,, 589 O.,6450«, 6870*, 
2 T2 1 8., 745 0,, 764 0 ., 779 0., 7920. • 8323., 8604., 8827., 90 18., 92 lb,» 94 05., 

39594 . . 9/8v.,, 1003o,, 10360., 0,» 419., 834, , l6lb, ,23.39, , 3032 , *3586, , 

44100. * 44 lu, 14630 . ♦4815. ,4985. ♦6150. ,5340, , 5594 , ,5990 . ,649b. ♦6995. * 
5 / 353 , ♦7620, ♦7840. ,8015. ,8172. ♦8597. ♦8901, ,91 50 . ♦9378, , 957 b. ♦9765., 
699b3.,10l8o.,104?.j.,10/30.,0.^ 419,, 834 ,,1 61 b ., 2339 ., ,3032 3596 , , 

74150. . 4485^,4735. , 4930 . ,5100. ♦5270. , 5480 , ,5724. ,005 J. ,6580. ♦7050. , 
074/0. ♦7760. ♦ 8000 . ♦8190. ♦8334. ♦88 11 . ^9144 , , 94 1 U . ♦ 964d . , 984 b. ♦ 

91 0 06U., 1 0220., 1 0500 ., 1 0740,, ll 040,, 0.,4l;:,, 834., 1616., 2339., 3032./ 
UAfA { TTO (LJ ♦L*2 '5,400) /36 19, ♦ 42 08, , 4555 , ,4^29, ,5040. ,5240. ♦ 5410, , 

15610 . . 5872. . 6195 . .665 0. , / 1 65 .♦ 7596 ., 7908 ., 8 1 50 ., 8355 ,, 8532 ., 9036 , , 



2’5'414. 

3 419, 
43to40, 
59378, 
611880 
7S3£!3. 
R87b0. 
91108 U 
14800* 
?O350. 
311810 
4<*328. 
50118. 
612080 
73604, 
B7810, 
9i220u 

uaTa 
] 01 00, 
29900. 
313960 

4004H, 
59490 , 

61.3960 

79483, 

80692. 

913770 

1C3J9, 

?80oe, 

313280 

4 834 , 
5/066, 
612290 
717220 
RD336, 
910^40 

ue TA 
11616. 
28302. 
314090 
4 419, 
5/080, 
612060 
719450 


♦ 965 / . » 992 7 . f 1 u 1 60 . » 1 038 1' . » 1 0590 , » 1 082 0 . » 1 1 1 0 0 . • 1 14 OU . » 0 • ♦ 

♦ 839, » 1616. ,2339, .3 032. »363o. «4209. ,4650. .4960. »520«i. ♦54 08, ♦ 
,5840, .6066, ,635 0, .6 790. ♦ 7250. .7722, ,8095. ,836 0. , 36oU, .6725. . 
.9 756. . 1006 f' , . 10300. . lOObO. . 10800. ♦ U 040 ,. 11300 ., 1 1 580 . , 

.. O.t 419.. 834. .1616.. 2339, .3032i. 3640. .4300,. 4750. tOObO., 
,5550, .6756. .5990. .6228, .6005. .6900, .7310, .7830. .8215. .8520, . 
,900 ,9630. , 10040, , 10320. . 10660, . 10920 , , 1 1 180. , U420. . 

♦ ,il9^'J«,12?90,,0,. 419.. H 34. .16l8.. 23 3 9. .3 032*. 3654. .4314., 
.5150 •♦ 0445 . .5683 •.0910 • tbiSo • .6394* ,b650 • . 7 045 . , 74R0* . 7956. , 
.8bH0, .6950. .9198, .9864, . 1 0 350 ,. 1 0890 , . 1 1 0 0 0 , . 1 1 290 , ♦ 11550,. 
,. 1210 i-',. 12380,. 12700*. 0.. 419*. 034, ,J616*. 2339,, 3 032, .3654,. 
,48flO, .0295. ,5598. .5080 . .611 0, .6350, ,6628. . 690 0 , , 725u . ♦ 7660, . 
. 855u • ♦ 691 0 . , 92 1 0 . . 9468 . . 1 02 30 • ♦ 1 0 740 . . 1 il5w . . 1 1 480 , ♦ 1 1 79 q , , 

, *12 J7y, ♦12^50,. 12980.. 1 330 0,. 0..4 19,, 834,. 161 6,.23J9,. 3032., 
,4416. ,4945. ,5390. ,0/24. .6015, ,6260. ,8540, ,o808. , 7100. . 7450. , 
.8270, .6690. ,9080. ,9415, ♦ 9 702, . 1 05 1 0 . , 1 1 0 70 , . 1 1 500 . . 1 !« 70 , . 

.. 1201 u .. 12«20 .. 131 10 .. 13400 ,. I3fl00i, 0.. 419., BjH , . l 6 1 6 . / 
(Tro(L> .L *4-1 ♦090/2339,. 3032 •.365‘».,44OO.,4980.,54 75..5839., 
.6450,. 6730., 7002.. 7 3 00,. 7640,. 8000,. 8433,. 6846. ,9243. .9600, , 
. 1 OH. -O.^ 114-0. ♦1 1880,. 12280,, 12630. ♦12970,. 13260,, 13610,. 
,,l43lo., 0., 419., 834, . I6l6. ,2339, ,3032. .3654. ,4489, .5010, , 
,5'’'76, .0 348, ,6650. .6950, .7281, ,7575. ,7910., 8270. ,8656. .9100, ♦ 
.9800,. 1 0l 6^>, » 11200.. 1 1860,. 12420 •♦12820,. 13? 7u.. 13600,. 

. ♦ 14310, . 14A90. . lOObO. .0. .419, , 834 , , 1 6 1 6 , . <;339 . , 3032 . . 3654 , , 
.5020, .0560. .6080, .6474, .6830. . 7l5o, ♦ 748b , , 7doo , , 8 140, . 85()0 , . 

♦ 9300, .970 0, . ICO 70. . 1040 9, . 1 1 5 0 0 , ♦ 1 2260 , ♦ 1 2i'-4n . . 1 330 0 . ♦ 

.♦I4i40,,l4'=:00.. 14880. ♦15270, .15660,. 0.. 419., 

♦3032, ♦3654. . 4497, .5o36.^ 5675, .6J69, ,6580 ,.6970,, 7340. . 7686, . 
.8300, . 8740 . ,91 08. .9025, . 99 3 0 , , 1 03 1 0 , , 1 oOb 3 , ♦ U s4o . . 1 2650 , . 

,. 1 3 bill, »i4'’ 70. ,147 0 0, .15100, . 1 5490^ . 1 5H9u . ♦ 1 62.'f0 , ♦ 0 , . 4 19,. 

♦ 1616. .2339. ♦3032. .3654 .♦40o5,.5n50, .5725 ,. 6266 ,. 8750 . ♦ 7188. . 
.7950, .8300. ,8b70. .9050, .9450, .9850. . 1 0260 . ♦ 1 0660 , . 1 1 052 . , 

, . 131^0, , 1 3 no, . 14000 , . 13010, . 15460, , i593o. ♦ 1 635o , . 1 6800 . . 

,. 0,. 419., 834, . 1616, .2339, ,3052, .3654, .4515, ,50 75, .5750 , . 

.6«48, .722 0 . . 7 750. .8 154, .804 0. .8940, .931 0, .9 720, .lOliu. « 

, . • . 1 1 ^4 0, . 12630, . 1 363 0, .14400, . 1 5o 5 , . 1 ioi q , ♦ 1 61 i ,j , / 

( 1 TO ( L ) . L *5.' 1 . 66 U ) / 1 660 0 . . 1 7 I,! 7 0 . , 1 702o . . 1 7 98 0,."', ,419. .834,, 
.2339, .3032, .3654, . 4425* .51 In, .5774. .6386, .6930. , 7440. .7905, . 
. 8 755,. wiOO., 907 0,. 9990,. 1 04 0 0, .10646,, 11230 ,, 11630. ♦130 10,. 

, ♦ 14940. , 10b4 0». 1624 0, , 1083 0, .1 7320. . 1 78o0 . . 1 83(; 0 . ♦ 1 8b20 , ♦ 0 . . 

, 834., 1616., 2339., 3o 32,. 3654., 4425,, 01 10,. 5 790. ♦ 6440,. 7u30., 
.809 0..859 0. .904 0..9490..99ij5,,in3 70,.l06uu..l I2l0.,llb40.. 

, ♦ I 3O4U , . 14 7 10. . 1 3660. ♦ 1 0480. . 17140. . 1 77 7y , . 1 8 JoO . ♦ 1 89 1 Q . . 

. , l998w,/ 
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— » bUBRUUTlMt InT ( y X»XT » ri » m f 2Z»2T »ML»^L» Y2t YTc:* Y3» YT3, IL» YT4) 

00<JBL£ lNTtKpuLATir.Ni ’jLnROOriNt. IF TL = l» THt StCONO iNOEPtNUtNT 

Variable is nut co^sFani with the fibst, (The mainge of the i>£cuNo 
Is NOT The SahE at each value ok the firs I.) 

U I MENS I ON Zl (NL ) .X r (HLtWL) »YT1 (ML »M ) » Y T2 (Ml. »NL ) « YT3 ( ML »NL) 
UIMENSIOn Y ( (ML.NL) 
iF(ZZ-ZTd) Ftt02,«01»8Ul 

801 lF(lL.F(i.l) 60 Tr 803 
lF(XX,Lr.XT(l»l) »(i0 TO 802 

803 00 800 I=1»NL 
Lsl 

LLsI-1 

1F(2Z ZT(I) )yo2» AOAifiOO 

800 Continue 

802 Y1s,2£;^2e + 3C 

hetukn 

902 HATlPs(ZZ-ZMLL) )/(ZT (L)-ZT(LL) ) 

IF(IL) lUitUlltilU 
80A 00 90A JaltML 
LM* J 
i-LMsJ-1 

iF (XX-XT ( J»L) ) 905»906»90‘+ 

904 continue 
60 TO 802 

905 «ATlOa(XX-XI (LLM.D ) / ( XT ( LM t L ) -X T ( LLM ♦ L ) ) 



n*Yri (LL'^tU.) ♦«ATIU«(YU (LM»Lj-YTl (LLMfL) ) 
T2 =Yt2<ULM,L) ♦NAtI0*(Y|2(LM,L)-YT2(LLM»L) ) 
t 3*YT3 (LLM»L) ♦h^ tIO« ( Y T3 (LM,L) - YT3 (U..M»L) ) 
rtsY T4(LLI^*L) ♦f^urTO*( Y r4(LM,L^-YT4 (LLMtL) ) 

9C6 YlsY I I (LM»L) 

Y2*Y 12 (LMtL) 

^3»Y I 3 (LM»L> 

YtsYr4 (Lm»L) 

HFUJWN 

1111 UO 9i, j Jsi,«iL 
LMsJ 
LLMsJ-1 

iF (XX-X) ( J»U ) 4] , 

<^13 CoNTiNUt 

00 8t>? 

912 YIsYH (LM»LL)-HArlP«(YTl (LMtLD-YTl (LM»U ) 
12aY 12 (LM*Lw) -PAtTP* ( YT2 (LM»LL) -YT2 (l.MfL) ) 
risYr3(LMiLi-)-kA i IP#(YT3 (LM,LL)-Yt3 (LM»L) ) 
Y£*yT 4 (LM*Li-) -WAtIP* { Y T4 (LM,LU -YT4 (LM»L) ) 
HFTURN 

911 KATIOs(AX-Xi (LLM.U )/ (XT (LM»L>-XT (LLPuLM 
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1113 UO 91 P J = l»f'iL 
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9l8 COiMTINUe 
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itssYl 2 (LLM*LL) ♦K/' T IU« ( TT2 (LM*LU - YT? (LLt'^tLL) ) 
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T AsYTCf (LLH»lL) ♦k TIU« ( YTh (LM»LU-YT4 (LLMfLL) > 
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RETURN 

t-i'iO 
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UlMtoSlON Z1 (NU .xr (NLfNL) »YT I (ML,M) 
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1 F (XX-X T(JtL))9l^, 917.918 

918 Continue 
00 TO B(>2 
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916 rat 10= XXX- XI (LLP.L) ) / ( XT (LM.L) -XT (LLM.U ) ) 

YB = Y11 (LLN.t-) ♦RAT1U*(YI1 (LN.L)-YTI (LLMyL) ) 

1114 UO 921 Jsl.ML 
CM = J 
CLMsJ-1 

IF (XX-XT(J.CL) )9i 9,92 0 .921 
921 CONTINUE 
OU TO 0(12 
92U Y4 = Y fl (LM.LU 
00 TO 925 

919 RATIO=(ax-XI (LLM.LL) ) / (XT (LM.LL) -XT (LLf“'»LL) ) 
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RF turn 

End 



DESCRIPTION OF OUTPUT 


The output for the real-gas code includes the flow condition from 
each region of the flow- field model, the geometry of the shock waves and 
the expansion waves , the heat-transfer distribution along the second 
wedge (i.e., the wing leading-edge), the boundary-layer velocity profile 
and temperature distribution profile, and boundary-layer parameters. The 
units for a particular parameter in any region will be the same as the 
free-stream parameter, iinless otherwise noted. The output for the free- 
stream flow includes : 

U1 - free-stream velocity (ft/sec) 

2 

PI - free-stream static pressure (Ibf/ft ) 

T1 - free-stream temperature (°R) 

3 

RHOl - free-stream density (slugs/ft ) 

PTl - free-stream stagnation pressure (lif/ft ) 

TTl - free-stream stagnation temperature (°R) 

GAMMAl - free-stream y 

HI - free-stream enthalpy (Btu/lbm) 

Z1 - free-stream molecular weight ratio 

Output for the flow condition in region "I", where I = 2,3, or 6 includes: 
U1 - the velocity 

PI/Pl - the static pressure ratio 

TI/Tl - the temperature ratio 

RHOI/RHOl - the density ratio 
PI - the static pressure 

TI - the temperature 

RHOI - the density 

PTI - the stagnation pressure 

TTI - the stagnation temperature 
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angle in region 4 and the total change of the Prandtl-Meyer 
angle in region 5 (should be zero so that flow in region 5e 
is parallel to the wing) 


Output of intersection points include: 

INITIAL POINTS ON BOW SHOCK WAVE, i.e., the origin of the coordinate 
system, or the nose. 

WING INTERSECTION POINT, i.e., intersection of the two wedges. 

BOW SHOCK: WING SHOCK INTERSECTION POINT, i.e., the intersection of 

the shock of the first wedge with the shock of the second wedge. 

Next Five Points - the intersection of the centered expansion fan 
with the wing leading edge. 

Last Five Points - the intersections of the reflected waves with the 


inboard shear layer. 



Output listed for the interaction region between the left running and right 
running expansion fan waves includes: 

U4IW - the velocity 
P4IW = - the static pressure 

where "I" = A, B, C, and D. 

The output for the subregions between the inboard shear layer and the shock 
waves between the shock layer and the free stream include: 

THETA - shock wave angle 

X - x-coordinate point of the intersection of the shock waves 

of two adjacent regions (in.) 

y - y-coordinate point of the shock waves intersection (in.) 

V - the velocity 

P - the static pressure 

T - the temperature 

R - the density 

H - the enthalpy (ft^/sec^) 

G - the effective y for the region 

The output for the viscous -region subroutines is tabulated xinder the title 
"OUTPUT DATA" includes: 

M - station number in x-direction 

PE(M) - Static pressure at the edge of the bovmdary layer 
THETAE(M) - temperature ratio (T /T ) at the edge of the boundary layer 
UE(M) - velocity at the edge of the boimdary layer 
URENOE(M) - unit Reynolds number at the edge of the boundary layer 
REXINT(M) - Reynolds number at the edge of the boundary layer inte- 
grated over the distance from the origin 
X(M) - distance along wing leading edge (ft) 

S(M) - transformed x-coordinate 



DS(M) - step-size of S(M) 

BETA(M) - velocity gradient at the edge of the boundary layer 
FLOW(M) - stream function at the wall used to indicate mass injection 
PTE - stagnation pressure at the edge of the boundary layer 

TTE - stagnation temperature at the edge of the boundary layer 

ALF - coordinate transformation parameter 

The next page of output is the initial boundary- layer profile for the first 
x-station. The output that follows is at each of the next (M-1) x-stations 
and Includes : 

M - station number in x-direction 

X(M) - distance along wing leading edge (ft) 

I - station number in y-direction 

YN(I) - transformed H-coordinate (see Ref. ll), n 

Y(I) - physical y- coordinate (ft) 

ETA(I) - transformed y-coordlnate 5 n 

F(I) - velocity profile in the boundary layer 

THETA(I) - temperature distribution profile in the boundary layer 
CS(I) - mass fraction of the stream species 
TAU - skin friction (Ibf/ft^) 

CF - skin friction coefficient 

TREC - recovery temperatxire , ®R 

QDOT - heat transfer rate (Btu / ft^-sec) 

STNO - Stanton number 

DELST - displacement thickness (ft) 

THMOM - momentum thickness (ft) 

DSTAR - displacement thickness with mass injection (see Ref. li) (ft) 
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P3 a 1.<:193?£*02 T3 = 4.07435E^U3 Hho3 a u73O3ae-05 Pt3 b b.79894£4.'j4 

TT3 a 9.b2UVt>E^03 H3 « 1.1740AE^03 23 s l.006«i*t400 

otLlA a 2,bUOOOE^01 theta • 3,o3806F^Ol GAMMA3 » l,2l42nE4ftO 

FLOW CUNOIIIOnS In REGION 

Ob HM/Pl Tb/ri RH06/RHO1 

l.lto2bot*f>4 7,U7399t4-ul 9,232b2E^OO 7.509b3F^oo 

P6 a b«b7792t4oi T6 » 4.52159E*03 HH06 » fl«3oOl9E-o6 PT6 c 1.784UE*04 



V,u75>1jE^03 


H6 « l,444<)Ott03 


Zb s 


U 6 = 

UKlft = 3,0Ul00t*01 T^^ETA 8 3.5>4064E^01 f^AHMAb = l,16?fl5E*nO 

Flu»» CoixUiTlUivS In T»^E EXPANSION FAN 

rtt'-rlUN 4 

VFLutlTie^ 

UA S i,^r2?h3t*04 UH s 1.226 i3E^04 uC a 1.22943E + 04 MU s i.?327]E*04 tiE a l,2.HHy9E*o4 

PWEbSUHES 

PA a i.iDr,bU4u2 PH a 1.0«534E*02 PC a l,0233Ut^02 PO s q.F.4658E*01 pE e P.n9S2lE*(»l 

TFMHtPflTUHf^ S 

lA a 4,VJA/WU*0J tH = 3*9S754E»U3 1C a 3,9447UE*03 TO s 1.90224E + 03 tF = 3,H6nHF*(;3 

OENblTIES 

whoa a i.ftSi^HE-flb PHOH a l,575b3E-o5 MHnc a l.bo250F-0? PHOO a 1.4327oE-05 HhOE = 1.3N64lF.n6 

TNlhALPIES 

riA a I.lb7«7t403 HH a l,l4l7ftE^03 hC a l.l255^E*o3 HD = i,10937E*0J ME a ),0M32nE*o3 

effective Gammas 

oA a l.«ii779c.*0C GH a 1.2?l4OF*Q0 OC a l*225o3E*00 GO a 1.2286/E*u0 r,r a 1.23?3?E + u» 

HtOlOH S 

vELUCITIE*i 


IIA 

= l.<i392bt4-04 

IJH 

* 1.242S3F*V4 

UC a 1.24570t*04 

UD a 

] ,24902E*04 

iiE = 1.2b?2NE*04 





PPEbSURES 




PA 

a d,b76b8E*01 

PH 

a 8,0fl897F^01 

PC a /.63016E401 

PD a 

7.19849E*01 

pF a N,S779?E+nl 





TFMPEkaTUNFs 




1 A 

a 3*8iH3it403 

TH 

a 3. 776R4F*03 

1C a 3,73569E*03 

TO a 

3.69487E^03 

TF a 3,63329F.4fl3 



OENSn lES 


= 1,30J<»4E-Ob RHUb 9 U243t)2E-Ob HHOC « l«l8b78E-05 PHOn c U1327SE-05 RMOF = l.nSUHF-05 


FnTHALPIES 

HA = i,u?7a3t403 HB » 1,060H7F^U3 MC ■ i.OA470E^03 HD a i,02ft53E4O3 mF = 1 . 0 l;»3ftE*03 


EFFFCTiVt Gummas 

Un a i .<: 3 h 0 ot.+()u GB a l,LM 9 hOF*lUi bC s l,? 434 y£ + 00 GD s 1 , 2 ^ 71?£*00 GF = l, 24 « 8 lE *00 


UnUUIF a-S*bB<J 40 fc'-O 4 


PUlr>^Tb up IMtNtST IN IHf FL 0 »^ FIeLO aNU ON THE BUUY (IN INCHES) 


INlTi^L HOInTS On HOW SHOCK waV^ 

K S y S U. 

WiMb iHTtBSFCUOH POINT 

> s u*OuUOi!t+OU Y = B.lOOOOF-nl 

BO'/; bHuCN-wlNb SHOCK iNfEHSFCTlnN POINT 
A s *V.Ji OHt *00 Y s 1 * 4 UHRIF* 0 ^> 

EAHANSln^. Tan TNfEKSECnoiM POINTS 
A s l.Y 31 o 2 t*('l Y s A.SoVBPF^n*’ 

> * |•yJ37lt♦fl^ Y a i#5Bjq6F + nfj 

*■ = l*03oO4t*nl Y s I*bHS45F«00 

A s Y s 1. 61 Oft OF* no 

A = l.UU 16 ut*(>i r a 1 . 62775 F + 00 

. = l*uov?f’t*0l Y = l*966n7F*0i^’ 

' s l*U/ul7L*(Jl r a l*99HhlF*0^» 

A = 1 *0 1 Y a 2 .U 3202 F+OC 

A = l'oH' 33 St^ol Y = 2 » 0 b 7 ? 4 F*ou 

A a l* 09 l 63 t* 0 l V a 2 ,10442 f*oO 



KLUt.. ALOnO [hE WiNli IN yHe wtGIUr- 

VKLUcWltS 

U4AW = l,^2ol3E*n4 IJ4GW = l.23?71t404 u4Ci«< a 1*23‘^2<>E^04 n40v , 1.2457«e*o4 

PHrSbOKFb 

= 1 ,!i«bJ4t*,l2 P4Bw O y,6485^iE^01 P4Cw e 8.57fcot>E*0l P4lJ'.» - 7-630lht*Oi 

ruow (.UNt)HlOHS IN brIOCKEC WEOlON OF EXPANSION fan 
PFtjiUi'J 4FS 

TMFIm s 3.0b272L*ui X = l,ll2b5E+ol Y * 2»42477E*00 

u = l*ilb>7oE + »^4 P 3 N72lS<JF*ol T * 4.b6flrt3F + 03 M 3 «», 3 « 22 l£-nX h s 3 • ^^74 iF-Mi ^ G = l,l5S??t + :io 

PE'JiOiJ SflS 

s J,^8372c*0l A 3 l,l?4d2F4-ol y s 2,52ul2E400 

(I = ].l?F2t^K>04 P 3 7.4i'ft^0F*0l T a 4,nU77^F*03 « s Q,2lHn?E-nX H s 3,^I<}4/F*n7 G 3 l . iS>»r»nF* lO 

Hh'JiON SHS 

TmMA 3 ^.M4J5E + 0i X 3 U13«^4'iF*nl y s 2.62j66E^0n 

U 3 l*i3t)ONt^0* p 3 7.2J278 f^oi r « 4.533blF*03 P s <J.ln333E,nft H = 3.62l7oF^O/ n s 1 . 1 T 

PEi^lON 6CS 

TrifciM a 3«tJ4ftO<*r>ni X 3 1.15374E*^1 Y 3 ^.73b34E*00 

(I = i • I 40 /'^F ♦ (i4 P 3 S»yy277F + ol 7 s 4 , 4b 12'’>E* 0 3 H s H,QX33lE,nb H » l,50b32E*07 G = l.ib^^IPF+oO 

PtGlUN snb 

TPMA s 3.b774nE*0l X = l,17ii77E*Ol Y 3 2.Hbyu3E*00 

U 3 l,i5078E^"4 P s 6.7b4?2F*ol T = 4.38SdSE*03 P s fl.Pl77<»E-oX H s .1 . 1 2rtF ♦(> 7 -i = j.l747nF*.)n 



OUTPUT Data 


^ Pt(M) IrtETrtC-(M) UE(M) UPE^UE<P) KhXlM(M) 

1 1.21937E^n2 4.279J4E-01 1.2l95iE^04 1 .52t>llF^0b U93l5lE*o4 

2 1,21937e^p2 ‘».27934£_(U i.2i951e^04 1.525ije 405 5.794B2fr4()3 

3 1*21937E^p2 4,?79J4E-o 1 U21951E404 1.525llE40b 9.ft57S4E403 

^ 1 .21937E4P2 4,279J4E-ni 1.2l95lE4u4 1 .525i i .352 i.6e4o‘^ 

5 1*21937E4'-2 ‘♦.27934E-01 1.21951E404 1 .525Uf4nb l,73H36E40*» 

b I,M85j4E4n2 ‘».ibbl7E-oi l*226iJE404 i,4ii55f>ub 1 • 74 <*p7e 4 n** 

/ i.0B534E4('2 1.2261JE4C4 1.4l35bE4‘b l*74997E*n4 

« 1.0«534F4^i 2 H,iHbi7E-'U 1.22613E4U4 1 . 4 1 355E4 '>b l,75b66E40** 
9 i.ii8534E402 •♦.i8bi7E-0l 1.226l3E4«.4 l.AlJSSF^'b l.?613SE*0*> 

10 l.o8534E^ot ‘♦•iabl7£-rl l,226lJE404 1 .41355 f*(i^» 1 • 76 7 c^>SE 4 f)4 
U 9.b46b8t*ol ‘♦.o^bbHE-nl 1,23271E404 1.3oo35fc>y5» 1.773l6E*n*^ 

12 9.b465ttE*0l •♦.n9«38E-ol l.23271F-*u4 1.30tt35F4i)b 1*779 p5E*0^ 

13 9.64bb8E4r'i 4.oybb8E-0l 1.2327iE*04 1.30b35F*0b 1.78493F40'^ 

14 9.646b8F40l ^,09bbbE-ol l,2J?71t4'J4 1.3Ub35F40b 1,790H1E*0^ 
lb 9.6463bE + tl ‘►•n9858E-0l 1.23271E404 I * 3i>b35F ♦ Ob 1.7966QE4p4 
lb 8.b7668E*oi ^.oluHiE-nl l.2392oE*fi4 l,2llbOE405 1.8nJi>5E404 
17 H.5766hE 4'»1 4,nlU‘+3£-Pl l.2392bE404 ).?llbOF*o5 l.8n9l7E4n** 
lb 8.5?66bE*nl H.01U43E-01 1.2392PE404 1.2ll50F4(ib Ubl529E4n‘» 

19 8*876bB£40i ‘».0lo43E-ni I .23929E*C'4 1.2ll50F40i> l.bpi^lE^O*^ 

20 8.b766bE4Pi ‘♦•nlO‘^3E-ol 1.2392oE+(;4 1.2ll50F*ob l.b27b3E*0‘» 

21 7.6JoioE + Pi J,92365E-01 l.P457bE4()4 i.l2264E*r)b l.b3420E40‘» 

22 7.6301bE4Ol J,92JobE-')l 1.2457b£404 l,i2264F40b Ub46b2£*o4 

23 7«b3oibE*Ol J*923bbE-ol l*2<^57bE*n4 l*l2264E*o5 l»b47n4E4(,4 

24 7*630l6E*nl J.923b5E-(il 1.24S7bE*04 l,]2264E40b 1.85J46E40^ 

2b 7,63oibE4nl 4,92305£_nl i.24S7be404 l.l22b4F*nb l,b59h8F*o^ 
2b b.b779£E*ol J,(Ubu9t-Ol 1 . 2b22e>£4 l.pr/6?F40b <:• 1 iny4fc'404 

2/ b.b7792E^n -».Hlbn9£-oi 1.2b?2«»£4u4 l.oW67£40b 2.39H7lE*o4 
2b 6#5?792E4i'l J.Hlb09£-i)l 1,2S22 oE 4Q4 i,ol767£*>b 2 • b4b<^8E*0‘* 

29 6.b?792E4ia J , H I b.;9£-o l l,2b22o£404 l,nl 767F4(ib 2.9o425F40*> 

30 6.b^7v2E*^'l J.Hlbo9E-nl 1.2b22bE*04 l.oUb7£4i>b 3.16202E*0‘» 


1 X(M) b(N) L)S(M) 8£TA(M) FLOW(M) 

1 l*2b647F-P2 J,69764 £-o 9 7,3952«E-0Q n, U. 

2 3.79941F-62 1.10929£-nb 7.39528£-09 0, U. 

3 6.3323bE-f<: i.84bb2E-:ib 7.39b28£-09 0. 0. 

A R«8bb29E-fi2 2,SdbJbE-u8 7.21823E-09 0, 0. 

b l.l39br£-fU J.29247E-08 3.b737o£-09 5*n0072f-02 0. 

b U2684HE-(‘1 J.3!!3l0E-nb l,06330£-lo t,67707t>00 U, 

7 i,272biE-ol J.3l373E-0b 1.06330E-10 0, 0. 

b l.276b4£-ol j,32^J7£-u 8 1.0633^»E-lO ft. 1). 

9 1.2fl0o7E-ftl J.33b00E-ftb 1.037b2£-ln ft* 0. 

10 U?b459£«nl J,34bl2E-Mb 1.04l6o£-lr) 1,72438 f*00 0. 

11 1.2b88b£-ol J.355O3E-0** l.O^lbbE-ir, \ €i6\y * nO 0. 

12 l,29335£-;ji J.3bb55E-0b l.o715d£-lo ft, 0. 

13 l.297b4£-oi J.37/26E-0B 1.0715b£-lo ft. 0. 

1*» 1.3o2-^4£-nl J.3b7'^8£-o8 l.ft457i)£-lo 0. 0, 

lb l,306b3£-0l b.39bl8E-0b l,0b37b£-ln 1. 71322^400 0. 

lb l,31lblE-ol ■>.4o9o5£.ob l.oB77uE-lo 1,6B627 f400 0 , 

17 i,316b6E-nl 3.4l993£-ob 1.0877u£-lo 0. 0. 

lb 1,32171E-(!1 3.430OIE-0B 1.08770£-lo 0. Jf 

19 l.32676E-nl J.4^169E-ob 1.06lbiE-lo o, 0« 

20 l.331blE-nl J,45204E-0b 1.0744bE-lft l,6Hb92F*cO 0. 

21 l*337l9E*f'l J*4b3l7E-nb lM1337E-lo i.62ob8F*00 0* 

22 1.34291E-C1 J.474J1E.08 1.11337E-10 ft. 0. 

23 1.34863E-01 J.4hb44E-0b l.llSa/E-lO ft. C. 



24 l*35435E-ri J,49b3aE-f'0 l»0796<iE-10 0 . U. 

2b 1.360' 7E-C1 b.bo7o3E-nti 2.219UE-09 8.22124 f-o2 0. 

26 1.36579E-PA J.94o*»u£-06 4,J336JE-09 4,70b5SF-02 0. 

27 i.619cBE-(U 4,373<6E-(;a 4.3336JE-09 o. 0. 

2a l.a7237E-f^i *^,807l2E-oa 4,333b3E-09 0. ' 0, 

29 2,i2567E-oi P,?40^9E-0« 4,333b3£-n9 0. 0. 


PTts fc.79894t:^n4 


TTEs 9.52096E^03- 


AlFb 8.45070F-01 



1 


F<n 


TMt I A < I ) 


CS(J) 


A n. f,4b723E-0i? l.OOOOOE^OO 

2 2.10499E-02 b.932/0E-02 l.OOOOUE^OO 

3 ^.AieS^E-Pa 1.05012E-01 l.OOOOOE^OO 
6.95 o 29E-0'^ i.21b/3E-0l l.OOOOOE^OO 

b 9.737blE-02 1.3^0to6e-0l l.OOOOOE^OO 
b 1.279/oE-Oi 1.57362E-01 l«t)OOO^E^OO 
/ i.6l63hE-oi i.76433£-01 l.OOOOOE^OO 
d l.9868bE*Ol 1.9foiu2E-0l 1*OOOOOE^OO 
V 2.3957lE-ri 2.1b236E-Ol 1,OOOOOE*00 
l\) 2,84769E-rl c.3bbi>bE-0l l.OOOO^E^UO 

11 3O48B6E-01 d.S6B49E-ol l.OOOOUE^OO 

12 3.904^4E-Ol 2,7ft432E-0l l.OOOOOE^OO 

13 4.b2239E-fil 2.947?oE-ol 1 . 0000 JE *00 
b.20794E-f>i J«1 o7o6E-o 1 l.OOOOOE^OO 

lb 5.9652lE-nl 3.23 i^ 7E-01 1 . f)000vJE*0O 
lb 6.7933bE-ol 3.304U1E-01 1*0000^E*00 
17 7*b7725E-"l 3.3U7u2E-ol 1*OOOUOE»00 
Itt B.b7307E-nl J.22673E-01 l.0000j£*00 
IV V.379o3E-''l J.07022E-01 l.OOOOOE^OO 

20 9.90375E-ri 2.R9473E-01 1.00000E*00 

21 l»OOOOOE^Ou 4.3J273E-01 l*9000uE^OO 


TAUb 2,bbbboE*01 CF= 2,2?l<s8E-02 


T«ECs h.bv^f>fct^03 QUor» l,5e773E*02 STnO* 9.94272F-0^ 


CONV(-»s-3tft7^b/E-Ol 

REJOICE fS^LVA I inr> IS nEAk 
C o Vbs-4* i083(^E-02 
REJOlCE*SAt.VAl ION IS NEAK 
COnVCjs 3.49228E-02 
REJOICE. SAlVAI ION IS NEAH 
CUNV(is-4.7?b0VE-02 
Rf JOlCe.SAuVAl ION IS NEAK 
Co- vbs 2.B57b?E-02 
REJOICE. SALVAIION IS NEAR 
C0NV(iB-1.87b75E-02 
REJOICE. SAuvAl ION IS n£Ah 
Cu-'-vGa 1.0h8«'^E-02 
REJOICE. SAL. '^aUon is NfAH 
CONV<j*-b.23ASbE-U3 
REJOICE. SALVAJ ION is near 
CU VOb 3.438bOE-03 
REJOlCt.SALVAt ION IS NEAR 
C0-'.V0*-1 .VnM4<>E-03 
REJOICE. SALVAI ION IS near 
CONV b« 1.0303VE-03 
REJOICE. SALVAI ION IS nEAk 
CuMVGB-5.bf>3bAE-04 
REJOICE. SAlVAUon is NEAR 
COMVG« 2.95724E-04 




Mb 

2 X (M) s 

,03799 

YN(I) 

til) 

El A(i) 

FID 

10 

*000 

0.00000 

0.00000 

0. 

2 

*050 

•uoool 

•06070 

1.48839E-0? 

3 

.100 

.OU002 

.12468 

3.17431E.02 

4 

.150 

•00003 

.19231 

5.08626E-02 

5 

.200 

.00005 

•264U5 

7.25324E-02 

b 

• 250 

•oonoti 

.34042 

9,71162E.0? 

7 

• 300 

•uool 0 

.42207 

1.24952E-01 

8 

.350 

•U00l4 

.5(j976 

1.56446E-01 

9 

.400 

.00018 

.60448 

1.92072E-01 

10 

• 450 

•00023 

•7o744 

2.32‘»61E-01 

11 

.500 

•00029 

•82022 

2.78408E-01 

12 

.55(1 

.000^*^ 

.94490 

3,3oU2:iE-01 

13 

.600 

• u 0046 

1.08426 

3.90O53E-01 

14 

.650 

.00057 

1.24229 

4.56779E-01 

15 

.700 

.00070 

1.42470 

5.35932E-01 

16 

.750 

• 0008 b 

1.64045 

8.22361E-01 

17 

.800 

•00106 

1.9C450 

7.17104E-01 

18 

,850 

.U0131 

2.24493 

S.16549E.01 

19 

.900 

•UUlbb 

2.72473 

9.12176E-01 

20 

.950 

.00219 

3.54H95 

9.85372E-01 

211 

.000 

•O0?9l 

4.72828 

1.00000E*00 


TMETAtl) CS(I) 

7.45723E-02 l.OOOOOE^OO 
Q-.fc44o7fc-o2 l.OO000E*0U 
1,20383e.oI I.OOOOOe^OO 
1.4563?E-UI l.OOOOOE^OU 
1.71697E-U1 UOUOOOE^OU 
l,*i*849?E-Ol l.OOOOOE^OO 
?»26lfeRE^Ol l.UOOOOE^OU 
2.5501RE-OX l.UOOOOE^OU 
2.S5120E-01 UOOOnOE*00 
3»162 o 3E-U1 1*OOOOOE*00 
3.47542E-UX l.UOOOOE*OU 
3.77907E-O1 UO0000E*0O 
4-0601SE-U1 l.OO000E*0O 
4.29847E-iJl l.OUOOOE»uO 
4.47375E-0X UUO000E*00 
4-b61fi3E-Ol l.OOOnOE^OU 
4.53745E-01 l.OOOOOE^UO 
4.3033OF-UI 1 ,OOOOOf*OU 
4.12447E-01 l.nooOOE*OU 
3*91944 E-j 1 l.OOOnOF^OO 
4.2793^E-01 l.UOOOOE^OO 


TAUa l.lbo97t*Ul 
TRtCa HWul23t*G3 
UELSTa 6.^2703t-04 


CE» 9.06957E-03 

OUOTs u35bb3E*r;2 
THMUMs 3.43955t-(i^ 


STnO= rt,25691F-03 

OSTAk= 6.92703E-n4 


C0NVO=-6.O?fc!31E-02 


rejoic£*salvai lorN is neak 
CO'.V tia 6,77D40E-02 
REJOlCEtSALVAl inn IS near 
CoimVGs* 4 •U473/E-03 
REJOlCEtSAL'^Al ion iS near 
CON iVvjs 2.U1940E-02 
PEJOlC£*SALVAnON IS NEAR 
Co^vOs-l •607‘^3E-03 
REJOIC£*SAlVAI ION is NEAR 
CONVG* 6.90395E-03 
PEJOlCEtSALVAI ION IS NEAR 
CONVG=-5.09363E-04 
REJOlCEtSALVAflON IS NEAR 
Cn. ufi= ,r'.?e9b9E-U3 



Ms A(M)a .13^^>8 


YN( I) 

Til) 

ETA(I) 

F(I) 

THETa ( 1 ) 

C 5 (I) 

10 

• 000 

u^uoooo 

0*00000 

% 

7 . 45723 E -02 

1 . 00000 E 40 O 

^ 

*050 

•00002 

•06070 

2 . 26329 E- 0 ? 

0 . 7970 ftE-o 2 

1 . 0 nnoOE *00 

3 

*100 

•UOOOb 

•12468 

4 . 83174 E -02 

1 . 23079 E-O 1 

l.OOOOOE^OU 

4 

• ISO 

• 00009 

•19231 

7 . 73631 E -02 

1 . 48821 E-uI 

l.OOOOOE^OO 

5 

*200 

•00013 

•26405 

1 . 10129 E -01 

1 . / 4711 E -01 

l.OOOOOE^OU 

b 

.250 

*00019 

*34042 

1 . 46729 E -01 

2 . 00697 E-*oI 

l.OOOOOE^OU 

7 

*300 

*00026 

.42207 

l,a 7 l 9 bE- 0 l 

2 . 269 IU -01 

l.OOooOE +00 

a 

*350 

*00035 

* 5 u 9?6 

2 . 31426 E -01 

2 . 53417 E -01 

l.OOOOOE^OO 

9 

.400 

•00045 

.60448 

2 . 79299 E -01 

2 . 80090 E -01 

l. 0 OO 00 E *00 

11 ) 

*450 

*00057 

• 7 o i ^<>4 

J. 30772 E -01 

3 . 06561 E -01 

l.OOOOOE^OU 

11 

*500 

*00072 


3 . 85931 E .01 

3 . 322 olE-ul 

l.OOOOOE^O^ 

12 

.550 

.ooobv 

.94490 

4 . 44965 E -01 

3 . 56140 E -01 

l.OOOOOE*OU 

13 

*600 

•OOllO 

1 •08428 

5 ,oai 28 E -01 

3 . 77302 E -01 

l.onooOE^oo 

14 

• 650 

.00134 

1.24229 

5 . 75729 E -01 

3 . 94428 E-UJ 

l.OOODOE^OO 

15 

.700 

*00164 

1 . 424/0 

6 , 47692 E -01 

4 . 06170 E -01 

l.onooOE*uo 

lb 

*750 

•U 02 UU 

1 . 6404 b 

/. 24 l 53 t -01 

4 . il 40 «E-Ol 

l.OOOOOE^OO 

17 

*bo 0 

.00244 

1.90450 

8 . 02945 E -01 

4 * 09335 E- 0 i 

l.OOOOOE^OO 

la 

.650 

*00299 

2.24493 

a, 80462 E -01 

3 . 99601 E -01 

1 . 00000 E *00 

1 ^ 

* 9 oO 

• U 0375 

2.72473 

9 . 46685 E -01 

3 . 63732 E -01 

l. 00000 E *00 

20 

• 950 

•O 050 U 

3.54495 

9 . 943 / 7 F .-01 

3 . 69144 E -01 

1 . 00000 E 400 

211 

• 000 

•O 06 B 4 

4.72628 

l. 0 O 0 O 0 E* 0 O 

4 . 01043 E -01 

1 . 00000 E* 0 U 


TAUC 7,3U4(S£ + U0 CFs 7.3U00E-03 


TREC* HW/t>0<!E*O3 
UELST* l•l84^2t-03 


Qi)OT= 5.HgM»^3E*f>l 
THMOMa 7.356 'i^3F-04 


StNOa 4 « 497 fl 0 F- 0 -» 
OSTAWa l,lfl46?F:-'^3 


C0NV&a-3.9o9i;oE-02 

REJOICE, SALVAIION iS NEAR 
CONVGb B*Hft402E-05 



APPENDIX C. -- SOLUTIONS FOR A SINGLE WEDGE CONFIGURATION 


It has been noted that the local heat-transfer rates along the *'wing 
leading-edge" (which were nondimensionalized using a current shuttle de- 
sign parameter) increased significantly with velocity- The velocity de- 
pendence which has been found in the present study of the Type VI shock- 
interference pattern had been found in previous studies of cones (work 
done at the Johnson Space Center) and of wedges (work done at the Univer- 
sity). Using procedures similar to those described previously, theoretical 
solutions have been obtained for a single wedge configuration with a 30° 
deflection angle. The solutions for the three flow conditions of the 
main text : 

Condition 1. - U = 1167 m/sec, P = 2.98 nmHg, T = 53°K 

Condition 2. - U = 4330 m/sec, P = 0.333 mmHg, T = 273°K 

Condition 3. - U = 7610 m/sec, P = 0.0268 mmHg, T = 195°K 

are reviewed briefly in this appendix. 

The heat-transfer rate at a point 0.30m (1.0 ft) from the apex is pre- 
sented in Fig. C.l as a function of the free-stream velocity. The local 
heating has been divided by the theoretical heat-transfer to the stagna- 
tion point of a sphere whose radius is 0.30m. The nondimensionalized 
heat-transfer rates increased by approximately 40% over the velocity range 
considered. The increase in heating with velocity for the single wedge 
is not as great as that observed for the double-wedge configuration which 
produces the Type VI shock-interference pattern. For a given velocity, 
the difference between the perfect-gas and the real-gas solutions is rela- 
tively small. 

Other flow parameters are presented in Table Cl. The pressure ratio 
across the oblique shock wave, the temperature ratio, the density ratio, 
the entropy increase, the shock-wave angle, and the effective gamma of the 


C-1 



shock-layer flow. The relatively high density ratio which occurs when the 
real-gas effects are accounted for results in a thinner shock-layer for the 
real-gas solution. Thus, the effective gamma and the shock-wave angle are 
smaller for the real-gas solution. However, for a given flow condition, 
the pressure ratio for the perfect-gas solution is approximately equal to 
that for the real-gas solution. 

Also presented are the theoretical value of the heat-transfer rate to 
the stagnation point of a sphere whose radius is 0.30m and the entropy in- 
crease across the normal shock for the assumed stagnation point. Note that 
the ratio 



i.e., the increase of entropy across the oblique shock-wave divided by the 
increase of entropy across a normal shock-wave, is approximately 0.5 for all 
three velocities. The constancy of this ratio suggests that the mechanism 
which produces the velocity dependence of the nondimensionalized heat- 
transfer rate is not related to the shock- induced entropy change (which 
relates to the shock-wave strength). This question is relevant because the 
heat transfer for the wedge flow (which passes through an oblique shock) is 
divided by the stagnation-point heat-transfer for the sphere flow-field 
(which passes through a normal shock wave) . 



Table C 1 


Flow parameters for a 30° single-wedge configuration 


Flow 

Condition 

Gas 


y\ 

Pj/Pi 

^^2 ®lV 

9 

s 

Y 

q * 

^t,R=0.3m 

2 

(watts/m ) 

^®t,ns ~ ^iV 

1 

Perf 

29.75 

5.925 

5.021 

— 

39.27 

1.400 

3.47x10^ 

— 

1 

Real 

29.58 

5.915 

5.000 

0.1111 

39.29 

1.397 

3.47x10^ 

0.216 

2 

Perf 

75.51 

13.556 

5.571 

— 

38.01 

1.400 

6.35x10^ 

— 

2 

Real 

70.65 

9.221 

7.510 

0.1666 

35.41 

1.163 

5.96x10^ 

0.338 

3 

Real 

284.21 

20.269 

11.530 

0.3420 

33.25 

1.216 

1.12x10^ 

0.696 



% Perfect-gas solutions 
O Real-gas solutions 



Appendix C : Figure 1. Effect of free-stream velocity on the dimensionless heat- transfer for a wedge 


whose deflection angle is 30°. T = 394°K 



